Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 435,336 G→T 100% Q100K (CAG→AAG)  dxs ← 1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273435,3360GT95.6% 974.2 / 30.2 293Q100K (CAG→AAG) dxs1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring
Reads supporting (aligned to +/- strand):  ref base G (5/8);  new base T (146/134);  total (151/142)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 295 total aligned reads displayed.

GACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCAGCCCGGAGG  >  CP009273/435187‑435482
                                                                                                                                                     |                                                                                                                                                  
     aGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGa                                                                                                                                               <  1:2475118/150‑1 (MQ=255)
      gACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGAt                                                                                                                                              >  1:2225746/1‑150 (MQ=255)
          gTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTg                                                                                                                                          <  1:509421/150‑1 (MQ=255)
          gTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTg                                                                                                                                          <  1:3519009/150‑1 (MQ=255)
           tGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGc                                                                                                                                         <  1:2657323/150‑1 (MQ=255)
                cgATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGAtt                                                                                                                                    <  1:3678654/150‑1 (MQ=255)
                 gATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAATGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGAttt                                                                                                                                   <  1:2564544/150‑1 (MQ=255)
                     tttGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTc                                                                                                                               <  1:83404/150‑1 (MQ=255)
                       tGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTcgc                                                                                                                             >  2:417377/1‑150 (MQ=255)
                        gCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTcgcg                                                                                                                            <  2:2729072/150‑1 (MQ=255)
                        gCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTcgcg                                                                                                                            >  1:1519700/1‑150 (MQ=255)
                        gCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTcgcg                                                                                                                            <  2:1534456/150‑1 (MQ=255)
                         ccTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTGTCTTACGGATGGTGCCGATTTTGTCgcgg                                                                                                                           >  2:2385507/1‑150 (MQ=255)
                             cttTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTc                                                                                                                       >  2:129638/1‑150 (MQ=255)
                              tttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTcc                                                                                                                      <  2:2010641/150‑1 (MQ=255)
                                tttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCgg                                                                                                                    <  1:4101809/150‑1 (MQ=255)
                                tttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCgg                                                                                                                    >  1:1099996/1‑150 (MQ=255)
                                 ttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGt                                                                                                                   <  2:3776350/150‑1 (MQ=255)
                                  tCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTc                                                                                                                  <  1:3894864/150‑1 (MQ=255)
                                   cGGCAGCAACCGCAATACCAAGTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCa                                                                                                                 <  1:2557599/150‑1 (MQ=255)
                                    ggCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCaa                                                                                                                >  2:2627237/1‑150 (MQ=255)
                                    ggCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTGTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCaa                                                                                                                >  1:3131307/1‑150 (MQ=255)
                                      cagcaACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCaaaa                                                                                                              >  1:2968943/1‑150 (MQ=255)
                                          aaCCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAAtttt                                                                                                          >  1:2192778/1‑150 (MQ=255)
                                             cGCAATACCCATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATg                                                                                                       >  1:2192564/1‑150 (MQ=255)
                                             cGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATg                                                                                                       >  1:4596954/1‑150 (MQ=255)
                                             cGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATg                                                                                                       >  1:2145158/1‑150 (MQ=255)
                                              gCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGc                                                                                                      <  2:3532202/150‑1 (MQ=255)
                                                aaTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCgg                                                                                                    >  1:1749779/1‑150 (MQ=255)
                                                aaTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCgg                                                                                                    >  1:3827733/1‑150 (MQ=255)
                                                aaTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCgg                                                                                                    <  1:3231587/150‑1 (MQ=255)
                                                 aTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGa                                                                                                   >  1:1946172/1‑150 (MQ=255)
                                                  tACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGAt                                                                                                  <  2:4442129/150‑1 (MQ=255)
                                                    ccAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATaa                                                                                                >  1:2913725/1‑150 (MQ=255)
                                                     cAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAg                                                                                               <  2:367230/150‑1 (MQ=255)
                                                     cAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAg                                                                                               >  1:1853812/1‑150 (MQ=255)
                                                     cAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAg                                                                                               >  1:1245862/1‑150 (MQ=255)
                                                      aaTTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGc                                                                                              <  2:3350032/150‑1 (MQ=255)
                                                           cGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGACCGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGAt                                                                                         >  2:3773726/1‑150 (MQ=255)
                                                           cGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGAt                                                                                         >  1:1316601/1‑150 (MQ=255)
                                                           cGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGAt                                                                                         >  1:3726304/1‑150 (MQ=255)
                                                            ggCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATg                                                                                        <  2:1930643/150‑1 (MQ=255)
                                                               aCTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGccc                                                                                     >  1:114356/1‑150 (MQ=255)
                                                                cTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGcccc                                                                                    >  2:1479273/1‑150 (MQ=255)
                                                                cTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGcccc                                                                                    <  1:4137302/150‑1 (MQ=255)
                                                                 tGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCca                                                                                   >  2:1228753/1‑150 (MQ=255)
                                                                  gATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCcac                                                                                  <  1:1483174/150‑1 (MQ=255)
                                                                  gATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCcac                                                                                  <  1:1659292/150‑1 (MQ=255)
                                                                  gATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCcac                                                                                  >  2:1103297/1‑150 (MQ=255)
                                                                  gATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCcac                                                                                  <  1:2530813/150‑1 (MQ=255)
                                                                   aTGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCcaca                                                                                 >  2:2298791/1‑150 (MQ=255)
                                                                     ggAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATc                                                                               <  1:1880145/150‑1 (MQ=255)
                                                                             atgaATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCa                                                                       >  2:3030888/1‑150 (MQ=255)
                                                                               gaATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAAt                                                                     >  2:4567819/1‑150 (MQ=255)
                                                                                aaTGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAAtt                                                                    <  1:4420149/150‑1 (MQ=255)
                                                                                aaTGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAAtt                                                                    <  1:2015552/150‑1 (MQ=255)
                                                                                 aTGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTg                                                                   >  2:1601064/1‑150 (MQ=255)
                                                                                    cccGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTc                                                                <  1:1059119/150‑1 (MQ=255)
                                                                                     ccGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCa                                                               <  2:2753679/150‑1 (MQ=255)
                                                                                     ccGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCa                                                               <  1:3486196/150‑1 (MQ=255)
                                                                                     ccGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAAATTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCa                                                               <  1:2340904/150‑1 (MQ=255)
                                                                                      cGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCCAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCaa                                                              >  2:4000958/1‑150 (MQ=255)
                                                                                      cGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCaa                                                              >  1:3431470/1‑150 (MQ=255)
                                                                                       gACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCaaa                                                             >  1:4778930/1‑150 (MQ=255)
                                                                                         cGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCCCATCCCAAATCAATTGGTCAAACg                                                           >  1:990753/1‑150 (MQ=255)
                                                                                         cGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACg                                                           <  2:3133220/150‑1 (MQ=255)
                                                                                           cTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACggg                                                         <  1:1701222/150‑1 (MQ=255)
                                                                                            ttAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACgggg                                                        >  2:3000992/1‑150 (MQ=255)
                                                                                             tAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGgt                                                       >  1:154347/1‑150 (MQ=255)
                                                                                             tAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGgt                                                       <  2:1668597/150‑1 (MQ=255)
                                                                                               aTACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGgtgt                                                     <  2:2428943/150‑1 (MQ=255)
                                                                                               aTACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGgtgt                                                     <  1:2838256/150‑1 (MQ=255)
                                                                                                    tCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAg                                                >  2:1792149/1‑150 (MQ=255)
                                                                                                      atatTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGAc                                              >  2:2084261/1‑150 (MQ=255)
                                                                                                       tatTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACa                                             <  1:2400741/150‑1 (MQ=255)
                                                                                                       tatTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACa                                             >  2:2544149/1‑150 (MQ=255)
                                                                                                        atTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACAt                                            <  1:1988907/150‑1 (MQ=255)
                                                                                                        atTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACAt                                            >  1:1128109/1‑150 (MQ=255)
                                                                                                         ttCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATa                                           >  2:281589/1‑150 (MQ=255)
                                                                                                         ttCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATa                                           <  1:2084261/150‑1 (MQ=255)
                                                                                                          tCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAg                                          <  2:2466159/150‑1 (MQ=255)
                                                                                                                 cgccgcGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGTTGCCCCACATCCCAAATCAATTGGTCCAACGGGGTGTTGTAGACATAGTGCAgcg                                   >  1:3935326/1‑150 (MQ=255)
                                                                                                                 cgccgcGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAgcg                                   >  1:4301790/1‑150 (MQ=255)
                                                                                                                 cgccgcGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAgcg                                   >  2:4310382/1‑150 (MQ=255)
                                                                                                                 cgccgcGCCACGGGAACGGGTGCAGACCGCCTGTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTAGGCGGATAAGCCTGCTGCCCCACATCCCAAATCACTTGGTCAAACGGGGTGTTGTAGACATAGTGCAgcg                                   <  1:3812867/150‑1 (MQ=255)
                                                                                                                    cgcgCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCa                                >  1:957061/1‑150 (MQ=255)
                                                                                                                      cgcCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACg                              >  1:1720281/1‑150 (MQ=255)
                                                                                                                       gcCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACgg                             >  2:813476/1‑150 (MQ=255)
                                                                                                                         cACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTc                           >  2:3859448/1‑150 (MQ=255)
                                                                                                                         cACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTc                           >  2:1449176/1‑150 (MQ=255)
                                                                                                                           cGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAg                         >  2:2350271/1‑150 (MQ=255)
                                                                                                                            gggAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGt                        <  2:3503412/150‑1 (MQ=255)
                                                                                                                            gggAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGt                        >  2:613791/1‑150 (MQ=255)
                                                                                                                             ggAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGtt                       <  2:1590805/150‑1 (MQ=255)
                                                                                                                              gAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTc                      >  2:3222176/1‑150 (MQ=255)
                                                                                                                              gAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTc                      >  2:471701/1‑150 (MQ=255)
                                                                                                                                aCGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGGTGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGa                    <  1:2298791/150‑1 (MQ=255)
                                                                                                                                       cAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGca             >  1:72123/1‑149 (MQ=255)
                                                                                                                                          aCCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCAg          >  2:4143987/1‑150 (MQ=255)
                                                                                                                                               cTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCAGCCCgg     >  1:3755004/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GACACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCAGCCCGGAGG  >  CP009273/435187‑435482

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