Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913826,9540AT39.7% 104.8 / 29.1 69I100N (ATT→AAT) ybhSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (36/5);  new base T (0/27);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 2.96e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCTTG  >  NC_000913/826814‑827097
                                                                                                                                               |                                                                                                                                                
gccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTcc                                                                                                                                            <  1:476073/147‑1 (MQ=255)
         cTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCAtttt                                                                                                                              >  1:823946/1‑151 (MQ=255)
              aCATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCg                                                                                                                          >  1:116422/1‑150 (MQ=255)
               cATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGAt                                                                                                                        >  1:554903/1‑151 (MQ=255)
                aTACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATc                                                                                                                       >  2:1091641/1‑151 (MQ=255)
                 tcccccTGTAAAAGGTTAGCGGTTTTCGGCTCACTGCGGTGGGTGATCCCCTGAATCGGTGCGGTGGCGTGGGCGCGCTCCATCTTTC‑CCGCAAA‑ACCCCCC‑GGAATACCCACCAGACCGCGATTTTTCC‑CCGCCTGCATTTTGGCGATc                                                                                                                       <  1:847302/148‑1 (MQ=255)
                 tcccccTGTAAAAAGTATGCGGTTTCCGCCTCACCGCGGTCGTTGACCCCCTAAATCGGTGGGGTGGGGTGGGCGCGCTCCTTTTGTC‑CCGCAAA‑ATCCCCC‑GGAATACCCCCCAGCCCGAAATTTTCCC‑CCGCCTGCATTTTGGCGAt                                                                                                                        <  1:897525/147‑1 (MQ=255)
                 tACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATc                                                                                                                       >  2:994168/1‑150 (MQ=255)
                      cTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTcc                                                                                                                 >  1:1231397/1‑151 (MQ=255)
                          aaaaaGTTGGGGGTATTCGGCTCCCTCCGGTGGGTGATCCCCTGAATGGGTGGGGTGGCTTGGGCGCCCTCCTTCTGTC‑CCGAAAA‑ACCCCCG‑GAAATACCCACCAGCCCGGGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCt                                                                                                                <  1:12679/146‑1 (MQ=255)
                          aaaaaGTTGGCGGTGTTCGGCTCCCTGGCGGTGATCATCCCCTGAATCGGTGCGGTGGTGTGGGCCCGCTCCATCGGTC‑CCCCAAA‑ACCCCCC‑GCAATACCCTCCAGGCCGCAAGTTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTg                                                                                                               <  1:619449/147‑1 (MQ=255)
                                ttAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTAt                                                                                                       >  1:549841/1‑151 (MQ=255)
                                     ggTATTCGCCTCTGCCCCGTCGTTGACCCCCTAAACCGGTGGGGGGGCTTGGGCGCCCTCCTTCTTTT‑CCACAAA‑CCCCCCC‑AAAATACCCCCCGGCCCGCTATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGc                                                                                                    <  1:628266/149‑1 (MQ=255)
                                       tATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAt                                                                                                >  2:317902/1‑151 (MQ=255)
                                           cggcccccgcccGGCGGTGACCCCCTGAATCGGGGGGGTGGCGTGGGCGCGCCCCTTCTGTC‑CCAAAAA‑ACCCACG‑GAAATACCCACCAGACCGCGATTTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtg                                                                                            <  1:776991/141‑1 (MQ=255)
                                              cTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCg                                                                                         >  1:954466/1‑151 (MQ=255)
                                                      cGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  1:1122734/1‑150 (MQ=255)
                                                      cGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  2:1174745/1‑150 (MQ=255)
                                                       gTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  1:256360/1‑149 (MQ=255)
                                                       gTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  1:1316187/1‑149 (MQ=255)
                                                         cGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAggg                                                                              >  1:410720/1‑151 (MQ=255)
                                                            tGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCATCGATGTAGGGCg                                                                            >  2:798762/1‑150 (MQ=255)
                                                             gATCCCCTGAATGGGGGCGGTGGCTTGGGCCCGCT‑CTTTTGTTCCCCCAAA‑ATCCCCG‑GAAATAAACCCCGGGCAGCAAATTTGCC‑CCGCCTGCATTTTGGCGGTCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGa                                                                           <  2:688105/150‑1 (MQ=255)
                                                             gATCCCCGGAATCGGGGCGGTGGCGTTGGCGCGCTCCATCTGTC‑CCGAAAA‑ACCCCCC‑GAAATAACCCCCAGCCCGCAAATTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGaa                                                                          <  1:420603/151‑1 (MQ=255)
                                                               tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAcc                                                                        >  1:1211475/1‑151 (MQ=255)
                                                                ccccTGATTCGGTGGGTGGGCGTGGCCGCGCTCAATCTGTT‑CCGCAAA‑ACCCCAC‑GGAAAACCCACAAGCCCGGGATTTTCCC‑CCGCCTGCATTTTGGCGAGCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGc                                                                             <  2:5938/143‑1 (MQ=255)
                                                                 aCCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCgg                                                                      >  1:840818/1‑151 (MQ=255)
                                                                  cctaacccggggcggggGGGTGGGCGCCCTCCATCTGTC‑CCGAAAA‑ACCCCCC‑GAAATACCCCCCAGACCGAGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGaa                                                                          <  2:1231333/139‑1 (MQ=255)
                                                                    tGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTc                                                                    >  2:1168568/1‑150 (MQ=255)
                                                                    tGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTc                                                                    >  2:998534/1‑150 (MQ=255)
                                                                      actggttggtggggCGTTGGCGCGCCCCATCTTTC‑CCGCAAA‑ATCCACC‑GGAATACCCCCCAGACCGCTATTTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                         <  1:89998/131‑1 (MQ=255)
                                                                      aaggggtggggggGCGTGGGCCCGCTCCATCTGTC‑CCAAAAA‑CCCCCCC‑GGAATACCCCCCAGCCCGAGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGt                                                                     <  2:1101848/143‑1 (MQ=255)
                                                                             ggggtgggtggggCGCGCTCCTTTTGTC‑CCGCAAA‑TCCCCCC‑GGATAAACCACC‑GCCCCCTGATTTCCCACCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAgggtg                                                                            <  2:250315/122‑3 (MQ=255)
                                                                             gcggtgGCGTTGGCGCCCTCCATCTTTC‑CCGAAAA‑ACCCCCC‑GAAATACCCCCCAGCCCGCGATTTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGtgtg                                                            <  1:395551/149‑1 (MQ=255)
                                                                               ggtggtGTTGGCGCCCTCCATCTGTT‑CCAAAAA‑ACCCACG‑GAAATACCCACCAGCCCGGGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTgg                                                           <  1:914265/148‑1 (MQ=255)
                                                                               ggtgGCGTTGGCGCGCTCCATCTGTT‑CCGCAGA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTg                                                        >  2:1135618/1‑151 (MQ=255)
                                                                               ggtgGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGt                                                         >  2:936813/1‑150 (MQ=255)
                                                                               ggtgGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTg                                                        >  2:1216255/1‑151 (MQ=255)
                                                                                     gTTGGCGCGCTCTATCTGTA‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATc                                                   >  2:329466/1‑150 (MQ=255)
                                                                                        ggcccccTCCTTCTGTC‑CCAAAAA‑ACCCCCC‑GAAATACCCCCCGGCCCGGAAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATc                                                   <  1:360621/141‑1 (MQ=255)
                                                                                        ggCGCGCCCCACCTTTC‑CCGCAAA‑CGCCCCC‑GAAATAACCACCCGGCCGCTATTTTTCG‑CCGCCTGCATTTTGGCGATCAGTTCCTGGCGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCa                                                 <  2:453793/149‑1 (MQ=255)
                                                                                        gcctccccccctttGTC‑CCGAAAA‑CCCCCCC‑GAAAAACCCCCCAGACCCCGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAg                                                <  2:903098/137‑1 (MQ=255)
                                                                                            cgcTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGcgccg                                           >  2:538180/1‑151 (MQ=255)
                                                                                             gcTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGcgccgc                                          >  1:1114343/1‑151 (MQ=255)
                                                                                                ccATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTc                                       >  2:808812/1‑151 (MQ=255)
                                                                                                 cttCTGTC‑CCGAAAA‑ACCCCCG‑GAAATACCCCCCAGCCCGCGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCtt                                        <  1:856325/147‑1 (MQ=255)
                                                                                                 cATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTc                                       >  1:74747/1‑150 (MQ=255)
                                                                                                  aTCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGc                                     >  2:1272684/1‑151 (MQ=255)
                                                                                                     tctttgcGCAAAGATCCCCC‑GTAATAACCACAAGCCGGCGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAgc                                               <  1:284511/134‑1 (MQ=255)
                                                                                                              aaa‑aTCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTcc                          >  1:919098/1‑151 (MQ=255)
                                                                                                              aaa‑aTCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTc                           >  1:1238418/1‑150 (MQ=255)
                                                                                                               aa‑aTCCCCCGGGAATACCCACCACCCCGAGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGc                                <  2:1075575/145‑1 (MQ=255)
                                                                                                                     ccccagaaaaaccccccaGCCCGCGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTaa                     <  2:325540/136‑1 (MQ=255)
                                                                                                                     cACC‑GAAATAACCCCCAGCCCGCGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTaaa                    <  2:116418/151‑1 (MQ=255)
                                                                                                                       cc‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAAt                  >  2:967460/1‑151 (MQ=255)
                                                                                                                       cc‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAAt                  >  1:908156/1‑151 (MQ=255)
                                                                                                                          gcaaTACCCCCCAGACCGCGATTTTTCT‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATcc                <  2:1008556/149‑1 (MQ=255)
                                                                                                                            aaTACCCCCCAGCCCGCATTTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATc                 <  2:948477/148‑1 (MQ=255)
                                                                                                                              tAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGAc            >  2:956493/1‑151 (MQ=255)
                                                                                                                                 cccccaGCCCGGGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGa             <  1:371404/144‑1 (MQ=255)
                                                                                                                                   acccGACCGCGAATTTTCC‑CCGCAGGCATTTTTGCGATCAGTTCCTGACGGTTATCGCTGTTAGTGGTGACGATGTAGGTTGAACCGGTCATGGGGTGGGTGATATCCACCGCCGCTTCGCTACGCTGTTCCAGCAAACTCCCGAACCGc        >  2:1107989/1‑150 (MQ=255)
                                                                                                                                   accaGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGc        >  1:333840/1‑150 (MQ=255)
                                                                                                                                      aGACCGGGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCa       <  1:1237019/148‑1 (MQ=255)
                                                                                                                                       gACCGCGATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAg      <  1:1179052/148‑1 (MQ=255)
                                                                                                                                       gACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCtt   >  2:1061086/1‑151 (MQ=255)
                                                                                                                                         ccGCTATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCTTg  <  2:1272626/150‑1 (MQ=255)
                                                                                                                                         ccGCAATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCa       <  2:763210/145‑1 (MQ=255)
                                                                                                                                         ccGCAATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAg      <  2:1041259/146‑1 (MQ=255)
                                                                                                                                         ccGCAATTTTCCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCtt   <  1:70541/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTT‑CCGCAAA‑ATCCACC‑GGAATAACCACCAGACCGCGAATTTTCC‑CCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCTTG  >  NC_000913/826814‑827097

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: ‑