Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,395	G→A	36.9%	V24V (GTC→GTT)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,395	0	G	A	36.9%	12.1 / 4.9	28	V24V (GTC→GTT)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base G (4/12);  new base A (4/6);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCC‑AGAT‑AG‑AAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGT  >  NC_000913/954185‑954595                                                                                                                                                                                                                              |                                                                                                                                                                                                            cacGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGGAAACACCGGCGGTAATTGCC‑AGAT‑AG‑AAAGACTTAA‑GCGGATGGTTCGTTGGTTTATATA‑CA‑ACCCCCCCTTCTG‑c                                                                                                                                                                                         >  2:74154/1‑228 (MQ=255)     gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAAGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGGGAT‑ATAG‑AAGACGAATGCGATTGAGATTAAAACACCGGCGGTAATCGCC‑AGAA‑AG‑AAAGTCTTTA‑GCGGGTGTTACGTTGGTTTATAAA‑CA‑CCCCCCCCCTCGG‑cccat                                                                                                                                                                                     >  2:483950/1‑225 (MQ=255)              aaGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATGGAGATGAAAACACCGGCGGTAATCGCC‑AGAT‑AG‑AAAGGATTAA‑ACGGAGGTTTTGTTGGTTTATAAA‑CA‑CCCCCCCCTTTGG‑CCACTTTGGCCCtt                                                                                                                                                                            >  2:207747/1‑228 (MQ=255)                      ggCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGCATTGA‑GCCTGTG‑CCAGGGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCC‑CGAT‑AG‑AAAGACTTAA‑AAGGCTGGTTCGTTGCTTTTTATA‑AA‑CACGCCCCTTCTG‑CCACTTTGTGCATTgtagactg                                                                                                                                                                    >  1:492620/1‑224 (MQ=255)                      gcctccgaaactaacaAAACGAATCAGCCCAAGAATGAAG‑AAGTGCCCGCCACCCATTTGCGACGGGCCG‑AGGGGAAGTGT‑GCC‑GTG‑GCCTTGGCTTGCTGTGAT‑AAAG‑AAGACAAATGAGATGGAAAGTAAACCACGGGGGTAAACCGCC‑AGAT‑AG‑AACCTGTTAT‑GGGGTTTTTCGGTTTTTTTATAGA‑CC‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAgg                                                                                                                                                                    <  1:126679/216‑1 (MQ=255)                        cTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGGTGAAAAAACCGGCGGTAAACGCC‑AGGT‑AG‑AAAGCATTAA‑GAGGAGGTTTCGTTGGTTTTTAGA‑AA‑ACCCCCCCTTCGGCCCCCTTTGGCCTTTgacacggaa                                                                                                                                                                   >  1:508555/1‑220 (MQ=255)                               aaaaacaCAAAGAATCAGCCCCAGAGGGAAGCAAATG‑CCGCAAACCGGTTTTGCCAGGACG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCC‑AGAT‑AG‑AAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CA‑CCCGCCTCCTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                            <  1:244140/224‑1 (MQ=255)                                  aacaCAAAGAATCACCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCAATGCCG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATACCC‑CGCT‑CG‑AAAAGCCTAC‑ACGGATGTTTCTCTTCGTCATAAA‑CA‑CACCCCCCCTCCG‑ACACTTTGTCCATTTGCGCAGGAATAaaaatgaa                                                                                                                                                        >  1:352527/1‑224 (MQ=255)                                     acaAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCC‑AGAT‑AG‑AAAGTCTTAA‑ACGGAGGTTTCGTGGCTTTATAAC‑CA‑CCCGCCCCTTCGG‑CCACTTTGGCCCtttctgaggggataaaaaaaacaa                                                                                                                                                      >  1:499310/1‑205 (MQ=255)                                             tcgaccacAGGGAGAAGCATATC‑CCGCCAACAAGTTGTGACTTGCC‑‑AAGGGGAATGG‑GGCTGGGTAGAGTGG‑TTGTTATGAA‑AGAG‑AAGAGTAATGCTTTTGGGTGGAAACCACCGTCGTTAACCGCC‑TTAT‑AG‑AAAGTCTTAA‑GCGGATGTTTGGTTCTTTTATAAA‑AA‑CTCGCCTCTTCGG‑CCCCTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGt                                                                                                                                             <  2:324801/221‑1 (MQ=255)                                                    gggagagaAAC‑AAGG‑CCCCCAACCTGTTTTTGCAACGCG‑AGGGGGTTTGTGGCCGGTG‑CCGGTGG‑TGGCTGTAAA‑AAAA‑AAAAGAAAGGGGTTTGAGGATAAAACCCGGGCGGTAATCCAC‑AAAA‑AA‑AAAGTTTTAG‑GGGGTTTTTTCGTTTCTTTAAAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCttt                                                                                                                                      <  1:406481/226‑1 (MQ=255)                                                     agagagAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGAATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATTGTC‑CGAA‑AGAAAAGACTTAG‑GGGGA‑GTGTCTTTGCTTTATAGA‑AA‑CCCCCCCCTTCGG‑ACACTTTGGCCAATGTAGCAGGAAATaaaaaaaccacaagtgtgtacggtg                                                                                                                                       >  1:324801/1‑205 (MQ=255)                                                               caatatgcccccctattttttTCCATGCCGAAAGTGTATTTG‑GACTGTG‑CGAGTGG‑TTG‑GGTGATAAAAG‑AACAGGATTCGTATGGGGCTGAACACCCAGACGGTAATCGCA‑AGAT‑AA‑AAGTTCTAAA‑GGGGGTTTTTCGTTGCTTTATAAA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTGCACACTAAc                                                                                                                             <  1:222402/211‑1 (MQ=255)                                                                   gcccgcaAACAAGTTTTGCCATGCGG‑AAGGGATTTGG‑CCCGGTG‑CCAGTGT‑TTGGTGTAAT‑ATAA‑AAGACGAATGCGATTGAGATAAAACCCCCGGGGGTAATCGCC‑AGAT‑AG‑AAAGTCTTAA‑GCGGATGTTTCTTTGCTTTATAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc                                                                                                                         <  2:508555/225‑1 (MQ=255)                                                                                          gccgggggcgaTTGT‑GCCTTGG‑CTAGTTG‑TTGTTG‑GAA‑AAAGAAAAACGAATGGTGTTGAAATGAAAACCCGGGGGGTACTCGCC‑AGAG‑AA‑AAAGTCTTAA‑GCGGTTGTTTCGTTGTTTTAAAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACgaga                                                                                                  <  1:501341/218‑1 (MQ=255)                                                                                            caaaagggatgtggccctggg‑ctAGTGG‑TTGCTGTGAT‑AAGG‑AAGGCGAAGGGGATTGAGATAAAAACCCCGGCGGTAAGCGCC‑AGAT‑AG‑AAATTCTTAA‑GGGGATGTTTCTTTGCTTTATAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  2:327978/213‑1 (MQ=255)                                                                                                     ttGT‑GCCTGTG‑CAAGTGG‑TGGCTTGGAT‑AGAG‑AAGCCGAAGGCGATTGAAAGGAAAACCCCGGCTGTAATCCCC‑ATAT‑AA‑AATGTCTAAA‑GGGGATGTTGCGTTGCTTGATAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                                                                       <  2:35485/228‑1 (MQ=255)                                                                                                                          gtgtgtggtaaaaaagaaagagAATGCGATAGAGAGGAACACCCGGGAGTTATTACAC‑AGAA‑AG‑AAAGTGTAAA‑GGGGATGTTTCTTTGCTTTAAAGC‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                       <  1:70687/207‑1 (MQ=255)                                                                                                                          cagccagtgaatgag‑aagaCGAATGCGATTGAGATGAAAACACCGGCGCTAAGCGCC‑AGAT‑AG‑AAAGTCTTAA‑GCGGATTTTTCGTGGTTTTATAAC‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                       <  2:617441/214‑1 (MQ=255)                                                                                                                               gtgAT‑AAAA‑AAGAAGAAGCGGTTTGAGTTAAAAACCCCGGCGTAAATC‑CA‑AGGTAAG‑AAAGTCATAA‑GGGTATGTTTCGTTGCTTTATAGA‑CC‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  2:399833/227‑1 (MQ=255)                                                                                                                                            gaCGAATGCGATTGAGATGAAAACACCGGCGGGAATCGCC‑AGAT‑TG‑AAAGTCTTAA‑GCGGATGGTTCGTTGCTTTATAGA‑CA‑CCCGCCTCTTCGG‑CCACCTTGGCCATTTCAGCAGGAAGTAAAAAATCAAAAGGGTTGTCCGCTTTCACCCTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCatattt                                                     >  2:570826/1‑226 (MQ=255)                                                                                                                                               gAATGCTATTGAAATAAAACCCCCGCGGGAAAACCCA‑AGAT‑AA‑AAACTCTAAA‑GGTTTTGTTCCGTTGCTTTATAGA‑CG‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGc                                                  <  2:499310/228‑1 (MQ=255)                                                                                                                                                 aaaaattttagagaaaaaCCCCGGGGGTAACCCCCAAAAT‑AG‑AAA‑TGTTAAGGAGGGTTTTGCGTTGTTTAACAGC‑CC‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  1:106802/214‑1 (MQ=255)                                                                                                                                                                 aaCACCGGCGGTAATCGCC‑AGAT‑AG‑AAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATCGGAAATTGGTATAGATACTATCTCCCCCTGCTTATAGGctgcc                                 >  2:106802/1‑223 (MQ=255)                                                                                                                                                                                         gaaagaattctttaacggttttttttGTTGTTTTAAAGA‑CC‑ACCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGc             <  2:492620/205‑1 (MQ=255)                                                                                                                                                                                            aaaatccaaaa‑GGGGATTTTTAGTTGATTTAAAAC‑CA‑CCCGTCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa        <  2:590052/219‑1 (MQ=255)                                                                                                                                                                                              gaagttctacgtcgtaggtctgTTGCTTTAAAAG‑CACCCCGCCTCTTCGG‑ACTCTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGAGGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa        <  2:588224/207‑1 (MQ=255)                                                                                                                                                                                                    ataaggggtttgtttgtttgtttATAAACCA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt  <  1:570826/208‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                            CACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCG‑AAGGGCATTGT‑GCCTGTG‑CCAGTGG‑TTGCTGTGAT‑ATAG‑AAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCC‑AGAT‑AG‑AAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CA‑CCCGCCTCTTCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGT  >  NC_000913/954185‑954595

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