Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,571,384	0	T	A	43.8%	37.4 / 12.7	33	Y64F (TAT→TTT)	yjiR	putative DNA‑binding transcriptional regulator/putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base A (14/0);  total (23/10)
Fisher's exact test for biased strand distribution p-value = 1.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTGCTCAATCC  >  NC_000913/4571239‑4571534                                                                                                                                                      |                                                                                                                                                       gCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAAcc                                                                                                                                                    <  1:25169/150‑1 (MQ=255)          tAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAAcgcg                                                                                                                                            <  1:333171/149‑1 (MQ=255)          tAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCgg                                                                                                                                           <  1:693297/150‑1 (MQ=255)          tAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCgg                                                                                                                                           <  2:70427/150‑1 (MQ=255)                  cAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGGTTTTCGTTGT‑TGGA‑CAAAAAAAACCCAACGCGGCTGCGGAg                                                                                                                                   >  2:549041/1‑150 (MQ=255)                      ggAACTGTCGCTACGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGGGTTTCGTTGT‑GGGG‑CAAAATAAACCCAACGCgggggggggggggtg                                                                                                                              >  1:296030/1‑138 (MQ=255)                        aaCTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTAGCGG‑CT‑GGGCTTTTCGGTGT‑TCGG‑CAAAAAAAACCGCAAGCGGCCGCGGAggg                                                                                                                                 >  1:466033/1‑144 (MQ=255)                              cGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GT‑GGGGTTTTCGTTGT‑GCGG‑CAAAAAAAACCGAAAGCGGCTGGGGGGGGGTGAGGtt                                                                                                                         >  1:140819/1‑148 (MQ=255)                              cGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCAt                                                                                                                      <  1:287429/151‑1 (MQ=255)                                cTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GC‑GGGCTTTTCGTTGT‑GCGG‑CAAAATAACCCGCACGCGGCTGCGGGGGGATGAGGtt                                                                                                                         >  1:301267/1‑146 (MQ=255)                                 tATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGGTTTTCGTTTG‑GGGACCAAAAAAACCCG‑ACGCCGCCGCGGGGGGATTAGCtt                                                                                                                         >  1:598666/1‑145 (MQ=255)                                 tATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCGTTTAGTTGTGGGGG‑CAAAAAAACCC‑AAAGCGGCCGCGGAGGGAAGAGCTTCAACGt                                                                                                                   >  1:150076/1‑151 (MQ=255)                                  aTGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCg                                                                                                                    <  1:531266/149‑1 (MQ=255)                                     cgcCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTTCGTTGT‑GCGG‑CAAAAAAAACCCAACGCGGCCGCGGGGTGGTGAGCTTCATCGTCTcc                                                                                                               >  1:430803/1‑151 (MQ=255)                                         tCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGGAATGGCGG‑GG‑GGGCCTTTCGGTAT‑GGGG‑CAAAAAAAAACGAAAGCGGCCGCGGGGGGATGAGCTTCATCggt                                                                                                                  >  2:179585/1‑142 (MQ=255)                                          ccAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGACCTGGCGGCCG‑GGGCCTTCCGTT‑T‑GGGA‑CCAAAAAACCCAGACGCGGCGGCGGGGGGGTGAGATTctttctcttccga                                                                                                            >  2:193387/1‑137 (MQ=255)                                                tatCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGG‑CTGGGGCTTTA‑GTTGT‑GCGA‑CAAAAAAAACCCAACGCGGCTGCGGGGGGATGAGGATCATCGTCTcccacgccttat                                                                                                     >  2:375874/1‑140 (MQ=255)                                                 atCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGAGATGAGCTTCATCGTCTCCAGCGTCTGATa                                                                                                    <  1:547083/150‑1 (MQ=255)                                                   ccAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CC‑GGGGCTTACCTTGT‑GCGG‑CAAAAAAAACCGCACGCGGCTGCGGGGGGGTGAGCTTCATCCTCTCCAGCGGCTGATa                                                                                                    >  1:30020/1‑148 (MQ=255)                                                      gCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctg                                                                                              <  2:401781/151‑1 (MQ=255)                                                       cACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GC‑GGGGTTTACGTTGT‑GCGG‑GAAAACAAACCCAACGCGGCCGCGGGGGGGTGAGCTTTATCGTTTCCACCGTCTGATACGcc                                                                                                >  2:578430/1‑148 (MQ=255)                                                                         aaTTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACAGGCGG‑CG‑AGGGTTTTTGTTGT‑GCGG‑CAAAAAAACACCGACGCGGGCGCGGGGTGGGGAGATTCATCTCCGCCAGCGTCTGATACCCCTTGTGCACgggggggg                                                                                >  2:386969/1‑140 (MQ=255)                                                                          aTTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CC‑GGGCTTTTCGTTGT‑GGGA‑CAAAAAAAACCCCACGCGGCGGCGGGGGGATGAGCTTCATCGTCTCCCGCGGCTGATACGCCTGCTGCACGGGGGTGATa                                                                              >  1:113603/1‑147 (MQ=255)                                                                            ttCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGAc                                                                         <  1:657483/150‑1 (MQ=255)                                                                                   gggCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGGTTTTCGGTGT‑GCGG‑CAAAAAAAACCCCACGCGGCTGCGGGGGGATGAGCTTTATCGTCTCCAGCGGCTGATACGCCTGCTGCCCGGGGGTGATACTTACGcc                                                                      >  1:457147/1‑146 (MQ=255)                                                                                         tGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GT‑GGGGTTTATGTTGG‑GGGC‑CAAAAAAAACCGGACACGGCGGCGGAGGGGTGGGGGTTTTTGTCCCCCGCGGCCGATACGGCCGCTGTCCGGGGCCGATACTGACGCCGGGCTCCCGAc                                                           >  2:172138/1‑151 (MQ=255)                                                                                               ggACGCGTCATCGGCGGTACTGGCGG‑CC‑GGGCTTTACGGTGT‑GCGG‑CAAAATAAACCCAACGCGGCCGCGGGGGGGTGGGCTTCATCGTCTCCCACGTCTGATACCCCCGCTTCACCGGGCTGGTAATGAAGCCGTGCTCCTGACTTaa                                                       >  1:642836/1‑149 (MQ=255)                                                                                                               gTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCa                                      >  1:304363/1‑150 (MQ=255)                                                                                                                 aCTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATAAGCCCGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAtt                                   >  2:387294/1‑151 (MQ=255)                                                                                                                   tGGCGG‑CT‑GGGCTTTACGCTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTtc                                 >  1:657129/1‑151 (MQ=255)                                                                                                                    ggcgg‑cT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAATtc                                 <  2:704867/150‑1 (MQ=255)                                                                                                                           t‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGATGATACTGACGCCGTGCTACTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGt                          >  2:473566/1‑151 (MQ=255)                                                                                                                                             cGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTg          >  1:60201/1‑151 (MQ=255)                                                                                                                                                      tAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTGCTCAATcc  >  2:140282/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       GCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CT‑GGGCTTTACGTTGT‑GCGA‑CAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTGCTCAATCC  >  NC_000913/4571239‑4571534

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