Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9170AC42.9% 164.1 / 18.9 93V43G (GTG→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (47/5);  new base C (0/39);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.32e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993767‑994064
                                                                                                                                                           |                                                                                                                                                   
ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCcca                                                                                                                                                     <  2:372944/151‑1 (MQ=255)
     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCaccacc                                                                                                                                                <  1:1085298/151‑1 (MQ=255)
      aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCg                                                                                                                                               >  1:1010221/1‑151 (MQ=255)
       gACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGc                                                                                                                                              >  2:381734/1‑151 (MQ=255)
        aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGcc                                                                                                                                             <  1:1648997/151‑1 (MQ=255)
         cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGcc                                                                                                                                             >  2:1278752/1‑150 (MQ=255)
          gcgtaaccagAAACATCATTCGCGTGTCTTCCGGAATTGCAGCCAGCGGTGTGGTGCCCGCTAACACACAGCCTGCGGTGGGCGTTACCAG‑ACCTG‑CCAGCAGGCGCCGCGGGGTTCTTTT‑G‑CCACCTCC‑GCTCCGGCCCACCGCCGCca                                                                                                                                            <  2:543763/141‑1 (MQ=255)
            gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCcaca                                                                                                                                          >  1:236944/1‑150 (MQ=255)
                 cTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                     >  1:1504217/1‑150 (MQ=255)
                      aCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                               >  1:451531/1‑151 (MQ=255)
                         tcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                            >  2:1467448/1‑151 (MQ=255)
                             tCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                         >  1:1023112/1‑150 (MQ=255)
                               ggttgtCTTCCTGATTTCCGCCCAGGGGGGTGGTCCCCGCAAACCCTCCGCCGGGGGTGGGCTTTCCCAG‑CCCGG‑CGAGCAGCCGCGGGGGGTTTTTTTT‑C‑CCCCCACC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                        <  2:497321/146‑1 (MQ=255)
                               gcgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                      >  1:141403/1‑151 (MQ=255)
                                         cTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                            >  2:1483258/1‑151 (MQ=255)
                                          tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                           >  1:268956/1‑151 (MQ=255)
                                          tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                           >  1:1583745/1‑151 (MQ=255)
                                            aaTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                         >  1:363144/1‑151 (MQ=255)
                                              tttcgcccaggGGTGGGGCCCCCCCTAACCCACCCCGTGCGGTTGGTTTTTCCAG‑CCCCG‑CCAGCGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                         <  1:1404455/143‑1 (MQ=255)
                                              ttccacccgggggTGGGGTCCCCGATAACCCCCCGCCGGGGTTGGGCTTTCCCGG‑CCCTC‑CCAGCGGGCGCAGGGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                           <  2:1211714/136‑1 (MQ=255)
                                                    ccgggggTGGGGTCCCCCCAACCCCACCGCCGGGGGTTGGCTTTCCCAG‑CCCTG‑CCAGCAGCCGCAGCGGGGTTTTTTT‑C‑CCACCCCC‑GCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                    <  1:955520/143‑1 (MQ=255)
                                                      agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                >  1:818009/1‑150 (MQ=255)
                                                       gggggggggtcCCCGCTAACCCACCGCCTGCGGTTGGCTTTTCCAG‑CCCCG‑CCAGCGGGCGCAGGGGGGTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                <  1:898666/142‑1 (MQ=255)
                                                         gttgtgGTCCCCCCAAACCCACCCCCTG‑GGTTGGGCTTTCCCG‑GCCCGCCCAGCGGGCGCGGCGGGTTTCTTTT‑G‑CCCCCACC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGata                                                                                             <  2:977494/148‑1 (MQ=255)
                                                         ggTGGGGCGCCCTCAAACAC‑GCGCCCGGGGTTTGGCTTTCCAGCCCCCG‑CCAGGGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                            <  1:1682457/151‑1 (MQ=255)
                                                          gtgtGGTCCCCGCTAACCCACCGCCGGGGTTTGGCTTTCCCAG‑CCCGC‑CCGGCGGGCGCGGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                           <  1:137687/151‑1 (MQ=255)
                                                            gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                         >  1:526859/1‑151 (MQ=255)
                                                              ggtcCCCCCAAACCCACCGCCTGGGGTTGGCTTTCCCAG‑CCCTG‑CCAGCGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                       <  1:1490134/151‑1 (MQ=255)
                                                               gtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                      >  2:1558699/1‑151 (MQ=255)
                                                                  cccccaaccccttCGCCTGGGGTGGGCTTTCCCCG‑CCCGG‑CCGGCGGGCCAAGC‑GGGTTTTTTT‑C‑CCCCCCCCGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                       <  2:1147888/135‑1 (MQ=255)
                                                                   ccGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                   >  1:1624704/1‑150 (MQ=255)
                                                                        accccaccGCCTGGGGTTGGCTTTCCCAG‑CCCTG‑CCAGCAGGCGCAGCAGGTTCTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               <  1:701954/142‑1 (MQ=255)
                                                                         acacATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                            >  1:1657392/1‑151 (MQ=255)
                                                                          cccgtCGCCGGCTGTGGGCTTTCCCGG‑CCCGG‑CCGGCAGCCGCGGCGGGGTACTTTT‑G‑CCCCCTCC‑GGCGCGCCCCTCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                            <  2:1354270/146‑1 (MQ=255)
                                                                              tCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                        >  2:298328/1‑150 (MQ=255)
                                                                               cGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      >  1:1354574/1‑151 (MQ=255)
                                                                                 ccgggggTTGGCTTTCCCAG‑CCCCC‑CCGGCGGGCGCAGCGGGGTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                       <  1:1655629/143‑1 (MQ=255)
                                                                                 ccTGGGGTTGGCTTTCCCAA‑CCCGG‑CCAGCGGGCGCGGCGGGGTCCTTTG‑C‑CCCCCCCC‑GCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                    <  1:1686237/151‑1 (MQ=255)
                                                                                 ccTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                    >  1:1460661/1‑151 (MQ=255)
                                                                                    gggttgggcttccccAG‑CCCTC‑CCAGCGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                        <  2:463629/131‑1 (MQ=255)
                                                                                      ggTTGGCTTTTCCCG‑CCCGC‑CCAGCGGGCGCAGC‑GGGTTCTTTTCC‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                                  <  1:381751/148‑1 (MQ=255)
                                                                                       gtgggCTTTTCCAG‑CCCCG‑CCAGCGGGCGCAGC‑GGGTTTTTTT‑GCCCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                    <  1:1737952/142‑1 (MQ=255)
                                                                                        ttGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                               >  1:1335189/1‑149 (MQ=255)
                                                                                        ttGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                               >  1:1581758/1‑149 (MQ=255)
                                                                                         tGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            >  2:1041039/1‑151 (MQ=255)
                                                                                          ggCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            >  1:684316/1‑150 (MQ=255)
                                                                                            cGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                         >  1:91682/1‑151 (MQ=255)
                                                                                              ttcccAG‑CCCTC‑CCAGCAGGCGCGGCGGGGTATTTTC‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                           <  1:1526524/144‑1 (MQ=255)
                                                                                              ttcccAG‑CCCGG‑CCGGCGGGCGGAGCGGGGTTTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                          <  1:909416/145‑1 (MQ=255)
                                                                                               ttcccaccccTC‑CCAGCAGGCGCGGC‑GGGTTCTTTTCC‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                          <  1:690024/140‑1 (MQ=255)
                                                                                                ttcccagcccccccAGCGGGCGCGGCGGGGTTCTTTC‑G‑CCCCCCCC‑GCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                              <  2:947453/132‑1 (MQ=255)
                                                                                                tccccaccccccccAGCAGGCGAAGCAGGGTTTTTTC‑C‑CCCCCCCC‑GCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                         <  2:48047/137‑1 (MQ=255)
                                                                                                cccAG‑CCCTG‑CCAGCGGGCGC‑GGGGGGTTCTTTTCC‑CCCCCCCC‑GCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                              <  2:904175/141‑1 (MQ=255)
                                                                                                 tccccgccccgccgccagGCACAGGAGGGTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                       <  2:777929/135‑1 (MQ=255)
                                                                                                  ccccccccagcgggcagGCCCAGCGTGTTTCTTTC‑C‑CCACCGCC‑GCCGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                      <  2:1018088/137‑1 (MQ=255)
                                                                                                  cccagcccgcccAGCGGGCGCAGCGGGGTTCTTTC‑C‑CCCCCCCC‑GCGCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                        <  1:310103/138‑1 (MQ=255)
                                                                                                   aG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                  >  2:649592/1‑151 (MQ=255)
                                                                                                    g‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                                 >  2:58631/1‑151 (MQ=255)
                                                                                                      ccccgcccggcggccgcAGCAGGTTCTTTTC‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                   <  2:114815/135‑1 (MQ=255)
                                                                                                      aCCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                                >  2:1160115/1‑151 (MQ=255)
                                                                                                        ccgcccggcgggCGCGGCGGGGTATTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                   <  1:1586535/137‑1 (MQ=255)
                                                                                                        ccccccAGCGGGCGCAGCGGGTTTTTTTC‑C‑CCCCCCCC‑GCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                                 <  1:22537/145‑1 (MQ=255)
                                                                                                        cTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                              >  2:924514/1‑151 (MQ=255)
                                                                                                         tG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGCATATGTAAATCCAGTTGGATCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                             >  2:208300/1‑151 (MQ=255)
                                                                                                            ccAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                           >  2:1166952/1‑151 (MQ=255)
                                                                                                               gcggccgcAGAGGGGTTCTTTG‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                         <  1:41627/145‑1 (MQ=255)
                                                                                                                  ggCGCAGCAGGTTACTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                           <  1:660494/145‑1 (MQ=255)
                                                                                                                  ggCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  2:750130/1‑149 (MQ=255)
                                                                                                                   gcggAGCAGGTTTTTTTC‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                            <  2:68382/143‑1 (MQ=255)
                                                                                                                    cgcAGCAGGGTATTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       <  1:1060672/147‑1 (MQ=255)
                                                                                                                    cgcAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                   >  1:510596/1‑151 (MQ=255)
                                                                                                                      cagcagGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                  >  2:712059/1‑150 (MQ=255)
                                                                                                                        gcagGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                               >  2:625114/1‑151 (MQ=255)
                                                                                                                          agGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                               >  1:1243525/1‑149 (MQ=255)
                                                                                                                           ggtttcTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                              <  1:1168309/144‑1 (MQ=255)
                                                                                                                             gTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                          >  1:771013/1‑151 (MQ=255)
                                                                                                                              tACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGGAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                         >  1:1711131/1‑151 (MQ=255)
                                                                                                                                  ttt‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                      >  1:1618942/1‑150 (MQ=255)
                                                                                                                                  ttt‑C‑CCCCCACC‑GGTGCGGCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                         <  1:1079013/147‑1 (MQ=255)
                                                                                                                                   tttcccccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        <  2:327208/138‑1 (MQ=255)
                                                                                                                                   tttcccccccccc‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        <  1:1016942/138‑1 (MQ=255)
                                                                                                                                    ttcccccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                      <  1:475920/140‑1 (MQ=255)
                                                                                                                                    t‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                   >  1:614399/1‑151 (MQ=255)
                                                                                                                                        cccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                 <  2:1605681/145‑1 (MQ=255)
                                                                                                                                        cccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                 <  1:901109/145‑1 (MQ=255)
                                                                                                                                        ccaccaCC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                 >  1:422816/1‑151 (MQ=255)
                                                                                                                                         cccccgccgcgtcGGCCCTCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                       <  2:360963/133‑1 (MQ=255)
                                                                                                                                         ccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  2:1749707/142‑1 (MQ=255)
                                                                                                                                           ccccc‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc              <  1:814547/148‑1 (MQ=255)
                                                                                                                                            cccc‑GCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc             <  2:302182/149‑1 (MQ=255)
                                                                                                                                            cacc‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc             >  1:1591389/1‑151 (MQ=255)
                                                                                                                                            cacc‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc             >  1:1693812/1‑151 (MQ=255)
                                                                                                                                             acc‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct            >  1:1719641/1‑151 (MQ=255)
                                                                                                                                                    gCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc       >  2:602701/1‑150 (MQ=255)
                                                                                                                                                     cgcccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa      <  1:529569/147‑1 (MQ=255)
                                                                                                                                                        ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt        <  2:406864/145‑1 (MQ=255)
                                                                                                                                                        cccACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  >  1:1284527/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                   
GGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993767‑994064

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