Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,826,247 A→T 94.8% Y551N (TAC→AAC)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,826,2470AT94.8% 679.3 / 20.2 221Y551N (TAC→AAC) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (6/6);  new base T (106/103);  total (112/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 221 total aligned reads displayed.

ATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTG  >  CP009273/3826104‑3826390
                                                                                                                                               |                                                                                                                                               
 tCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGccc                                                                                                                                            >  1:1242366/1‑148 (MQ=255)
    tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACg                                                                                                                                         >  2:1947165/1‑148 (MQ=255)
      cGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGc                                                                                                                                       <  2:660141/148‑1 (MQ=255)
         gAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGcaca                                                                                                                                    >  1:1698794/1‑148 (MQ=255)
         gAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGcaca                                                                                                                                    >  2:2862507/1‑148 (MQ=255)
           aCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACg                                                                                                                                  <  1:34570/148‑1 (MQ=255)
            cTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGa                                                                                                                                 <  1:35899/148‑1 (MQ=255)
             tCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGAt                                                                                                                                >  1:1787339/1‑148 (MQ=255)
               catcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATaa                                                                                                                              >  1:187395/1‑148 (MQ=255)
                 tcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGa                                                                                                                            <  1:3139492/148‑1 (MQ=255)
                 tcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGa                                                                                                                            <  1:764171/148‑1 (MQ=255)
                    tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGAttt                                                                                                                         <  1:2000040/148‑1 (MQ=255)
                     catGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                                                                                                        <  2:2958510/148‑1 (MQ=255)
                       tGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCt                                                                                                                      >  2:1986321/1‑148 (MQ=255)
                       tGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCt                                                                                                                      >  2:2606947/1‑148 (MQ=255)
                       tGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCt                                                                                                                      >  2:217215/1‑148 (MQ=255)
                        ggCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTa                                                                                                                     <  2:2329113/148‑1 (MQ=255)
                        ggCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTa                                                                                                                     >  2:1628272/1‑148 (MQ=255)
                        ggCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTa                                                                                                                     >  2:3366989/1‑148 (MQ=255)
                         gCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTAc                                                                                                                    >  1:581969/1‑148 (MQ=255)
                         gCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTAc                                                                                                                    >  1:582148/1‑148 (MQ=255)
                               atcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTa                                                                                                              >  1:2429051/1‑148 (MQ=255)
                                tcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTaa                                                                                                             >  1:1001939/1‑148 (MQ=255)
                                tcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTaa                                                                                                             >  1:2799794/1‑148 (MQ=255)
                                tcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTaa                                                                                                             >  1:2757522/1‑148 (MQ=255)
                                 catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTaaa                                                                                                            >  2:678396/1‑148 (MQ=255)
                                  atcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAAc                                                                                                           <  1:127747/148‑1 (MQ=255)
                                   tcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACg                                                                                                          <  2:82508/148‑1 (MQ=255)
                                    cGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACgg                                                                                                         >  1:2616767/1‑148 (MQ=255)
                                    cGGCGTACCCCGCGCGGTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACgg                                                                                                         <  2:2422114/148‑1 (MQ=255)
                                     ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGc                                                                                                        <  2:44654/148‑1 (MQ=255)
                                         tACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATg                                                                                                    <  1:695041/148‑1 (MQ=255)
                                          aCCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATgg                                                                                                   <  2:1453351/148‑1 (MQ=255)
                                                cgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGgag                                                                                             >  1:1487834/1‑148 (MQ=255)
                                                  cgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGagag                                                                                           >  1:973064/1‑148 (MQ=255)
                                                      cGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaa                                                                                       <  1:389592/148‑1 (MQ=255)
                                                         aCGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaacca                                                                                    >  2:3405730/1‑148 (MQ=255)
                                                          cgcgCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaaccaa                                                                                   >  2:2632312/1‑148 (MQ=255)
                                                           gcgcAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaaccaaa                                                                                  <  1:109739/148‑1 (MQ=255)
                                                           gcgcAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaaccaaa                                                                                  <  2:2756222/148‑1 (MQ=255)
                                                             gcAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAAcg                                                                                <  2:1741791/148‑1 (MQ=255)
                                                              cAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAAcgc                                                                               >  2:3183220/1‑148 (MQ=255)
                                                               aGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAAcgcg                                                                              >  2:1694401/1‑148 (MQ=255)
                                                                   tCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCAcc                                                                          >  2:206839/1‑148 (MQ=255)
                                                                      cGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGc                                                                       >  2:1247382/1‑148 (MQ=255)
                                                                      cGAGACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGc                                                                       <  2:1860212/148‑1 (MQ=255)
                                                                       gATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCg                                                                      <  1:1166456/148‑1 (MQ=255)
                                                                         tACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGtt                                                                    <  2:890744/148‑1 (MQ=255)
                                                                           cAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTg                                                                  >  2:2041571/1‑148 (MQ=255)
                                                                           cAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTg                                                                  >  2:2042383/1‑148 (MQ=255)
                                                                            aGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACCGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGt                                                                 <  2:1201321/148‑1 (MQ=255)
                                                                             gATACGGCATCATTCGGCATATCAGTTGCGTGAAGAAGCGCACCACATCACATGACTCATCAACGTTGGCCCACGGCACACGAAAAGATAAGATAACGAGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCTCACCAGCCtttcaa                                                                >  2:2340803/1‑145 (MQ=255)
                                                                                aCGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAAc                                                             <  1:2045676/148‑1 (MQ=255)
                                                                                 cGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACg                                                            <  2:925788/148‑1 (MQ=255)
                                                                                 cGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACg                                                            <  2:1849998/148‑1 (MQ=255)
                                                                                  ggCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGt                                                           <  1:930542/148‑1 (MQ=255)
                                                                                   gCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTg                                                          <  1:1913055/148‑1 (MQ=255)
                                                                                   gCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTg                                                          <  2:3282622/148‑1 (MQ=255)
                                                                                         tGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCgg                                                    <  2:2033346/148‑1 (MQ=255)
                                                                                            ggCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGc                                                 <  1:1708682/148‑1 (MQ=255)
                                                                                              cATGTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCgg                                               <  1:2095427/148‑1 (MQ=255)
                                                                                                tttCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTa                                             <  2:2696215/148‑1 (MQ=255)
                                                                                                 ttCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTAt                                            >  2:220261/1‑148 (MQ=255)
                                                                                                    aGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTAtttt                                         <  2:3401223/148‑1 (MQ=255)
                                                                                                       tGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCaa                                      >  1:1021736/1‑148 (MQ=255)
                                                                                                       tGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCaa                                      <  2:1906785/148‑1 (MQ=255)
                                                                                                       tGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCaa                                      <  2:133948/148‑1 (MQ=255)
                                                                                                       tGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCaa                                      <  1:3085900/148‑1 (MQ=255)
                                                                                                       tGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTCTCaa                                      >  2:3426393/1‑148 (MQ=255)
                                                                                                        gCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCaaa                                     <  2:3249643/148‑1 (MQ=255)
                                                                                                         cGTGAAGAAGCGCCCCACATCACAGGACTCATAATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGATATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAg                                    <  2:1597872/148‑1 (MQ=255)
                                                                                                         cGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAg                                    <  2:1699402/148‑1 (MQ=255)
                                                                                                           tGAAGAAGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAgcc                                  >  1:2256641/1‑148 (MQ=255)
                                                                                                            gaagaaGCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAgccg                                 >  2:1762961/1‑148 (MQ=255)
                                                                                                               gaagCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGccgccg                              >  2:1438086/1‑148 (MQ=255)
                                                                                                                aagCGCACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGa                             >  1:1850909/1‑148 (MQ=255)
                                                                                                                   cgcACCACATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGatat                          <  1:625346/148‑1 (MQ=255)
                                                                                                                      accacATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGt                       >  2:2210961/1‑148 (MQ=255)
                                                                                                                       ccacATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTg                      >  1:1070253/1‑148 (MQ=255)
                                                                                                                       ccacATCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTg                      >  2:3471023/1‑148 (MQ=255)
                                                                                                                        cacaTCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                     <  2:2994255/148‑1 (MQ=255)
                                                                                                                        cacaTCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                     >  1:2566322/1‑148 (MQ=255)
                                                                                                                        cacaTCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                     >  1:2321786/1‑148 (MQ=255)
                                                                                                                        cacaTCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                     >  1:334756/1‑148 (MQ=255)
                                                                                                                        cacaTCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                     <  1:2189592/148‑1 (MQ=255)
                                                                                                                        cacaTCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                     <  1:2189564/148‑1 (MQ=255)
                                                                                                                        cacaTCACAGGACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                     <  1:176999/148‑1 (MQ=255)
                                                                                                                          ccacacaGGCCACATCATCGTTGGCCCACGGCACACGATAAGATTTTCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCt                   <  2:1544578/145‑1 (MQ=255)
                                                                                                                                 ggACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCTTGTAAACGGCTATGGCTGGAGAGCATACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCGGCCGATATCGTGGCTCCAGAAg            <  2:582148/148‑1 (MQ=255)
                                                                                                                                 ggACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAg            <  1:1334359/148‑1 (MQ=255)
                                                                                                                                 ggACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAg            <  1:781856/148‑1 (MQ=255)
                                                                                                                                 ggACTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAg            <  2:581969/148‑1 (MQ=255)
                                                                                                                                  gACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGc           <  1:3183220/148‑1 (MQ=255)
                                                                                                                                   aCTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGcc          >  2:2930890/1‑148 (MQ=255)
                                                                                                                                   aCTCATCATCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGcc          >  2:3288115/1‑148 (MQ=255)
                                                                                                                                      catcatCGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCaaa       <  2:2104782/148‑1 (MQ=255)
                                                                                                                                        tcatcGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAAc     >  2:3031320/1‑148 (MQ=255)
                                                                                                                                         catcGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAAcc    <  1:1247382/148‑1 (MQ=255)
                                                                                                                                          atcGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCt   >  2:76154/1‑148 (MQ=255)
                                                                                                                                          atcGTTGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCt   >  1:606048/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTG  >  CP009273/3826104‑3826390

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: