Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,3840TA40.0% 43.3 / 12.6 30Y64F (TAT→TTT) yjiRputative DNA‑binding transcriptional regulator/putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (6/12);  new base A (12/0);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 3.34e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTGC  >  NC_000913/4571240‑4571527
                                                                                                                                                  |                                                                                                                                                
cTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GTGGGGTTTTCGTTGTGC‑GGCAAAAAAAacc                                                                                                                                             >  1:127793/1‑150 (MQ=255)
   ttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGG‑GACCAAAAAAccc                                                                                                                                             >  1:291466/1‑147 (MQ=255)
      cTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAACTGCGG‑CCCGGGGTTTCGTGGTGC‑GGCAaaaaaaaaccc                                                                                                                                           >  2:696555/1‑141 (MQ=255)
       ttAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCCCGGGCTTT‑CGGTGTGC‑GGCAaaaaaaaac                                                                                                                                             >  2:562536/1‑140 (MQ=255)
              aaCAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGGGG‑GTGGGGTTTTCGTTGTGG‑GACAAAAAAAACCGAAAGCGGCTGC‑gg                                                                                                                              >  1:366264/1‑150 (MQ=255)
                  aTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑TTGGGCTTTACGTTGTGC‑GACAAAATATCCCGAACGCGGCTGC‑GGAGtga                                                                                                                         <  2:261773/151‑1 (MQ=255)
                     gAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGtgatga                                                                                                                      <  1:906928/151‑1 (MQ=255)
                       aCTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GCGGGCGTTTCGTTTTGC‑GGCAAAAAAACCCGAACGCGGCGGC‑GGGGtggt                                                                                                                        >  2:259025/1‑145 (MQ=255)
                              cTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGGGG‑GTGGGGGTTTCGTTTTGG‑GGCAAAAAAAACCCCACACGGCTGC‑gggggggttgggctc                                                                                                                 >  1:576023/1‑134 (MQ=255)
                               tATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CGGGGCCTTTCGTTGTGG‑GACAAAAAAAACCCCACGCGGCCGC‑GGAGGGGTGAGGtt                                                                                                                  >  2:213792/1‑145 (MQ=255)
                                     ccTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GTGGGGGTTTCGGTGTGCGGACAAAAAAACCC‑AAAGCGGCCGC‑GGGGGGGTTGGCTTCATCCTcccc                                                                                                        >  1:753875/1‑149 (MQ=255)
                                         cAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGTCGCTGCCCCGGACGCGTCATCGGCGGTCCTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAg                                                                                                      <  1:910575/147‑1 (MQ=255)
                                          aGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GCGGGCCTTTCGTTTTGG‑GGCAAAAAAAACCCGACGCGGCCGC‑GGGGGGGTGAGGTTCCTCGTCTCCcgcg                                                                                                    >  1:646369/1‑148 (MQ=255)
                                                      aCCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGGACCGGCCG‑CCGGGCCTTACGCTTTGC‑GGCAAAAAAAACCCAACGCCGCCGG‑CGGGGGGGGAGCCTCATCGTcccccccctc                                                                                                  >  2:688181/1‑132 (MQ=255)
                                                          gATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGGACTGGCCG‑CTGGGCGTTACGTTCTGC‑GCCAAAAAAACCCGAAACCGGCTGCGGGAGT‑ATGGGCTTCTTCCTCTCCAGCCGCTGGTACCGCCGCTGCAc                                                                                 >  2:683509/1‑151 (MQ=255)
                                                                aCTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑GCGGGGTTTTCGTTGTGG‑GGCCAAAAAAACCCAACGCGGCCGG‑GGGGGGGTGAGCTTCATTCTCTCCAGCGGCTGATAAGCCCGCTGCACgggg                                                                             >  1:549264/1‑147 (MQ=255)
                                                                   gggTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGAt                                                                        <  1:596348/151‑1 (MQ=255)
                                                                          ttCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGAc                                                                  <  2:694471/150‑1 (MQ=255)
                                                                          ttCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGAc                                                                  <  2:616783/150‑1 (MQ=255)
                                                                                             ggACGCGTCATCGGCGGTACTGGCGG‑CTGGGGTTTACGCTGTGC‑GCCAAAAAAAACCCCACGCGGGCGC‑GGGGGGATGAGGTTCATCCTCTCCCGCGCCTGATACGCCTGCCCCCCGGGGGTGGTACTTACCCCCCCCTCCCGACttta                                                >  2:508061/1‑147 (MQ=255)
                                                                                                 gcgTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTgcg                                          <  1:276047/151‑1 (MQ=255)
                                                                                                            ggTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGc                                >  1:218094/1‑150 (MQ=255)
                                                                                                               aCTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAtt                            <  1:905988/151‑1 (MQ=255)
                                                                                                               aCTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAtt                            <  1:121711/151‑1 (MQ=255)
                                                                                                                         tGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGt                   >  1:66819/1‑151 (MQ=255)
                                                                                                                          gggCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCg                    <  2:536735/149‑1 (MQ=255)
                                                                                                                                tACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGAt             >  1:345314/1‑150 (MQ=255)
                                                                                                                                  cGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATAc           >  2:653318/1‑150 (MQ=255)
                                                                                                                                   gatgtgC‑GCCAAAATAACCCGAACGCGGCGGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATAc           >  2:911374/3‑149 (MQ=255)
                                                                                                                                      gtgC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGcc       <  1:619273/150‑1 (MQ=255)
                                                                                                                                       tgC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGccc      >  2:600313/1‑150 (MQ=255)
                                                                                                                                           gACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTGc  <  2:526834/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGG‑CTGGGCTTTACGTTGTGC‑GACAAAATAACCCGAACGCGGCTGC‑GGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTGC  >  NC_000913/4571240‑4571527

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