Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,313,545 C→T 45.8% intergenic (+23/+17) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD complex/predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,313,5450CT45.8% 2.6 / 24.2 25intergenic (+23/+17)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD complex/predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (4/9);  new base T (7/4);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTC  >  W3110S.gb/1313475‑1313594
                                                                      |                                                 
gtggtgAATATCCTGTTTAAAATTAACGTCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                   >  1:2172711/1‑71 (MQ=255)
gtggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                   >  1:2388109/1‑71 (MQ=255)
     gAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAg                                              >  1:2365984/1‑71 (MQ=255)
        tatCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCt                                           <  1:3599729/71‑1 (MQ=255)
        tatCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGGCAAAAGCCTCCGACCGGAGGCt                                           <  1:1312857/71‑1 (MQ=255)
               tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                        <  1:92652/67‑1 (MQ=255)
               tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACt                                    <  1:592235/71‑1 (MQ=255)
                  aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTAt                                  <  1:3733294/70‑1 (MQ=255)
                    aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAcc                                                  <  1:1753052/52‑1 (MQ=255)
                       tAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGCCCGGAGGCttt                                         >  1:3838885/1‑58 (MQ=255)
                         aCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAcc                                                  <  1:71880/47‑1 (MQ=255)
                             caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGt                      <  1:699770/71‑1 (MQ=255)
                              ttaattaacggcactcctcagccAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                        <  1:2514529‑M2/29‑1 (MQ=255)
                              accaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGt                      <  1:397665/70‑1 (MQ=255)
                               taattaacggcactcctcagccAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGa                                      >  1:1351664‑M2/23‑53 (MQ=255)
                                aattaacggcactcctcagccAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGActaca                                 >  1:847068‑M2/22‑55 (MQ=255)
                                 attaacggcactcctcagccAAGTCAAAAGCCTCCGGGCGGAGGCttt                                         <  1:3574222‑M2/28‑1 (MQ=255)
                                  ttaacggcactcctcagccAAGTCAAAAGCCTCCGGTCGGa                                               >  1:1586356‑M2/20‑41 (MQ=255)
                                  ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgc                 >  1:2247100/1‑71 (MQ=255)
                                    gAAATTCAGTAAGCAGAAAGTCAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgcgc               >  1:1941081/1‑70 (MQ=255)
                                    aacggcactcctcagccAAGTCAAAAGCCTCCGGTCGGa                                               >  1:3681104‑M2/18‑39 (MQ=255)
                                    aacggcactcctcagccAAGTCAAAAGCCTCCGGTCGGa                                               >  1:1675932‑M2/18‑39 (MQ=255)
                                    aacggcactcctcagccAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGa                                      >  1:586721‑M2/18‑48 (MQ=255)
                                          actcctcagccAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTa                                   <  1:3440577‑M2/34‑1 (MQ=255)
                                              ctcagccAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGActaca                                 >  1:745026‑M2/8‑41 (MQ=255)
                                              ctcagccAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTa                                   <  1:3661786‑M2/34‑1 (MQ=255)
                                              ctcagccAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTa                                   <  1:3554717‑M2/34‑1 (MQ=255)
                                                 cAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTc  <  1:1031040/71‑1 (MQ=255)
                                                          aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt        <  1:691256/56‑1 (MQ=255)
                                                                      |                                                 
GTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTC  >  W3110S.gb/1313475‑1313594

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