Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6890GT64.7% 12.8 / 6.5 35intergenic (+43/+22)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (9/3);  new base T (14/8);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAG  >  NC_000913/2151567‑2151839
                                                                                                                             |                                                                                                                                                       
gggtttatcctgttGATCTTTTATTGGTCCCCTCTAAAAAAGAAAAAGCCGACGAGATAAGAGCAAAATTAGTTCGCTAATCAGCCTTTATTTTTTTTTTTAAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑CAACCGCCATGAGAAactcactc                                                                                                                             <  1:131338/138‑1 (MQ=255)
  ggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTCTT‑CCAACGCCATGAGAAACTCAATCTg                                                                                                                           >  1:286579/1‑150 (MQ=255)
  ggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑CAACCGCCATGAGAAACTCACTCTg                                                                                                                           >  1:93806/1‑150 (MQ=255)
   gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTCTT‑TAACCGCCATGGGAAAATCACTCTGc                                                                                                                          >  2:102193/1‑150 (MQ=255)
   gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTCTT‑CAACCCCCCTGAGAAAATCACTCTGc                                                                                                                          >  2:253057/1‑150 (MQ=255)
         ccATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCCT‑TAAACCCCCTGAGAAACTCACCCTGccccccg                                                                                                                    >  2:100006/1‑146 (MQ=255)
               aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTT‑TTTTCTT‑CAACCGCCCTGAGAAAATACCTCTTCCCCCCGCGCCAGt                                                                                                             >  1:241593/1‑151 (MQ=255)
               aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTT‑TTTTCTT‑CAACCGCCATTAGGAACTCACTCTGCCCACAGCGCCCGt                                                                                                             >  1:401672/1‑151 (MQ=255)
                tcttttATTGCTCCCCTTACAAAAAGAAGAGGCCGACGAGCTAAGAGAAAATTTGTTTCGTTAATCAGCCTTTATTTTTTTTGTTTAAC‑CTGGCGTTTAAT‑TTT‑TTGTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTg                                                                                                            <  1:263105/148‑1 (MQ=255)
                tCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTT‑TTTTTTT‑CAACCCCCCTGAGAAACTCACCCTGGCCCCAGCGCCaagg                                                                                                            >  2:269141/1‑148 (MQ=255)
                      aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑CAACCGCCATGGGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACg                                                                                                      >  2:278486/1‑151 (MQ=255)
                       ttGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTTTA‑TTT‑TTGTCTT‑CAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGt                                                                                                     >  1:405956/1‑151 (MQ=255)
                       ttGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTTAA‑TTG‑TTGTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCCGTGCCAACGt                                                                                                     >  2:392900/1‑151 (MQ=255)
                         ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑CAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGttt                                                                                                   >  1:97001/1‑151 (MQ=255)
                              ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTTTT‑CAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCCGTGCCCACGTTTTttca                                                                                              >  2:376484/1‑151 (MQ=255)
                              ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTatca                                                                                              >  2:69222/1‑151 (MQ=255)
                                 aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑TAACCCCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAAAGTTTTttcatca                                                                                           >  2:141541/1‑151 (MQ=255)
                                 aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGATGAA‑TTG‑TTGTCTT‑CAACCCCCATGAGAAAATCACCCTGCCCACAGCCCCAGTGGCCAAGGTTTatcaata                                                                                           >  1:144888/1‑148 (MQ=255)
                                   caaaaaGAAAAAGCGTACGAGCAAAGAGAAAAATTAGTCCGCAATACCGCCTTTATATTTATGGTTAAAA‑CGGGGGTTAAAA‑TTG‑TTTTCTC‑AACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatcatc                                                                                         <  1:142958/149‑1 (MQ=255)
                                       aaaaagaagCCGAAGAGTTAAGAGCAAATTTGTTTCTCTATCCCGCCTTTATATTTTTTGTTTAAC‑CGGGCGTTGAAA‑TTG‑TTTTTTC‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTa                                                                                     <  2:206525/148‑1 (MQ=255)
                                          aagaacCCGACGAGCAAAGAAAAAAATTTTTTCTTAAATCAGCCTTTATTTTTTTTGTTTAAA‑CGGGCGGTTAAA‑TTT‑TTTTCTT‑AAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc                                                                                    <  2:276003/149‑1 (MQ=255)
                                                      aGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTTAA‑ATG‑TTTTCTT‑TCACCGCCCTGAGAAACCCCCTCTGCCCCCAGCGCCCGTGCCCACGTTTTTTCCTCATCCCTACGCTCGTTTAcccc                                                                       >  2:376289/1‑148 (MQ=255)
                                                            gagCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTCTT‑CAACCCCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAAAGTTTTTTCATCCTCCGTACCCTCGTTTTCCCCCggg                                                                   >  2:381229/1‑148 (MQ=255)
                                                              gCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTTTT‑CAACCGCCATGAGAAACTTACTCTGCCCACAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGa                                                              >  1:192584/1‑151 (MQ=255)
                                                               cAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTTAA‑TTT‑TTTTCTT‑CAACCCCCATGAGAAAATCACTCCGCCCACAACGCCAGTGCCAACGTTTTATCATCATCCCTACCCTCGTTTTCCCCCGGGCTGGa                                                              >  1:181065/1‑150 (MQ=255)
                                                                 aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTCTT‑CAACCGCCATGAGAAAATTACTCTGCCCACAGGGCCAGTGCCAAAGTTTTATCTTCATACGTAAGCTCCCTTACCGCCCGGCCGGACaa                                                           >  1:129556/1‑151 (MQ=255)
                                                                            cTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAAATGCTTTa                                                >  1:323399/1‑151 (MQ=255)
                                                                              aacccccctttATTTTTATTG‑TTAAA‑CCGGGGGTTGAATTTT‑TTTTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAc                                              <  2:397831/142‑1 (MQ=255)
                                                                              aaacacccGTTATATTTTTTG‑TTAAA‑CGGGGCGGTGAA‑TTGTTTAGTTT‑AAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAc                                              <  2:301361/145‑1 (MQ=255)
                                                                               aacAGCCTTTATTTTTTTTTTTAAAA‑CGGGCGTTTAAA‑TTT‑TTTTTTC‑CACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACa                                             <  2:201233/149‑1 (MQ=255)
                                                                                ccAGCCGTTATTTTTATTTTTAAAC‑CTGGCGTTTAAT‑TTT‑TTTTCTC‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACaa                                            <  1:89018/150‑1 (MQ=255)
                                                                                              ttttGTTAAAC‑CGGGGGTTTAAA‑TTG‑TTTTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTGCCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                              <  2:219330/149‑1 (MQ=255)
                                                                                               tttttttAAACCTGGCGGTAAAG‑TTT‑TTGTCTC‑AAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                              <  1:187640/147‑1 (MQ=255)
                                                                                                  gTTTAAA‑CTGGCGGTTGAA‑TTG‑TTTTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACCGCGCCAGTGCCAACGTTTTATTATCATCCGTACGCTCGTTTACCGCCGGGCCGGAGAAAAACTGCTTTAACAATTCTGGTAATAAATCGAg                          >  1:381247/1‑151 (MQ=255)
                                                                                                    ttAAA‑CCTGGGGTTTAA‑TTT‑TTTTTTTCCAACCGCCATGAGAAACTCGCTCTGCCCTCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgt                         <  2:129584/151‑1 (MQ=255)
                                                                                                                             gtCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAg  >  1:288675/1‑151 (MQ=255)
                                                                                                                             |                                                                                                                                                       
GCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGGTTGAA‑TTG‑TTGTCTT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAG  >  NC_000913/2151567‑2151839

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