Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,151,689	G→T	31.7%	intergenic (+43/+22)	yegJ → / ← yegK	uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,689	0	G	T	31.7%	16.4 / 4.2	31	intergenic (+43/+22)	yegJ/yegK	uncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (10/10);  new base T (6/4);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

AATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTT‑ACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCT  >  NC_000913/2151492‑2151927                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  aaaggggaagggAGGGAAAAAACAAAAAAAAAGTAATTGAACACGGGTTTTATACGCG‑GGGGAAAAAATTTAACGCCGGTTTTCCATTTGACCTTTTATTGTCCCCCAATCAAAAAAAAAAAGCGGACGAGTAAAAAGAAAATTTAGTTCGCTAA‑TCAGCCGTTATATTTTTTTTTAAAA‑CTGGCGG‑TTAAA‑TTGTTGTT‑TC‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCa                                                                                                                                                                                                       <  1:1037307/229‑1 (MQ=255)             agtaaaaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATCCGCGGGGGAAAAAAATTT‑GCGGCGTTTATACCATTGATCCTTTATTGGTCCCCTACCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTATTTCGTTAA‑TCCGCCGTTATTTTTTTTGTTTAAA‑CGGGGGG‑TTAAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                            <  1:1060838/233‑1 (MQ=255)             aGTGAAAAACAAAGAAAGAAGATATTTCAGACTGGATGTATCCGCGGGGGGATAAAATTT‑ACGGCGTTTATCCCTTTGATCTTTTATTGGTCCCCTACCCAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTTTTGTTTAAA‑CTGGGGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                            <  1:1001456/238‑1 (MQ=255)                   aaaaaaaaaaaaaaGATATTGCAGACTGGTTTTTATCGCGGGGTGAAAAAATTG‑ACGCCGTTTATCCCATTTCTCTTTTATTGGTCCCCTACCAAAAAGAAAAAGCCGACGAGCTAAGAGCAAATTTATTTCTCTAA‑CCAGCCTTTATATTTATTTTTTAAC‑CGGGCGT‑TGGAA‑TTTTTTTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                                                                                                                     <  2:446524/231‑1 (MQ=255)                                                                       tt‑ACGGCGTTTATACCATTGATCTTTTATGGGTCCCCTACCAAAAAAAAGAGCCCAACGAGTAAAAAGAAAAATTAGTCCGTTAG‑TCACCAGTTATATTTATTTTTTAAC‑CGGCGGG‑TTAAT‑TTGTTTTC‑TT‑A‑AACCGCCATGAGAAACTCACTCTGCCCAAAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                                 <  1:260894/239‑1 (MQ=255)                                                                                   acaCCTTTGATCTTTTATTGGTCACCTACACAAAAAAAGAAGCCGAGGAGAAAAGAGAAAATTTTTTTCGATAA‑TCAGCTGTTATTTTTATTGTTTAAA‑GGGGCGT‑TTGAA‑TTGTTGTC‑TA‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                                                                                                                      <  1:1032405/237‑1 (MQ=255)                                                                                          tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttg                                                                                                               >  1:233142/1‑239 (MQ=255)                                                                                          tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑ATGTTTTG‑TT‑C‑GAACCCCACGAGAAACTCACCATTCCCACCACGCCAGCGCCCACCTTTTATCATCATCCGTACACTCCTTTCCCGCCGGGCGGGACAAAAACAGCTTCAACAAATAAGGAGATAAATCGatgtg                                                                                                               >  2:805371/1‑235 (MQ=255)                                                                                              ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgtt                                                                                                             >  1:116902/1‑237 (MQ=255)                                                                                                  tATGGGCCACCTACCAAAAAGAAGAACCCGACGAGCTAAGAGAAAAATTAGTTCGCTAA‑CCAGCCTTTATATTTATTGTTTAAA‑CGGGCGT‑TTGAT‑TTGTTTTC‑TC‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCg                                                                                                       <  1:559861/239‑1 (MQ=255)                                                                                                      ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTTTTGGC‑TT‑T‑CACCGCCCTGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGCTTTTTCATTCTACGTACGCTCGCTTTCCGCCCGGCTTGACAAAAACTGCTTTAACAATTCAGGGAATAAATCCAGGTGTTTCTTCggttc                                                                                                   >  2:1020073/1‑235 (MQ=255)                                                                                                       gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑TTGTTGTC‑TT‑C‑TACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAAAGTTTTATCATTATCCGTACGCTCGTTTACCGACGGGCTGGGCACAAAACGCTTTCACAATTCAGGTAATAAATCGAGTTGTTCCTTCGtt                                                                                                     >  2:765178/1‑236 (MQ=255)                                                                                                       gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGGA‑TTGTTTTT‑TT‑C‑AACCGCCATGAGGAACTCACTCTGCCCACAGCGACAGCGCCAACGTTTTTTCATCATCCGTACGCCCGTTTACCGCCGGGCTGGGCCACAACTGCTTTTACAATTTCGGGAATAAAACGAGTTGGTCCttgggtgct                                                                                                  >  2:446165/1‑231 (MQ=255)                                                                                                       gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGGA‑ATTTTTTT‑TT‑T‑AACCGCCATGAGAAACTCACTCTTCCCCCAGCGCCCGTGCCCACCTTTTATCATTATCCGGAAGCTCCTTTACCACCCGGCCGGGCAAAAACTTCTTTTACAATTCAGGGAATAAATCGCGGTTGTCCTGCGGtgct                                                                                                  >  2:21152/1‑239 (MQ=255)                                                                                                        tCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTTTTGTT‑TT‑C‑AACCGCCCTGAGAAACTTACTCTGCCCACCACGCCAGTGCCAACCTTTTTTCATCATCCGTACGCCCGTTTTCCGCCCGGCCGGACAAAAACTGCTTTTACACATTAGGTAATAAATCCAGTTGTTTCtttgtttct                                                                                                  >  2:866416/1‑230 (MQ=255)                                                                                                            ttaaaaaaaaGAAAAAGCCGAGAAGACAAAAGAAAAATT‑GTTCCGTAACTCGGCCTTTAT‑TTTTTTTTTTAAA‑CCGGGGT‑GTGTA‑TTTTTTTC‑TTCC‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccgc                                                                                             <  2:116902/233‑1 (MQ=255)                                                                                                                 aaaaaaaAGAAGCCGCCGGGCAAAGAAAAAAATTAGTTCGCAAA‑ACACCCTTTAT‑TTTTTTTTTTAAA‑CCGGGGGTTTGAA‑TTTTTGTC‑TC‑A‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTc                                                                                        <  1:805371/238‑1 (MQ=255)                                                                                                                         gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTGA‑TTTTTGTC‑TT‑T‑AACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACCCCCCTTTACCCCCGGGCCGGGCAAAAACTCGTTTAACCATTCCGGGAATTAATCGAGTTGGTCCTTCGTTTCTGCCCCCCGTCCATTGgaa                                                                                >  2:860063/1‑239 (MQ=255)                                                                                                                         gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTTTC‑TT‑T‑AAACGCCATTAGAAACTCACTCTGCCCACAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCAGGCTGGGCAAAAACTGCTTTAACAATTCAGGTAATAAATAGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaa                                                                                >  1:746787/1‑239 (MQ=255)                                                                                                                          aagCGGACGAGTAAAGAGAAAATTTATTTCGTTAA‑CCAGCTGTTAT‑TTTTTTTTTTAAA‑AGGGGGG‑TTTAATTTGTTGTC‑TC‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaa                                                                               <  1:217792/239‑1 (MQ=255)                                                                                                                             cgaaggaGCAAAGAGAAAATTTGGTCCCCTAA‑TCAGCGGTTATATTTTTTGTTAAAA‑C‑GGGGG‑GTGAA‑TTTTTGTCTTT‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                                              <  1:765178/232‑1 (MQ=255)                                                                                                                                             aaaaTAAGCTCGCAAA‑CCAGCGTTTATATGTATTTTTTAAA‑CGGGGGG‑TTGAT‑TTGTGGTC‑TT‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGcac                                                            <  1:1126355/239‑1 (MQ=255)                                                                                                                                                  tcgtccgttAA‑TCATCTGTTATATTTCTTTTTAAAA‑CGGGGGG‑TTAAT‑TTGTT‑TC‑TT‑CTAACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                       <  1:426931/231‑1 (MQ=255)                                                                                                                                                    gTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑GTTTTTTC‑TT‑C‑AAACGCCCTTAGACAGTCAACCTGCCCACACCGGCCGTGGCCACGGTTTATCATCATCCCTAAGCCGCTGTACCCCCGGGCTGGACCAACACATCTTTTAAAATTCCAGTAATAAATCCAGTTTGTCCTATGTGTCTGGCGCCCGGCGCTTTAAAAAAAGGGTAAAGAAAGGCACAttcgga                                                     >  2:1142029/1‑234 (MQ=255)                                                                                                                                                    gTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTAACGCCGGGGCGGACAAAAACTGCTTTAACAATTTCGGTAATAAATCGCGTTGTTTCTGCGTTGCTTCCGCCAGTCCCTTGAAAAACGGGGGAAAAAAAAGCACATgagga                                                     >  1:218977/1‑239 (MQ=255)                                                                                                                                                                   cGTTATA‑TTATTTTTTAAACCGGGCGG‑TTGAT‑TTGTTGTC‑TA‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                         <  1:453539/237‑1 (MQ=255)                                                                                                                                                                     ttATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑TTTTTTTC‑TT‑C‑AAACGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCCACGTTTTTTTCTTATCCGTACACTCGTTTACCGCCGGGCCGGGCAAAAACTGCTTTAACAATTCAGGGAATAAATAGCGTTGTTCCTGCGTTGCTGCCGCCAGACCATTGGAAAACGGGGTAAAGAAAGGCACATGAGGGGGATTATCCAGCctat                                     >  2:191291/1‑239 (MQ=255)                                                                                                                                                                            tattGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTTCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAAGCACATGAGGAGAATTATCCAGCATATTTAACg                               >  1:1051334/1‑238 (MQ=255)                                                                                                                                                                               tGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCCCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTTCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGGGGAGAATTATCCAGCATATTTAACGCCAg                           >  1:152446/1‑239 (MQ=255)                                                                                                                                                                                         tGGCGG‑TTAAT‑TTTTTGTC‑TT‑C‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGa                   <  1:583138/238‑1 (MQ=255)                                                                                                                                                                                            cGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTaaag                                                                    >  1:120927/1‑186 (MQ=255)                                                                                                                                                                                                           gtC‑TT‑C‑ATCCGCCATGAGAAACTCCCTCTGCCCACAGCTCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAtgct  >  2:460266/1‑239 (MQ=255)                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  AATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTT‑ACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCT  >  NC_000913/2151492‑2151927

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