Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,051,238 T→C 74.1% intergenic (+286/+147) insB → / ← cspH IS1 transposase InsAB'/stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,051,2380TC74.1% 45.4 / 13.9 27intergenic (+286/+147)insB/cspHIS1 transposase InsAB'/stress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (7/13);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTTGTGGAA  >  W3110S.gb/1051175‑1051302
                                                               |                                                                
ttGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGc                                                           <  1:1344787/71‑1 (MQ=255)
      cGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTtat                                                     <  1:543226/71‑1 (MQ=255)
      cGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTtat                                                     <  1:1711098/71‑1 (MQ=255)
      cGAAAGCTATTAATTATTGTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTtat                                                     <  1:1046732/70‑1 (MQ=255)
         aaGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTTAAATGGCATTTATAAc                                                  <  1:1417191/71‑1 (MQ=255)
              attaattaTTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTTTCAAATGGCATTTATAACTAt                                               >  1:2157103/1‑69 (MQ=255)
                 aattaTTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGt                                          >  1:2606504/1‑71 (MQ=255)
                   ttattGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGc                                        <  1:714523/71‑1 (MQ=255)
                   ttattGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGc                                        <  1:2277848/71‑1 (MQ=255)
                               atatAGTTTCAAGTGGTACTATAACTGGCTGTCAAAt                                                              <  1:1947693/37‑1 (MQ=255)
                               atatAGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGc                            >  1:2385202/1‑71 (MQ=255)
                                tataGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAAc                                                  <  1:443777/48‑1 (MQ=255)
                                 ataGTTTCAAGTGGTACTATAACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCt                           >  1:591418/1‑70 (MQ=255)
                                     tttCAAGTGGTACTATAACTGGCTGTCAAATGGCAtt                                                        <  1:900021/37‑1 (MQ=255)
                                       tCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATg                     >  1:2614056/1‑70 (MQ=255)
                                       tCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATg                     >  1:1180582/1‑70 (MQ=255)
                                       tCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATg                     >  1:1057123/1‑70 (MQ=255)
                                         aaGTGGTACTATAACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTa                         <  1:1302494/64‑1 (MQ=255)
                                         aaGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGAt                  <  1:2461396/71‑1 (MQ=255)
                                              gTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGt              <  1:2399605/70‑1 (MQ=255)
                                                 cTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTa                         <  1:2704600/56‑1 (MQ=255)
                                                   ataACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGata                                 >  1:2224388/1‑46 (MQ=255)
                                                    taACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTaaa                       >  1:2166804/1‑55 (MQ=255)
                                                    taACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTaaa                       >  1:1923466/1‑55 (MQ=255)
                                                          gCTGTTAAATGGCATTTATAACTATTAGGTGCAgaga                                   >  1:1333243/1‑37 (MQ=255)
                                                          gCTGTTAAATGGCATTTATAACTATTAGGTGCAgaga                                   >  1:623665/1‑37 (MQ=255)
                                                          gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTTGTGGaa  <  1:310844/70‑1 (MQ=255)
                                                               |                                                                
TTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTTGTGGAA  >  W3110S.gb/1051175‑1051302

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: