Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,689	0	G	T	48.9%	11.8 / 23.5	46	intergenic (+43/+22)	yegJ/yegK	uncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (14/9);  new base T (8/14);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGC  >  NC_000913/2151543‑2151836                                                                                                                                                      |                                                                                                                                                       ataCGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑tt                                                                                                                                                   <  2:517655/151‑1 (MQ=255)    cgcgcgGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑Caa                                                                                                                                               <  1:762276/151‑1 (MQ=255)      cgcgGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TA‑CAA‑cc                                                                                                                                            <  2:754550/151‑1 (MQ=255)                            gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGGATTGTTGTC‑TT‑TAA‑CCGCCATGAGAAACTCACTCTTcc                                                                                                                      >  1:24027/1‑151 (MQ=255)                             ttATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAAATCACTCTGccc                                                                                                                     >  1:805796/1‑151 (MQ=255)                               ataCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTTTTGTC‑TT‑TCA‑CCGCCATGAGAAAATCACTCCGCCca                                                                                                                    >  2:622248/1‑150 (MQ=255)                               ataCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAAATTTTGTT‑TT‑CAA‑CCGCCCTGAGAAACTCACTCTGCCcac                                                                                                                   >  1:1184288/1‑151 (MQ=255)                                  ccttttaTCTTTTATGGTTCCCCAATCAAAAAAAAAAAGCCGACGGCATAAGAGAAAAATTA‑TTTCGCTACTCACCCTTTATATTTATTTTTAAAC‑CG‑GGCGTTTAATTTTTTTTC‑TA‑AAG‑CCGCCATGAGAAACTCACTCTg                                                                                                                        <  2:861587/137‑1 (MQ=255)                                  catttgtgcttttATGGGCCCCCTATCAAAAAGAAGAAGCCGCCGATCAAAGATAAAAATTG‑GTTCGCTAATCAGCCTTTATATTTATTTTTAAAC‑CG‑GGCGGTTGAATTGTTTTT‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAg                                                                                                                 <  1:998013/140‑1 (MQ=255)                                   ctttGATCCTTTATGGGTCACCTACCAAAAAGAAGAGGCCGACGGCAAAAGAGCAAAATTA‑TTCCGTTAACCACCCTTTATATTTTTTGTTAAAG‑CG‑GGGGGTTAAATTGTTGTC‑TA‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAgc                                                                                                                <  1:104682/148‑1 (MQ=255)                                   cATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑CCGCCCTGAGAAACTCACTCTGCCcaca                                                                                                                  >  1:1005164/1‑148 (MQ=255)                                      tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑ACGCCATGAGAAACTCACTCTTCCCACAGCGcc                                                                                                             >  2:921939/1‑150 (MQ=255)                                                ttGGTCCCCTATCAAAAAAAAAAAGCCGAGGAGCAAAGAAAAAAATTA‑TTTCGCTAACCAGCCGTTATTTTTTTTGTTAAAC‑CG‑GGCGGTTAATTTGTTTTC‑TC‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACg                                                                                                   <  1:497112/150‑1 (MQ=255)                                                 tGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTTAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                 >  1:152180/1‑151 (MQ=255)                                                   gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTTTC‑TT‑CAA‑ACGCCCTGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGtttt                                                                                               >  1:299290/1‑151 (MQ=255)                                                      aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTTTTTTC‑TT‑TAA‑CCGCCCTGAGAAACTCACTCTGCCCCCAGCGCCCGAGCCAACGtttt                                                                                               >  2:1119493/1‑148 (MQ=255)                                                               aaaGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTTAATTGTTTTC‑TT‑TTA‑CCGCCTTGAGGAAAACACTCTGCCCCCAGCCCCCATGCCAACATTTTTTCTTCAAcc                                                                                     >  2:803345/1‑149 (MQ=255)                                                                aagaagaggCCGACGAGCAAAGAGAAAAATGA‑GTTCGCTAATCAGCCTTTATTTTTTTTTTTAAAC‑CG‑GGCGGTTGATTTGTTGTC‑TC‑AAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTa                                                                                  <  1:714065/151‑1 (MQ=255)                                                                              gaGCTAAAAGCAAAATTG‑TCGCGTTACTCAGCCGTTATATTTATTTTTTAAC‑CG‑GGCGGTTGAAGTGTTGTC‑TA‑CAT‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTa                                                                        <  2:406929/147‑1 (MQ=255)                                                                                aaaaaaaagaaaaTTATGTTCTTTAATCATCAGTTATTTTTATTTTTTAAC‑CG‑GGGGGTTTATTTTTGGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccg                                                                     <  1:299528/140‑1 (MQ=255)                                                                                  tAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTTTC‑TT‑TCA‑CCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTTTCATCATCCGTACGCTCGTTTAACGCCggg                                                                >  2:892666/1‑151 (MQ=255)                                                                                  tAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCggg                                                                >  2:335767/1‑151 (MQ=255)                                                                                   aagaaaaaaaTTT‑TTTCGCTAATCAGCTGTTATTTTTATTTTTAAAC‑CG‑GGGGGTTGATTTGTTTTC‑TA‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCgg                                                                 <  2:996433/143‑1 (MQ=255)                                                                                   aaaagaaaaaTTG‑GTTCATTAAGCAGCTTTTATTTTTATTGTTTAAA‑CG‑GGCGGTAAATTTGTTTTC‑TA‑AAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccgccg                                                                  <  1:304487/142‑1 (MQ=255)                                                                                   aaGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTTAATTGTTTTT‑TT‑TAA‑CCGCCATTAGAAAATCACTCTGCCCCCAGCGCCAGTGGCAAAGTTTTTTCATCATCCCTACCCCCGTTTACCGCCcgg                                                                >  2:399005/1‑150 (MQ=255)                                                                                     gagCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTTAATTTTTTTC‑TT‑TAA‑ACGCCATGAGAAAAACACTCTTCCCACAGCGCCAGTGCCCACGTTTTATTATCATCCCTACGCTCGTTTACCGCCggggcg                                                             >  2:161046/1‑148 (MQ=255)                                                                                      aaaaaaaTTG‑CTTCGATAACCAGCGTTTATTTTTTTTTTTAAAA‑CG‑GGTGGTTAATTTTTTTTC‑TT‑CAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTg                                                             <  2:1090387/147‑1 (MQ=255)                                                                                         aaaaTTA‑GTCCCCTAACCGCCCTTTTTTTTTTTTGTTAAAG‑CG‑GGGGTTATAATTTTTCTC‑TC‑CAA‑CCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGAc                                                          <  2:444470/150‑1 (MQ=255)                                                                                          aaatttttctCGCTAGTCAGCCTTTATATTTATTGTTTAAC‑CG‑GGCGGTTGATTTGTTGTT‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGAc                                                          <  2:179236/141‑1 (MQ=255)                                                                                          aaaTTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGATTAATTTTTGTC‑TT‑TAA‑CCGACATTAGAAAATCACTCTGCCCACCCCGCCCTTGCCAACCGTTTATAATCATCCGTACGCCCGTTTTTCGCCGGGGTgg                                                            >  2:1095914/1‑147 (MQ=255)                                                                                                    cGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTTTTTTC‑TT‑CAA‑CCGCCTTGAGAAACTCCCTCTGCCCGCAGCGCCAGTGCCAACGTTTTCTTATCATCCGTACGCTCGTGTAACGCCGGGGGGGACAAAAACTGCtt                                               >  1:627274/1‑151 (MQ=255)                                                                                                      cTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGTATTGGTGTT‑TT‑CAA‑CCCCCATGAGAAACTCACTCTGCCCACCGCGCCACAGCCAAAGATTTATCATCATACGGACGCGCGTTTTCCGCCGGGCTGGACAAAAACCGCtttt                                             >  1:13833/1‑150 (MQ=255)                                                                                                         aTCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑ACGCCATTAGAAAATCACTCTGCCCACAGCGCCAGTGCCAAAGGTTTTTTATCATCCCTACGCTCGGTTACCGCCCGGCTGGACAAAAACTGCtttta                                            >  2:213706/1‑147 (MQ=255)                                                                                                           tcagctgttttttttttttttAAA‑AG‑GGGGGTTTAATTTTTTTT‑CT‑AAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACa                                          <  1:68988/139‑1 (MQ=255)                                                                                                              cgtttATTTTTTTTTTTTAAA‑CT‑GGGGTTTAATTTTTTTTC‑TA‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTa                                             <  1:1041513/140‑1 (MQ=255)                                                                                                                ttttttttttttttttAAACCG‑GGCGGTTAATTTTTTTTT‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAAtt                                       <  2:804405/142‑1 (MQ=255)                                                                                                                       tttttttttAAA‑CCGGGCGTTTAATTTTTTTTC‑TC‑CAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                                 <  2:979282/141‑1 (MQ=255)                                                                                                                       tttttttttAAA‑CC‑GGGCGTTTAATTTTTTTTTTT‑CAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGt                                  <  1:1197100/140‑1 (MQ=255)                                                                                                                       tttttGTTAAAC‑GT‑GGGGGTTAATTTTTTTTC‑TC‑CAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaataa                             <  2:751722/147‑1 (MQ=255)                                                                                                                             ttttaaaacggtggtttGATTTTTTGTT‑TA‑AAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCg                         <  1:985114/135‑1 (MQ=255)                                                                                                                              ttAAA‑CT‑GGCGGTTGAATTGTTTTC‑TT‑TAA‑ACGCCATGAGAAAATCCCTCTGCCCACAGCGCCCGTGCCAACGTTTTATCATTATTCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTAAAAATTCAGGGAATTAATCGAgt                      >  2:864909/1‑150 (MQ=255)                                                                                                                                   accgggCGGTTGATTTTTTTTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgt                      <  2:542205/143‑1 (MQ=255)                                                                                                                                            ttGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtg         >  1:500472/1‑151 (MQ=255)                                                                                                                                             ttgaTTTTTTTC‑TTCCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGt           <  2:722932/146‑1 (MQ=255)                                                                                                                                             tGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgc        >  1:1031826/1‑151 (MQ=255)                                                                                                                                               atttgtttTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtg         <  1:600484/146‑1 (MQ=255)                                                                                                                                                aTTGTTGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctgc     >  2:1033109/1‑151 (MQ=255)                                                                                                                                                   gttgtC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccgc  >  1:535145/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       ATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGC  >  NC_000913/2151543‑2151836

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