Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,930,795 G→T 6.6% A48S (GCT→TCT)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,930,7950GT 6.6% 1111.5 / 34.5 288A48S (GCT→TCT) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base G (149/120);  new base T (8/11);  total (157/131)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 291 total aligned reads displayed.

CCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTG  >  CP009273/3930648‑3930944
                                                                                                                                                   |                                                                                                                                                     
ccGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCt                                                                                                                                                     <  2:853706/150‑1 (MQ=255)
 cGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCtg                                                                                                                                                    >  2:1463330/1‑150 (MQ=255)
    atatGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTgg                                                                                                                                                 >  2:2722839/1‑150 (MQ=255)
     tatGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgc                                                                                                                                                >  1:1751351/1‑150 (MQ=255)
      atGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgct                                                                                                                                               <  1:2045539/150‑1 (MQ=255)
      atGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgct                                                                                                                                               >  1:1992057/1‑150 (MQ=255)
      atGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgct                                                                                                                                               >  1:192035/1‑150 (MQ=255)
           aaaCGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTgg                                                                                                                                          >  1:4073223/1‑150 (MQ=255)
           aaaCGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTgg                                                                                                                                          <  1:435230/150‑1 (MQ=255)
           aaaCGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTgg                                                                                                                                          >  2:3203390/1‑150 (MQ=255)
              cGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCg                                                                                                                                       >  2:2369157/1‑150 (MQ=255)
               gcaggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGt                                                                                                                                      <  1:3512982/150‑1 (MQ=255)
                caggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCTAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTTGGCGTa                                                                                                                                     >  2:2066122/1‑150 (MQ=255)
                caggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTa                                                                                                                                     <  1:1727776/150‑1 (MQ=255)
                caggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTa                                                                                                                                     <  1:794439/150‑1 (MQ=255)
                caggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTa                                                                                                                                     <  1:832166/150‑1 (MQ=255)
                 aggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGTCTGTGGCTGCTGGGCGTAg                                                                                                                                    >  1:2130201/1‑150 (MQ=255)
                 aggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGTCTGTGGCTGCTGGGCGTAg                                                                                                                                    >  1:2299219/1‑150 (MQ=255)
                 aggcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAg                                                                                                                                    >  2:1227673/1‑150 (MQ=255)
                    cagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCgg                                                                                                                                 >  2:3174571/1‑150 (MQ=255)
                     agCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGt                                                                                                                                >  2:2834303/1‑150 (MQ=255)
                     agCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAAGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGt                                                                                                                                <  1:2421144/150‑1 (MQ=255)
                      gCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTg                                                                                                                               >  2:1690068/1‑150 (MQ=255)
                             tgttCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAatat                                                                                                                        <  2:2482385/150‑1 (MQ=255)
                                  ttGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGt                                                                                                                   >  1:3050438/1‑150 (MQ=255)
                                    ggCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGttt                                                                                                                 <  1:2384168/150‑1 (MQ=255)
                                     gcagcaTTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTa                                                                                                                >  1:224684/1‑150 (MQ=255)
                                     gcagcaTTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTa                                                                                                                <  2:4052675/150‑1 (MQ=255)
                                      cagcaTTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTAt                                                                                                               <  1:1009086/150‑1 (MQ=255)
                                               tGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTAc                                                                                                      <  2:1154532/150‑1 (MQ=255)
                                                 cTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACgg                                                                                                    <  2:1476126/150‑1 (MQ=255)
                                                 cTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACgg                                                                                                    >  1:112481/1‑150 (MQ=255)
                                                 cTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACgg                                                                                                    <  2:907012/150‑1 (MQ=255)
                                                   gACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGt                                                                                                  <  1:1638670/150‑1 (MQ=255)
                                                    aCCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGTCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGtg                                                                                                 <  1:1228589/150‑1 (MQ=255)
                                                    aCCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGtg                                                                                                 >  2:188027/1‑150 (MQ=255)
                                                      cacaTTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGtgat                                                                                               >  2:1662679/1‑150 (MQ=255)
                                                           tCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGTCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAg                                                                                          <  2:708735/150‑1 (MQ=255)
                                                            cTTTATCATCACTATTTTCATACTCCAGTCGATACCGTTACCGGGAACCTCTATCAGGTTGCATTTGTCGGCATTGGCGCGAATCAGGATGTTGCTGCTGGGGGTAGCGGTGCGAATATCGCGTTTCTTGCCTGTACGGGtgatgtctgc                                                                                         <  2:646013/150‑6 (MQ=255)
                                                            cTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGTCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGc                                                                                         <  1:3476015/150‑1 (MQ=255)
                                                               aaTCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCAtt                                                                                      >  1:4133833/1‑150 (MQ=255)
                                                               aaTCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGNAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTANGGGTGATGACAGCAtn                                                                                      <  2:2974224/150‑2 (MQ=255)
                                                                 tCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTgg                                                                                    <  1:2929067/150‑1 (MQ=255)
                                                                  cTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGt                                                                                   <  2:2682793/150‑1 (MQ=255)
                                                                     cAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTgag                                                                                <  2:3471421/150‑1 (MQ=255)
                                                                      aaTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGTCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTgaga                                                                               >  1:1605044/1‑150 (MQ=255)
                                                                      aaTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTgaga                                                                               >  1:2671624/1‑150 (MQ=255)
                                                                      aaTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTgaga                                                                               <  2:793797/150‑1 (MQ=255)
                                                                       aTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGagag                                                                              >  1:809000/1‑150 (MQ=255)
                                                                       aTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGagag                                                                              >  1:919021/1‑150 (MQ=255)
                                                                          ttttCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGt                                                                           >  1:3987235/1‑150 (MQ=255)
                                                                            ttCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTc                                                                         >  1:3710269/1‑150 (MQ=255)
                                                                            ttCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTc                                                                         <  1:2731448/150‑1 (MQ=255)
                                                                            ttCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTc                                                                         >  1:3564528/1‑150 (MQ=255)
                                                                               cTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGcc                                                                      >  1:1550950/1‑150 (MQ=255)
                                                                                tACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCca                                                                     <  2:193577/150‑1 (MQ=255)
                                                                                  cTCCACGCGAAACCGTAACCGGTAACCACTCTCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCcagc                                                                   <  1:3458098/150‑1 (MQ=255)
                                                                                    ccAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCcagcag                                                                 >  1:1917244/1‑150 (MQ=255)
                                                                                      aGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCt                                                               >  2:301912/1‑150 (MQ=255)
                                                                                      aGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCt                                                               >  2:590060/1‑150 (MQ=255)
                                                                                       ggCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTc                                                              <  1:3133918/150‑1 (MQ=255)
                                                                                          gAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGcc                                                           >  2:396191/1‑150 (MQ=255)
                                                                                          gAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGcc                                                           <  2:970609/150‑1 (MQ=255)
                                                                                           aaaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCa                                                          >  2:446197/1‑150 (MQ=255)
                                                                                           aaaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCa                                                          >  2:2100283/1‑150 (MQ=255)
                                                                                           aaaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCa                                                          <  1:332254/150‑1 (MQ=255)
                                                                                           aaaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCa                                                          >  2:548631/1‑150 (MQ=255)
                                                                                            gaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCAc                                                         >  2:4159692/2‑150 (MQ=255)
                                                                                            aaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGTCGCGACTCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATAACGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGAGTTCGCCAGCAGCTCGCCAc                                                         >  2:1539829/1‑150 (MQ=255)
                                                                                              ccGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTg                                                       <  1:3244342/150‑1 (MQ=255)
                                                                                               cGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGa                                                      <  2:1756250/150‑1 (MQ=255)
                                                                                                 tAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAATGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCGGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATa                                                    >  2:1561318/1‑150 (MQ=255)
                                                                                                 tAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATa                                                    >  2:2177583/1‑150 (MQ=255)
                                                                                                  aaCCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATaa                                                   >  2:161435/1‑150 (MQ=255)
                                                                                                  aaCCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATaa                                                   >  2:160356/1‑150 (MQ=255)
                                                                                                    ccGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACa                                                 <  1:650079/150‑1 (MQ=255)
                                                                                                     cGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                <  2:2548745/150‑1 (MQ=255)
                                                                                                        tAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGa                                             >  2:3752495/1‑150 (MQ=255)
                                                                                                        tAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGa                                             >  2:317035/1‑150 (MQ=255)
                                                                                                          aCCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGata                                           >  1:3803878/1‑150 (MQ=255)
                                                                                                           ccACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGatat                                          <  1:1320696/150‑1 (MQ=255)
                                                                                                              cTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTAc                                       >  1:2928966/1‑150 (MQ=255)
                                                                                                                aTCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTc                                     >  2:2313380/1‑150 (MQ=255)
                                                                                                                    ggTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGt                                 >  2:2908849/1‑150 (MQ=255)
                                                                                                                     gTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTc                                >  2:3669881/1‑150 (MQ=255)
                                                                                                                         cATTTGGCGGCAAAGGCGCGAATCAGTCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCg                            <  1:3837150/150‑1 (MQ=255)
                                                                                                                         cATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCg                            <  2:2251297/150‑1 (MQ=255)
                                                                                                                           tttGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTg                          <  2:3846221/150‑1 (MQ=255)
                                                                                                                            ttGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGa                         <  1:2522125/150‑1 (MQ=255)
                                                                                                                             tGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGAt                        <  1:3147203/150‑1 (MQ=255)
                                                                                                                              ggcggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATc                       >  1:1346615/1‑150 (MQ=255)
                                                                                                                                cggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCaa                     >  1:3410846/1‑150 (MQ=255)
                                                                                                                                 ggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCaaa                    >  1:1504852/1‑150 (MQ=255)
                                                                                                                                 ggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCaaa                    >  2:2441105/1‑150 (MQ=255)
                                                                                                                                       ggCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGaa              >  2:745217/1‑150 (MQ=255)
                                                                                                                                       ggCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGaa              >  2:3331660/1‑150 (MQ=255)
                                                                                                                                           cgAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCaa          >  2:3760997/1‑150 (MQ=255)
                                                                                                                                                cAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgt     <  2:270082/150‑1 (MQ=255)
                                                                                                                                                  ggCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgtgt   <  2:3589320/150‑1 (MQ=255)
                                                                                                                                                   gCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGtgtg  >  2:4157826/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTG  >  CP009273/3930648‑3930944

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: