Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,534,2340GT46.7% 33.3 / 14.0 30E71* (GAG→TAG) nhoAN‑hydroxyarylamine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/13);  new base T (14/0);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 4.79e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAA  >  NC_000913/1534090‑1534375
                                                                                                                                                 |                                                                                                                                              
gATACGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTTTTGTTTTTAg                                                                                                                                              >  1:660551/1‑147 (MQ=255)
 aTACGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTTAGCAg                                                                                                                                           >  2:335700/1‑149 (MQ=255)
        tCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTTTGCAGAATGGCggg                                                                                                                                  >  2:134936/1‑149 (MQ=255)
          gTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGtt                                                                                                                                <  1:402455/151‑1 (MQ=255)
               ttGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTTAGTAGAAAGGCGGGGTTGAg                                                                                                                            >  1:461663/1‑150 (MQ=255)
                  cACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTTCTGTTTTGGGGCGAATGGCGGGTTTGAGgggg                                                                                                                        >  1:581805/1‑151 (MQ=255)
                     cTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTTAGGAGAATGGCGGGTTTGAGCggggg                                                                                                                      >  1:307443/1‑148 (MQ=255)
                       gAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGtt                                                                                                                    <  1:694023/150‑1 (MQ=255)
                                                gAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTTCTGTTTTTTGGAGAAAGGGGGGTTTGAGCGGGGGTTTCGCGAGCT‑GGGGGTTTACGTTCg                                                                                           >  2:292221/1‑149 (MQ=255)
                                                gAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCa                                                                                         <  2:193018/151‑1 (MQ=255)
                                                  aaaCCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTTCTGTTTTTAGCAGAAAGGCGGGGTTTTGCGGGGGTTTCGCGCGCT‑GGGGGTTTAAGTTTGCCGc                                                                                       >  2:141802/1‑151 (MQ=255)
                                                     ccTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTTCTGTTTTTAGGAGAAAGGGGGGGTTGAGCGGGGGTTTCGCGGGGT‑GGGGGTTTAAGGTCGCacgct                                                                                     >  1:487708/1‑146 (MQ=255)
                                                        cGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTa                                                                                 <  2:608091/151‑1 (MQ=255)
                                                           cGTTTTGCTGCCGAGAGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTTTGGAGAAAGGGGTGGTTTTGGGGGGGTTTTGCGAGCT‑GGGGGTTTAAGTTTGCCGCTTGTTTggg                                                                              >  2:698771/1‑151 (MQ=255)
                                                            gTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTTCTCTTTTTTGGAGAAAGGGGGGTTTTAGGGGGGGGTTCGCGGGCT‑GGGGGTTAAAGTTTGCGGCCTggttgtgg                                                                             >  1:455142/1‑144 (MQ=255)
                                                            gTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTTTGCAGAAAGGGGGGGTTTAGGGGGGGTTACGCGAGGT‑GGGGGTTTAAGTTCGCCGCTTGTTTgggg                                                                             >  2:550572/1‑150 (MQ=255)
                                                               ttGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGgcgc                                                                           <  2:711877/150‑1 (MQ=255)
                                                                 gctgcCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGt                                                                         <  1:508957/150‑1 (MQ=255)
                                                                                 cAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTAATGTTTTTAGCAGAAAGGGGGGTTTGAGGGGGGGGTTTGCGAGGT‑GGGGGTTTAAGTTTGCCGCTTGTTTGGGGGCGGAGTGttttttaa                                                             >  2:28092/1‑140 (MQ=255)
                                                                                                tcGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTAccgccg                                          <  2:582083/150‑1 (MQ=255)
                                                                                                tcGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑CTTATGTTTTTAGGAGAATGGGGTGGTTTGGCGGGGGGTTAGCGAGGG‑GGGGGTTTAAAATCTCAGCGTGTTATGGGGCGGAGGGGTTTTAAAATCGCCCGCATTAacg                                             >  2:582099/1‑147 (MQ=255)
                                                                                                    cGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCAc                                     <  1:228398/151‑1 (MQ=255)
                                                                                                        agagAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCAt                                 <  2:267582/151‑1 (MQ=255)
                                                                                                          agaAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCg                               <  1:311951/151‑1 (MQ=255)
                                                                                                              aCTGGTGATAGCCCGTCGTGGCGGTTTACT‑TTTTTTGCCGAAAGGGGGGGTTGGGCGGGGGGTTCGCGGGGTGGGGGTTTAAGGTT‑GCCGCTTGTTTGGGGGCGGGGGGTTTTCAAAACCGCCCGCACTTCCCCCCCGCCCCCCTCGtt                             >  1:217330/1‑149 (MQ=255)
                                                                                                                    gATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGtt                      <  1:593637/150‑1 (MQ=255)
                                                                                                                         cccGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGACGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGa                <  1:152473/151‑1 (MQ=255)
                                                                                                                                 ggCGG‑TTACTGTTTTGAGCAGAATGGCGTTTTTGAGCGGGTGTTACGCGAGCT‑GGAGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGAAACTGGaa         <  1:693203/150‑1 (MQ=255)
                                                                                                                                    gg‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAgg     >  2:238149/1‑151 (MQ=255)
                                                                                                                                         aCTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGaaa  >  2:326603/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GATACGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGG‑TTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCT‑GGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAA  >  NC_000913/1534090‑1534375

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