Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,686	0	G	T	41.3%	35.6 / 15.1	46	intergenic (+40/+25)	yegJ/yegK	uncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base T (6/13);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 7.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGC  >  NC_000913/2151538‑2151836                                                                                                                                                        |                                                                                                                                                          gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACTAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAAttgtt                                                                                                                                                         <  2:296292/151‑1 (MQ=255) gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAAttgtt                                                                                                                                                         <  1:967868/151‑1 (MQ=255)      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)                                                                                                                                                        |                                                                                                                                                          GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTA‑GTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CT‑GGCGGTTGAATTGTTGTC‑TT‑CAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGC  >  NC_000913/2151538‑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: ‑