Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,356	C→A	31.3%	L37L (CTG→CTT)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,356	0	C	A	31.3%	52.0 / 3.2	30	L37L (CTG→CTT)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base C (8/12);  new base A (2/8);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCGTTGCTTTATAG‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCC  >  NC_000913/954126‑954548                                                                                                                                                                                                                                          |                                                                                                                                                                                                   cAGTTTTTCGCCAACTGACCCCAGGTTATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGa                                                                                                                                                                                             <  1:129106/238‑1 (MQ=255) cAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGa                                                                                                                                                                                             <  1:520653/238‑1 (MQ=255)         cGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAAGGCGATTGAGAGGAAAACACCG‑GC‑GGGAATCTCCAGATAGaaaaagataa                                                                                                                                                                                    >  2:603136/1‑232 (MQ=255)                     ctcggTGTTCCGCCCCCTCCCCTTAGAAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCTCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATTCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACGCCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTcc                                                                                                                                                                       <  2:104260/236‑2 (MQ=255)                     ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCg                                                                                                                                                                       <  1:689889/239‑1 (MQ=255)                     ccAGGGGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCg                                                                                                                                                                       <  1:814709/239‑1 (MQ=255)                     ccAAGTTATGCCCCCTCTCGCCTTAGCAACAACAAGTAACACGGTGGAAGTAAAGAGATCCGCTCCGCAGACAACACAAATAATCAGCCCCATAGAGAAGCAAATTCCGACAACCAGTTTTGCCATGCCGAAGGGCAGTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACCCCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCg                                                                                                                                                                       <  2:151319/239‑1 (MQ=255)                         gggttGCGCCCCCTCCCCTAAGAAAAAAAAATCAACACGGGGGAAATAAAGAGATCGGCTCCGCAGAAAACCCAAAGAACTAGCCCCAAGAAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GGCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACGCCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCGTTGc                                                                                                                                                                   <  2:701664/235‑1 (MQ=255)                                                  caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGATTGCGGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GTCAATTGCCCGAAAGACAGCCTCAAGC‑GGAGGTTTCACTGTTTTTTAGAACACAC‑CCTCTTCGGc                                                                                                                                          >  1:492339/1‑238 (MQ=255)                                                     caaTCAACGCGGGGGAAGAAAAGAGATCGGCCCCGAAGACAACAAAAAAAATCAGCCCAAGAGAAAAGCAAATGCCGCCAACCAGTTTTGCCATGCAGAAGGGCATTGT‑GCCTGTGCCAGGGGTGGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCGTTGTTTTATAG‑ACACCCGCCTCTTCGGCCACt                                                                                                                                      <  1:603136/239‑1 (MQ=255)                                                     caaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCg                                                                                                                                                                       <  1:168199/207‑1 (MQ=255)                                                        tCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GCGGGAAATCCCAAGTAAGAAGGCTTTAAGG‑GGATG‑TTCGTTGGTTTATAG‑GCACCCCCCCCTTGCGCCACTTTg                                                                                                                                   >  1:354358/1‑239 (MQ=255)                                                                gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGCGATGAAAACACCG‑GCGGGAAATCCCCAAATAAAAAGTTTTAAGG‑GGGGATTTATTGGTTTTATAT‑AAACCCGCCCTTTAGGCCCATTTGGcatttc                                                                                                                            >  2:104355/1‑234 (MQ=255)                                                                gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATAGCCAGATAGAAAGGCTTAAAC‑GGGTGTTTTGTTGGTTTTTAG‑ACACCCCCCCCTTTGGCCCCTTTTGCCATTgc                                                                                                                           >  1:66297/1‑239 (MQ=255)                                                                    aaGTAAAGAGGCCGGCCCCCAAAACACTAAAAAAAACCACCCCAGGGGAGAACAAAATGCCGCCAACAGGTTTTCCCAGGCCGAAGGGGATGGG‑GCTGGTGCCGGGGGTTGTGGTAAAATGGAAAACGAATGGGATTGAGATAAAACCCCGG‑GC‑GGTAATCGCAAGATAGAAAGTTTAAAGC‑GGAGTTTTCGTTGCTTTAAAA‑ACCCCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                       <  2:57672/239‑1 (MQ=255)                                                                       tAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGGGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GA‑GGTAAGAGCACGATAGAAAGTCTTTAGA‑CGATGGTTTGTTGCTTTATTGAACACCC‑CCTCTTCGGCCACCTTGGGCCGTGTAGCAGGa                                                                                                                    >  2:189378/1‑239 (MQ=255)                                                                              aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTGATCGCCCATGCGAAAGGCTTTAAG‑GGGTGTTTGGTTGCTTTATAG‑AACCCCACCCCCTCGGGCACTTTGGGCCTTTCAGCAGGAGGGaaaa                                                                                                             >  2:118397/1‑239 (MQ=255)                                                                              aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATTGAAAAGCTTAAGC‑GGATGTTTCGTTGCTTTATAG‑AAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                             >  1:329878/1‑239 (MQ=255)                                                                               tCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAAAACCG‑GC‑GGTAATCGCCCGAAAGAAAGACTTAAGCGGGAGGTTTCG‑TGCGTTTAAG‑ACACCCGCCCCTTCGGCCCCTTTGGGCCTTGTACCAGGAAAGaaaaa                                                                                                            >  2:380114/1‑238 (MQ=255)                                                                                             aacaCAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGa                                                                                                                                                                                             >  1:116476/1‑146 (MQ=255)                                                                                             aacaCAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGa                                                                                                                                                                                             <  2:116476/146‑1 (MQ=255)                                                                                                                                      ccggTTTGGCCAGGCCGAAGGG‑ATTGTGGCCGGTCCCAGGGGTTGCTGTAATATAGAAGACAAAGGCGTTTGAGATGAAAACCCCG‑GG‑GGTAATCGCAAGAGAAAAAGTCTTAAGC‑GGATGTTTCTTTGCTTTATAG‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgacga                                                     <  1:853610/236‑1 (MQ=255)                                                                                                                                                cATGCGGAGGGGCATTGT‑GCTTTGTCCGGTGTTTGGTATAATATAGAAGAAGAATGCGATTGAGATGAAAACCCCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAGGC‑GTATTTTTCGTTGCTTTAAAG‑ACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                           <  1:93523/239‑1 (MQ=255)                                                                                                                                                    gcGAGGGGCTTTGT‑GCCTGTGCCAGGGTTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCACAAGATGAAAATTCTTAAGC‑GATTTTTTCGTTGTTTTATAA‑CCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                        <  1:118397/237‑1 (MQ=255)                                                                                                                                                                     ctgtgccAGTGTTTGCTGGTAAATAGAAGACAAAGGCGTTTGAGATAAAACCCCCG‑GC‑GGTAACCGCAAGATAGAAAGTCTAAGGC‑GGATGTTTCGTTGCTTTAGAG‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                       <  1:380114/239‑1 (MQ=255)                                                                                                                                                                     ctgggcaAGGGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAACCCCCG‑GC‑GGTAACCGCAAGATAGAAAGTTTAAAGC‑GGATGTTTCGTTGCTTAAAAG‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                       <  1:104355/239‑1 (MQ=255)                                                                                                                                                                           cAGGGGTTGGTGTAATAAAAAAGGCGATTGCGATGGAAATGAAAACCCCG‑GC‑GGTAACCGCAAGATAGAAAGTCTAAAGG‑GGTTGTTGCGTTGCTTTAAAG‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                  <  1:116694/238‑1 (MQ=255)                                                                                                                                                                                    gtgtgATAAAGAAGAGAAAGGCGATTGAGATGAAAACCCCG‑GG‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGTTTTTCCGTTGCTTTAAAG‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCatat        <  1:741800/238‑1 (MQ=255)                                                                                                                                                                                          taaagaagaCGAATGCGATTGAGAT‑AAAACCCCGTGG‑GTTAACCGCAAGATAGAAAGTCTTAAGG‑GGTTGTTTCGTGTCTTTATAA‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc  <  1:405251/236‑1 (MQ=255)                                                                                                                                                                                          taaaaaagaCAAAGGCGATTGATAGAAAAACCCCG‑GC‑GGTAACGGCAAGATAGAAATTCTTAAGC‑GGTTGTTTCTTTGCTTTAAAG‑ACCCCCGCCTCTTCGGCCACTTTGGCCATTACAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc  <  1:586636/234‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                   CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCG‑GC‑GGTAATCGCCAGATAGAAAGTCTTAAGC‑GGATGTTTCGTTGCTTTATAG‑ACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCC  >  NC_000913/954126‑954548

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑