Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,056,018 T→A 84.6% F138I (TTC→ATC)  yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,056,0180TA84.6% 1003.2 / 158.6 350F138I (TTC→ATC) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base T (22/32);  new base A (161/135);  total (183/167)
Fisher's exact test for biased strand distribution p-value = 7.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 365 total aligned reads displayed.

TTATTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTT  >  CP009273/4055870‑4056166
                                                                                                                                                    |                                                                                                                                                    
 tATTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCttcttc                                                                                                                                                    <  2:710362/150‑1 (MQ=255)
 tATTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCttcttc                                                                                                                                                    <  2:3559647/150‑1 (MQ=255)
  aTTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCg                                                                                                                                                   <  2:2209346/150‑1 (MQ=255)
  aTTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTGACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCg                                                                                                                                                   <  1:5243772/150‑1 (MQ=255)
     tGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGaaa                                                                                                                                                >  1:4221722/1‑150 (MQ=255)
     tGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGaaa                                                                                                                                                <  1:2138688/150‑1 (MQ=255)
      gacgacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAAc                                                                                                                                               >  1:2296200/1‑150 (MQ=255)
      gacgacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAAc                                                                                                                                               >  2:4053883/1‑150 (MQ=255)
        cgacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAg                                                                                                                                             <  1:3648944/150‑1 (MQ=255)
         gacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGTGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAgg                                                                                                                                            >  2:82120/1‑150 (MQ=255)
         gacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAgg                                                                                                                                            <  1:3317314/150‑1 (MQ=255)
          acAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAggg                                                                                                                                           <  2:3351809/150‑1 (MQ=255)
             aTGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGt                                                                                                                                        >  2:4900090/1‑150 (MQ=255)
             aTGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGt                                                                                                                                        >  1:3036050/1‑150 (MQ=255)
               gCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTgg                                                                                                                                      <  2:3213644/150‑1 (MQ=255)
                ccGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGc                                                                                                                                     <  2:5441479/150‑1 (MQ=255)
                   gggTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCAtt                                                                                                                                  >  2:4605623/1‑150 (MQ=255)
                   gggTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCAtt                                                                                                                                  <  1:1212023/150‑1 (MQ=255)
                      tATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTc                                                                                                                               <  2:3016029/150‑1 (MQ=255)
                       aTTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCg                                                                                                                              >  2:5641703/1‑150 (MQ=255)
                         tGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGAt                                                                                                                            <  1:3255161/150‑1 (MQ=255)
                         tGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGAt                                                                                                                            <  2:326133/150‑1 (MQ=255)
                           gggCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATAc                                                                                                                          <  2:5058405/150‑1 (MQ=255)
                                 ccTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCgtggt                                                                                                                    >  2:1922798/1‑150 (MQ=255)
                                  ctctTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCgtggtg                                                                                                                   >  2:934242/1‑150 (MQ=255)
                                  ctctTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCgtggtg                                                                                                                   >  1:4329911/1‑150 (MQ=255)
                                  ctctTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCAGTGTCGATACCgtggtg                                                                                                                   >  2:5465874/1‑150 (MQ=255)
                                    ctTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGc                                                                                                                 >  1:1413602/1‑150 (MQ=255)
                                     tttGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCa                                                                                                                <  2:5317680/150‑1 (MQ=255)
                                      ttGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCAt                                                                                                               <  2:1221662/150‑1 (MQ=255)
                                        gTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCAttt                                                                                                             >  1:4714384/1‑150 (MQ=255)
                                          cTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTc                                                                                                           <  1:1899759/150‑1 (MQ=255)
                                           tGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCt                                                                                                          <  1:293397/150‑1 (MQ=255)
                                           tGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCt                                                                                                          >  1:2265212/1‑150 (MQ=255)
                                            gggCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTg                                                                                                         >  2:3670311/1‑150 (MQ=255)
                                                 cTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTg                                                                                                    <  2:4314069/150‑1 (MQ=255)
                                                 cTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCAGTTTCTGCGTTg                                                                                                    >  2:2895562/1‑150 (MQ=255)
                                                   aTTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGc                                                                                                  >  2:3786265/1‑150 (MQ=255)
                                                   aTTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGc                                                                                                  >  1:1794346/1‑150 (MQ=255)
                                                       ggCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTa                                                                                              >  2:2670392/1‑150 (MQ=255)
                                                        gCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTAt                                                                                             <  2:2469209/150‑1 (MQ=255)
                                                         cGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTAtt                                                                                            >  2:4561706/1‑150 (MQ=255)
                                                         cGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTAtt                                                                                            >  2:5599911/1‑150 (MQ=255)
                                                          gTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTAttt                                                                                           <  1:3997286/150‑1 (MQ=255)
                                                              cTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCt                                                                                       >  1:4888848/1‑150 (MQ=255)
                                                                  cGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGtt                                                                                   >  1:4512931/1‑150 (MQ=255)
                                                                  cGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGtt                                                                                   >  1:887692/1‑150 (MQ=255)
                                                                   gACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGttt                                                                                  <  1:2743571/150‑1 (MQ=255)
                                                                   gACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGttt                                                                                  >  1:1013081/1‑150 (MQ=255)
                                                                       tgatgaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCa                                                                              <  2:4485462/150‑1 (MQ=255)
                                                                           gaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTgg                                                                          <  1:4643316/150‑1 (MQ=255)
                                                                           gaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTgg                                                                          <  1:4424547/150‑1 (MQ=255)
                                                                           gaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTgg                                                                          <  2:4356116/150‑1 (MQ=255)
                                                                           gaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTgg                                                                          >  2:4977396/1‑150 (MQ=255)
                                                                              cTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTCTTTACCTGGTTTGGCAGTGGCtt                                                                       >  1:5566519/1‑150 (MQ=255)
                                                                               tGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCgta                                                                      >  2:2858291/1‑150 (MQ=255)
                                                                                 gCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCttatt                                                                    <  1:219572/150‑1 (MQ=255)
                                                                                   tggTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAg                                                                  >  1:1500269/2‑150 (MQ=255)
                                                                                    ggTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAgg                                                                 <  2:2017635/150‑1 (MQ=255)
                                                                                     gTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGt                                                                >  2:3434129/1‑150 (MQ=255)
                                                                                        ttGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGtttt                                                             <  2:2383698/150‑1 (MQ=255)
                                                                                        ttGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGtttt                                                             <  2:2834649/150‑1 (MQ=255)
                                                                                          gCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTaa                                                           <  2:5131419/150‑1 (MQ=255)
                                                                                              aaGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCa                                                       >  1:5313612/1‑150 (MQ=255)
                                                                                               aGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAg                                                      <  2:5320414/150‑1 (MQ=255)
                                                                                                  cGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCa                                                   <  2:1822696/150‑1 (MQ=255)
                                                                                                   ggTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCAt                                                  <  2:1708207/150‑1 (MQ=255)
                                                                                                   ggTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCAt                                                  >  2:1073239/1‑150 (MQ=255)
                                                                                                    gTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATc                                                 <  1:413356/150‑1 (MQ=255)
                                                                                                    gTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATc                                                 <  1:1624943/150‑1 (MQ=255)
                                                                                                     tAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCt                                                >  1:634114/1‑150 (MQ=255)
                                                                                                      aaGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTg                                               <  1:2417814/150‑1 (MQ=255)
                                                                                                      aaGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTg                                               <  1:342384/150‑1 (MQ=255)
                                                                                                         tCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGtt                                            <  2:3073635/150‑1 (MQ=255)
                                                                                                          cGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTc                                           >  1:374043/1‑150 (MQ=255)
                                                                                                           gATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCt                                          >  1:1839706/1‑150 (MQ=255)
                                                                                                           gATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCt                                          <  1:3471241/150‑1 (MQ=255)
                                                                                                            aTTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTa                                         <  2:1526029/150‑1 (MQ=255)
                                                                                                             ttGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTAt                                        <  1:208405/150‑1 (MQ=255)
                                                                                                                  cAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCt                                   <  1:3900126/150‑1 (MQ=255)
                                                                                                                    aTAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTa                                 >  2:2528389/1‑150 (MQ=255)
                                                                                                                    aTAGCGAACAACAAGGTCGGTGGTTTGGCTTCATCGAACCAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTa                                 <  1:3523328/150‑1 (MQ=255)
                                                                                                                     tAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTAt                                >  1:3950953/1‑150 (MQ=255)
                                                                                                                       gCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATtg                              >  1:5483501/1‑150 (MQ=255)
                                                                                                                           acaacaAGGTCGGGTGTGTGGCTTCTTCGAAACCGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATGTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTAGTCCCTTATTGTGAt                          <  2:3870881/150‑1 (MQ=255)
                                                                                                                              acaaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGc                       >  1:3400908/1‑150 (MQ=255)
                                                                                                                               caaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGcc                      >  1:2699391/1‑150 (MQ=255)
                                                                                                                               caaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGcc                      <  1:4612213/150‑1 (MQ=255)
                                                                                                                               caaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGcc                      >  2:1338928/1‑150 (MQ=255)
                                                                                                                                  ggTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTa                   >  2:4728617/1‑150 (MQ=255)
                                                                                                                                   gTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAg                  <  2:858540/150‑1 (MQ=255)
                                                                                                                                   gTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAg                  >  1:4100089/1‑150 (MQ=255)
                                                                                                                                   gTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAg                  <  1:4627458/150‑1 (MQ=255)
                                                                                                                                    tCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAgg                 >  2:5502812/1‑150 (MQ=255)
                                                                                                                                         tgttTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGAAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCatta            >  2:1483882/1‑150 (MQ=255)
                                                                                                                                           tttGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCattatt          <  1:5208494/150‑1 (MQ=255)
                                                                                                                                            ttGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCattatta         >  2:1656387/1‑150 (MQ=255)
                                                                                                                                                cttcatCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTAtttt     >  2:1973562/1‑150 (MQ=255)
                                                                                                                                                 ttcttcGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTc    <  1:2610497/150‑1 (MQ=255)
                                                                                                                                                 ttcatcGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTCGGCCTGATTATTTTc    >  1:3860229/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTATTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTT  >  CP009273/4055870‑4056166

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: