Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,413,523 G→A 15.2% S47F (TCC→TTC)  ydaO ← predicted C32 tRNA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,413,5230GA15.2% 119.0 / 11.8 46S47F (TCC→TTC) ydaOpredicted C32 tRNA thiolase
Reads supporting (aligned to +/- strand):  ref base G (24/15);  new base A (3/4);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACG  >  AP009048/1413455‑1413588
                                                                    |                                                                 
aaTTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGgag                                                                 <  1:1590973/71‑1 (MQ=255)
    gattgGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGAAGAgg                                                              >  1:558701/1‑70 (MQ=255)
      ttgGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCa                                                            >  1:148710/1‑70 (MQ=255)
        ggCGCGCTTTGCTGCAAATTGCGCAGCATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAAc                                                         >  1:3230226/1‑71 (MQ=255)
         gcgcgcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAAcc                                                        <  1:3544605/71‑1 (MQ=255)
         gcgcgcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAAc                                                         >  1:3065043/1‑70 (MQ=255)
          cgcgcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCa                                                       <  1:385461/71‑1 (MQ=255)
             gcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCatg                                                     <  1:3436608/70‑1 (MQ=255)
              cTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCatgat                                                   >  1:1891968/1‑71 (MQ=255)
                ttGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGc                                                 >  1:2982630/1‑71 (MQ=255)
                 tgctgcAAATTGCGCAGAATCTCAAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACAATGATGc                                                 >  1:1413717/1‑70 (MQ=255)
                  gctgcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCat                                                      >  1:1346353/1‑64 (MQ=255)
                      cAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCg                                            >  1:3461837/1‑70 (MQ=255)
                       aaaTTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGga                                                                  >  1:2614015/1‑47 (MQ=255)
                       aaaTTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAACCCATGATGCGATCGcc                                          >  1:2224166/1‑71 (MQ=255)
                        aaTTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCct                                         <  1:4097311/71‑1 (MQ=255)
                         aTTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCctt                                        <  1:2820600/71‑1 (MQ=255)
                          ttGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGcc                                          >  1:3100232/1‑68 (MQ=255)
                          ttGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCcttc                                       <  1:2869855/71‑1 (MQ=255)
                            gcgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCcttct                                      <  1:4517296/70‑1 (MQ=255)
                            gcgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCaaa                                                          >  1:3669681/1‑50 (MQ=255)
                                 gAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGAGATCGCCTTCTTCAAt                                 >  1:1650448/1‑70 (MQ=255)
                                  aaTCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCttcttc                                    <  1:88456/66‑1 (MQ=255)
                                   aTCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCa                               >  1:2180245/1‑70 (MQ=255)
                                      tcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCata                             >  1:4611225/1‑69 (MQ=255)
                                      tcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATAtt                           >  1:2287359/1‑71 (MQ=255)
                                      tcCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCg                                                <  1:3471228/50‑1 (MQ=255)
                                        cAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTg                          >  1:150542/1‑70 (MQ=255)
                                           cATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGt                      <  1:1316951/71‑1 (MQ=255)
                                             tGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGt                      >  1:1375779/1‑69 (MQ=255)
                                               gTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGc                  <  1:3629170/71‑1 (MQ=255)
                                                 ataGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCaa                >  1:348866/1‑71 (MQ=255)
                                                  taGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAAt               <  1:4137564/71‑1 (MQ=255)
                                                    gCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATg              <  1:3601729/70‑1 (MQ=255)
                                                     cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCata                             >  1:3402552/1‑54 (MQ=255)
                                                     cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCaa                <  1:4338655/67‑1 (MQ=255)
                                                        tCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCtt          <  1:4206126/70‑1 (MQ=255)
                                                        tCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCt           <  1:586502/69‑1 (MQ=255)
                                                        tCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTc         <  1:762504/71‑1 (MQ=255)
                                                        tCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCt           <  1:293102/69‑1 (MQ=255)
                                                          tttACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATg              >  1:3915365/1‑64 (MQ=255)
                                                          tttACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCg        >  1:2973483/1‑70 (MQ=255)
                                                          tttACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATATGATTGAAgaa                     >  1:1745376/1‑55 (MQ=37)
                                                            tACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGcc      >  1:2820244/1‑70 (MQ=255)
                                                            tACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGcc      >  1:1171672/1‑70 (MQ=255)
                                                             aCCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGccc     <  1:326326/70‑1 (MQ=255)
                                                             aCCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCa    >  1:1763566/1‑71 (MQ=255)
                                                               ccccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACg  >  1:4499022/1‑71 (MQ=255)
                                                                   ggAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAg                   >  1:3770679/1‑50 (MQ=255)
                                                                    |                                                                 
AATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACG  >  AP009048/1413455‑1413588

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: