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

GCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAG  >  NC_000913/2151567‑2151839
                                                                                                                             |                                                                                                                                                       
gggtttatcctgttGATCTTTTATTGGTCCCCTCTAAAAAAGAAAAAGCCGACGAGATAAGAGCAAAATTAGTTCGCTAATCAGCCTTTATTTTTTTTTTTAAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAactcactc                                                                                                                             <  1:131338/138‑1 (MQ=255)
  ggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCCAACGCCATGAGAAACTCAATCTg                                                                                                                           >  1:286579/1‑150 (MQ=255)
  ggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTg                                                                                                                           >  1:93806/1‑150 (MQ=255)
   gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTTAACCGCCATGGGAAAATCACTCTGc                                                                                                                          >  2:102193/1‑150 (MQ=255)
   gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCCCCCTGAGAAAATCACTCTGc                                                                                                                          >  2:253057/1‑150 (MQ=255)
         ccATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCCTTAAACCCCCTGAGAAACTCACCCTGccccccg                                                                                                                    >  2:100006/1‑146 (MQ=255)
               aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCGCCCTGAGAAAATACCTCTTCCCCCCGCGCCAGt                                                                                                             >  1:241593/1‑151 (MQ=255)
               aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCGCCATTAGGAACTCACTCTGCCCACAGCGCCCGt                                                                                                             >  1:401672/1‑151 (MQ=255)
                tcttttATTGCTCCCCTTACAAAAAGAAGAGGCCGACGAGCTAAGAGAAAATTTGTTTCGTTAATCAGCCTTTATTTTTTTTGTTTAACCTGGCGTTTAATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTg                                                                                                            <  1:263105/148‑1 (MQ=255)
                tCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTTTTCAACCCCCCTGAGAAACTCACCCTGGCCCCAGCGCCaagg                                                                                                            >  2:269141/1‑148 (MQ=255)
                      aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGGGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACg                                                                                                      >  2:278486/1‑151 (MQ=255)
                       ttGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTTATTTTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGt                                                                                                     >  1:405956/1‑151 (MQ=255)
                       ttGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCCGTGCCAACGt                                                                                                     >  2:392900/1‑151 (MQ=255)
                         ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGttt                                                                                                   >  1:97001/1‑151 (MQ=255)
                              ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCCGTGCCCACGTTTTttca                                                                                              >  2:376484/1‑151 (MQ=255)
                              ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTatca                                                                                              >  2:69222/1‑151 (MQ=255)
                                 aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCCCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAAAGTTTTttcatca                                                                                           >  2:141541/1‑151 (MQ=255)
                                 aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGATGAATTGTTGTCTTCAACCCCCATGAGAAAATCACCCTGCCCACAGCCCCAGTGGCCAAGGTTTatcaata                                                                                           >  1:144888/1‑148 (MQ=255)
                                   caaaaaGAAAAAGCGTACGAGCAAAGAGAAAAATTAGTCCGCAATACCGCCTTTATATTTATGGTTAAAACGGGGGTTAAAATTGTTTTCTCAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatcatc                                                                                         <  1:142958/149‑1 (MQ=255)
                                       aaaaagaagCCGAAGAGTTAAGAGCAAATTTGTTTCTCTATCCCGCCTTTATATTTTTTGTTTAACCGGGCGTTGAAATTGTTTTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTa                                                                                     <  2:206525/148‑1 (MQ=255)
                                          aagaacCCGACGAGCAAAGAAAAAAATTTTTTCTTAAATCAGCCTTTATTTTTTTTGTTTAAACGGGCGGTTAAATTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc                                                                                    <  2:276003/149‑1 (MQ=255)
                                                      aGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAAATGTTTTCTTTCACCGCCCTGAGAAACCCCCTCTGCCCCCAGCGCCCGTGCCCACGTTTTTTCCTCATCCCTACGCTCGTTTAcccc                                                                       >  2:376289/1‑148 (MQ=255)
                                                            gagCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCCCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAAAGTTTTTTCATCCTCCGTACCCTCGTTTTCCCCCggg                                                                   >  2:381229/1‑148 (MQ=255)
                                                              gCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCATGAGAAACTTACTCTGCCCACAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGa                                                              >  1:192584/1‑151 (MQ=255)
                                                               cAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTCTTCAACCCCCATGAGAAAATCACTCCGCCCACAACGCCAGTGCCAACGTTTTATCATCATCCCTACCCTCGTTTTCCCCCGGGCTGGa                                                              >  1:181065/1‑150 (MQ=255)
                                                                 aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAAATTACTCTGCCCACAGGGCCAGTGCCAAAGTTTTATCTTCATACGTAAGCTCCCTTACCGCCCGGCCGGACaa                                                           >  1:129556/1‑151 (MQ=255)
                                                                            cTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAAATGCTTTa                                                >  1:323399/1‑151 (MQ=255)
                                                                              aacccccctttATTTTTATTGTTAAACCGGGGGTTGAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAc                                              <  2:397831/142‑1 (MQ=255)
                                                                              aaacacccGTTATATTTTTTGTTAAACGGGGCGGTGAATTGTTTAGTTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAc                                              <  2:301361/145‑1 (MQ=255)
                                                                               aacAGCCTTTATTTTTTTTTTTAAAACGGGCGTTTAAATTTTTTTTTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACa                                             <  2:201233/149‑1 (MQ=255)
                                                                                ccAGCCGTTATTTTTATTTTTAAACCTGGCGTTTAATTTTTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACaa                                            <  1:89018/150‑1 (MQ=255)
                                                                                              ttttGTTAAACCGGGGGTTTAAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTGCCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                              <  2:219330/149‑1 (MQ=255)
                                                                                               tttttttAAACCTGGCGGTAAAGTTTTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                              <  1:187640/147‑1 (MQ=255)
                                                                                                  gTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACCGCGCCAGTGCCAACGTTTTATTATCATCCGTACGCTCGTTTACCGCCGGGCCGGAGAAAAACTGCTTTAACAATTCTGGTAATAAATCGAg                          >  1:381247/1‑151 (MQ=255)
                                                                                                    ttAAACCTGGGGTTTAATTTTTTTTTTCCAACCGCCATGAGAAACTCGCTCTGCCCTCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgt                         <  2:129584/151‑1 (MQ=255)
                                                                                                                             gtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAg  >  1:288675/1‑151 (MQ=255)
                                                                                                                             |                                                                                                                                                       
GCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAG  >  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: