Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,032,491 C→T 63.3% intergenic (‑357/‑70) serT ← / → hyaA tRNA‑Ser/hydrogenase 1, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,032,4910CT63.3% 21.3 / 28.0 30intergenic (‑357/‑70)serT/hyaAtRNA‑Ser/hydrogenase 1, small subunit
Reads supporting (aligned to +/- strand):  ref base C (7/4);  new base T (13/6);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GTGCGCCAGTGCAGAAGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGCAAGGAG  >  W3110S.gb/1032423‑1032546
                                                                    |                                                       
ttGCGCCAGTGCAGAAGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTACAg                                                       >  1:1635761/2‑71 (MQ=255)
gTGCGCCAGTGCAGAAGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTACAg                                                       >  1:2996735/1‑71 (MQ=255)
gTGCGCCAGTGCAGAAGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTACAg                                                       >  1:2985612/1‑71 (MQ=255)
               aGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGcc                                         >  1:86047/1‑70 (MQ=255)
               aGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGcc                                         >  1:3305825/1‑70 (MQ=255)
               aGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGcc                                         >  1:2502990/1‑70 (MQ=255)
               aGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGAACGTCGCTTGcc                                         >  1:1606507/1‑70 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:1405305/1‑48 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:903297/1‑48 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:705141/1‑48 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:416452/1‑48 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:3354631/1‑48 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:1229673/1‑48 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:1483248/1‑48 (MQ=255)
                          tCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTATAGGGa                                                    >  1:1982183/1‑48 (MQ=255)
                          tCGATTTCAGCATCTCACGTGAAGCGATGGTTTGCCTAGCTATAGGGa                                                    >  1:1614000/1‑48 (MQ=37)
                                          aCGTGAAGCGATGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTc                     >  1:2669692/1‑63 (MQ=255)
                                          aCGTGAAGCGATGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTc                     >  1:2768312/1‑63 (MQ=255)
                                          aCGTGAAGCGATGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTc                     >  1:1777630/1‑63 (MQ=255)
                                          aCGTGAAGCGATGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTc                     >  1:216780/1‑63 (MQ=255)
                                               aaGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGc        <  1:1807474/71‑1 (MQ=255)
                                               aaGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGc        <  1:3187878/71‑1 (MQ=255)
                                               aaGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGc        <  1:1923304/71‑1 (MQ=255)
                                               aaGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGc        <  1:2791160/71‑1 (MQ=255)
                                               aaGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGc        <  1:2454641/71‑1 (MQ=255)
                                                aGCGATGGTTTGCCTTGCTATAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGc        <  1:2467005/70‑1 (MQ=255)
                                                     tGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGCAaggag  <  1:3317984/71‑1 (MQ=255)
                                                     tGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGCAaggag  <  1:1496557/71‑1 (MQ=255)
                                                     tGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGCAaggag  <  1:1426864/71‑1 (MQ=255)
                                                            ccTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGGCGCAaggag  <  1:2303121/64‑1 (MQ=255)
                                                                    |                                                       
GTGCGCCAGTGCAGAAGGATGAGCTTTCGTTTTCAGCATCTCACGTGAAGCGATGGTTTGCCTTGCTACAGGGACGTCGCTTGCCGACCATAAGCGCCCGGTGTCCTGCCGGTGTCGCAAGGAG  >  W3110S.gb/1032423‑1032546

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: