Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,313,823 T→C 36.0% E349G (GAA→GGA)  yojL ← predicted thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,313,8230TC36.0% 21.2 / 21.5 25E349G (GAA→GGA) yojLpredicted thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base C (2/7);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCG  >  W3110S.gb/2313756‑2313876
                                                                   |                                                     
aaGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGc                                                    >  1:1335273/1‑71 (MQ=255)
aaGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGc                                                    >  1:403800/1‑71 (MQ=255)
      aGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGAc                                                >  1:115510/1‑69 (MQ=255)
                    ccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAAc                                <  1:545190/71‑1 (MQ=255)
                    ccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAAc                                <  1:50208/71‑1 (MQ=255)
                    ccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAAc                                <  1:345775/71‑1 (MQ=255)
                    ccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAAc                                <  1:1694043/71‑1 (MQ=255)
                    ccGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAAc                                <  1:1000452/71‑1 (MQ=255)
                           ccATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                          <  1:1853719/70‑1 (MQ=255)
                           ccATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                          <  1:201625/70‑1 (MQ=255)
                           ccATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                          <  1:1918629/70‑1 (MQ=255)
                           ccATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                          <  1:1596134/70‑1 (MQ=255)
                                aCGCAAAAACCAACGATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                    >  1:492487/1‑71 (MQ=255)
                                aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                    >  1:1785343/1‑71 (MQ=255)
                                aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                    >  1:1778620/1‑71 (MQ=255)
                                aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                    >  1:1445709/1‑71 (MQ=255)
                                aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                    >  1:1446200/1‑71 (MQ=255)
                                aCGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATc                    >  1:1181717/1‑71 (MQ=255)
                                 cGCAAAAACCAACAATCTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATcc                   <  1:153674/71‑1 (MQ=255)
                                 cGCAAAAACCAACAATCTTGCGCTTTAATTTCTTCCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATcc                   <  1:181567/71‑1 (MQ=255)
                                      aaaCCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGtt              <  1:776012/71‑1 (MQ=255)
                                      aaaCCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGtt              <  1:778167/71‑1 (MQ=255)
                                                 cTTGCGCTTTAATTTTTTCCGCTGACAAGGAAGCTTTTAAACTgtggt                          <  1:961845/48‑1 (MQ=255)
                                                    gcgcTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCg  <  1:1944057/69‑1 (MQ=255)
                                                    gcgcTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCg  <  1:1961842/69‑1 (MQ=255)
                                                                   |                                                     
AAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAATCTTGCGCTTTAATTTTTTTCGCTGACAAGGAAGCTTTTAAACTGTGGTGACATCCAGGTTTTAAAGCTATCG  >  W3110S.gb/2313756‑2313876

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: