Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,634,018 T→C 22.7% V46A (GTA→GCA)  yigF → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,634,0180TC22.7% 39.1 / 9.2 22V46A (GTA→GCA) yigFconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (5/12);  new base C (2/3);  total (7/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTG  >  W3110S.gb/3633953‑3634083
                                                                 |                                                                 
aCATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGAc                                                              <  1:1445881/71‑1 (MQ=255)
           cTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAAc                                                    <  1:190650/70‑1 (MQ=255)
                         gCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAAttt                                     <  1:1252842/71‑1 (MQ=255)
                         gCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAAttt                                     <  1:1829963/71‑1 (MQ=255)
                          cAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTc                                    <  1:517735/71‑1 (MQ=255)
                           aaCACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTc                                    <  1:146197/70‑1 (MQ=255)
                                                   tAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:1215936/1‑70 (MQ=255)
                                                   tAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:817563/1‑70 (MQ=255)
                                                   tAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:644635/1‑70 (MQ=255)
                                                   tAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:1310395/1‑70 (MQ=255)
                                                   tAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:1119163/1‑70 (MQ=255)
                                                   tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:1860471/1‑70 (MQ=255)
                                                   tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:273609/1‑70 (MQ=255)
                                                      aGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATAtatt        <  1:162739/71‑1 (MQ=255)
                                                      aGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATAtatt        <  1:150459/71‑1 (MQ=255)
                                                      aGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATAtatt        <  1:466155/71‑1 (MQ=255)
                                                          ttAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGa              <  1:1459711/61‑1 (MQ=255)
                                                          ttAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGa              <  1:771189/61‑1 (MQ=255)
                                                          ttAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGa              <  1:931995/61‑1 (MQ=255)
                                                          ttAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGa              <  1:970203/61‑1 (MQ=255)
                                                          ttAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCT‑‑TGTTCCTGGATATATTTATTTg  <  1:1689612/71‑1 (MQ=255)
                                                          ttAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTCCTTTCTTCTGTTCCTGGa              <  1:1780600/61‑1 (MQ=255)
                                                                 |                                                                 
ACATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTG  >  W3110S.gb/3633953‑3634083

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