Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,217,289 A→T 100% H433L (CAC→CTC)  ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,217,2890AT98.0% 647.2 / 2.1 202H433L (CAC→CTC) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base T (101/97);  total (101/101)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 202 total aligned reads displayed.

GAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGA  >  CP009273/3217144‑3217433
                                                                                                                                                 |                                                                                                                                                
gAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTcg                                                                                                                                                >  2:119470/1‑148 (MQ=255)
gAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTcg                                                                                                                                                <  1:244845/148‑1 (MQ=255)
gAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTcg                                                                                                                                                >  2:3082855/1‑148 (MQ=255)
   cTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCgcgc                                                                                                                                             >  2:3050668/1‑148 (MQ=255)
     gtgtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAg                                                                                                                                           >  2:2516477/1‑148 (MQ=255)
           tatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAAc                                                                                                                                     >  1:2633060/1‑148 (MQ=255)
           tatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAAc                                                                                                                                     >  1:1985553/1‑148 (MQ=255)
            atCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAAcc                                                                                                                                    <  2:462032/148‑1 (MQ=255)
             tCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCa                                                                                                                                   <  2:1781243/148‑1 (MQ=255)
              cTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCAt                                                                                                                                  <  1:1890177/148‑1 (MQ=255)
               tACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATc                                                                                                                                 <  1:2105178/148‑1 (MQ=255)
                   gCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTc                                                                                                                             >  2:2561679/1‑148 (MQ=255)
                    ccTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCg                                                                                                                            >  1:3421299/1‑148 (MQ=255)
                    ccTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCg                                                                                                                            <  2:572635/148‑1 (MQ=255)
                     cTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGa                                                                                                                           >  2:3410273/1‑148 (MQ=255)
                      tGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGat                                                                                                                          <  2:1236659/148‑1 (MQ=255)
                       gTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatc                                                                                                                         <  2:223809/148‑1 (MQ=255)
                         ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcat                                                                                                                       <  2:2325874/148‑1 (MQ=255)
                            tgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcat                                                                                                                    <  2:3451628/148‑1 (MQ=255)
                             gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcatc                                                                                                                   <  2:2967435/148‑1 (MQ=255)
                             gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcatc                                                                                                                   >  2:1910595/1‑148 (MQ=255)
                             gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcatc                                                                                                                   >  2:1710273/1‑148 (MQ=255)
                                ggCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTgg                                                                                                                <  2:1996341/148‑1 (MQ=255)
                                 gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGt                                                                                                               <  2:3557296/148‑1 (MQ=255)
                                 gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGt                                                                                                               <  2:1512402/148‑1 (MQ=255)
                                 gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGt                                                                                                               <  2:2477345/148‑1 (MQ=255)
                                     aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                           <  2:2147849/148‑1 (MQ=255)
                                      aaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTg                                                                                                          >  2:1923493/1‑148 (MQ=255)
                                      aaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTg                                                                                                          <  1:83841/148‑1 (MQ=255)
                                        ccGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTggg                                                                                                        >  1:3202077/1‑148 (MQ=255)
                                         cGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGc                                                                                                       >  2:1616954/1‑148 (MQ=255)
                                          gACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCa                                                                                                      <  2:605308/148‑1 (MQ=255)
                                           aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCaa                                                                                                     >  1:1791509/1‑148 (MQ=255)
                                           aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCaa                                                                                                     <  1:2368363/148‑1 (MQ=255)
                                            cGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAAt                                                                                                    <  2:703438/148‑1 (MQ=255)
                                            cGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAAt                                                                                                    <  1:3261293/148‑1 (MQ=255)
                                              tCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGa                                                                                                  >  1:2748045/1‑148 (MQ=255)
                                               cGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGaa                                                                                                 <  1:2347186/148‑1 (MQ=255)
                                               cGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGaa                                                                                                 >  2:3493873/1‑148 (MQ=255)
                                                 aaTCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATc                                                                                               <  1:60351/148‑1 (MQ=255)
                                                 aaTCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATc                                                                                               >  1:1117237/1‑148 (MQ=255)
                                                   tCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCg                                                                                             >  2:2913945/1‑148 (MQ=255)
                                                     gCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGc                                                                                           <  1:2932893/148‑1 (MQ=255)
                                                        cGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGTAAAACCATCCGTCGATCATCCTCTGGTCGCTGGGCAATGAATCCGGCTAt                                                                                        >  1:3317384/1‑148 (MQ=255)
                                                        cGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTAt                                                                                        >  2:2617689/1‑148 (MQ=255)
                                                        cGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTAt                                                                                        >  1:1549999/1‑148 (MQ=255)
                                                          gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATgg                                                                                      <  1:3529966/148‑1 (MQ=255)
                                                              tGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGt                                                                                  <  1:2758374/148‑1 (MQ=255)
                                                               gCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATAATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTa                                                                                 <  2:3516619/148‑1 (MQ=255)
                                                                 tAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAAc                                                                               <  1:1220003/148‑1 (MQ=255)
                                                                 tAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAAc                                                                               <  2:3239308/148‑1 (MQ=255)
                                                                 tAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAAc                                                                               <  1:1085562/148‑1 (MQ=255)
                                                                  aaTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACa                                                                              <  2:2456659/148‑1 (MQ=255)
                                                                   aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  2:1239927/148‑1 (MQ=255)
                                                                      tCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcg                                                                          <  1:257122/148‑1 (MQ=255)
                                                                          cGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcgcgcg                                                                      <  2:901426/148‑1 (MQ=255)
                                                                            atatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGAt                                                                    >  1:173279/1‑148 (MQ=255)
                                                                                 aGCCGTATTACCGACTATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCTGCTATGGCTGTAACATCCGCGCGATGTAcc                                                               <  1:902208/148‑1 (MQ=255)
                                                                                 aGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTAcc                                                               <  1:1005722/148‑1 (MQ=255)
                                                                                  gCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCa                                                              >  2:1889447/1‑148 (MQ=255)
                                                                                     gTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgc                                                           >  1:688110/1‑148 (MQ=255)
                                                                                      tATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgcg                                                          >  2:2636988/1‑148 (MQ=255)
                                                                                       aTTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgcgg                                                         >  1:2605738/1‑148 (MQ=255)
                                                                                        ttACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgcggc                                                        >  1:644718/1‑148 (MQ=255)
                                                                                           ccGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaa                                                     <  1:2554435/148‑1 (MQ=255)
                                                                                             gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  1:3050668/148‑1 (MQ=255)
                                                                                               cgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCAGATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGTCTATGGATGTAACATCCGCGCGATGTACCATGCGGCGAAAgcg                                                 <  1:622643/148‑1 (MQ=255)
                                                                                                gaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAgcgc                                                <  2:1550406/148‑1 (MQ=255)
                                                                                                gaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAACGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAgcgc                                                <  1:2315744/148‑1 (MQ=255)
                                                                                                 aTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCt                                               <  2:1027079/148‑1 (MQ=255)
                                                                                                  tCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTg                                              >  1:722163/1‑148 (MQ=255)
                                                                                                   ccGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTgg                                             <  1:2137700/148‑1 (MQ=255)
                                                                                                      cAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATg                                          >  2:504841/1‑148 (MQ=255)
                                                                                                      cAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATg                                          >  2:1486531/1‑148 (MQ=255)
                                                                                                       aGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGa                                         <  1:2192934/148‑1 (MQ=255)
                                                                                                         tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                       >  1:772399/1‑148 (MQ=255)
                                                                                                           ggAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcg                                     >  2:2483502/1‑148 (MQ=255)
                                                                                                             aaaaaGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcgcg                                   <  2:1117237/148‑1 (MQ=255)
                                                                                                             aaaaaGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcgcg                                   >  2:2370558/1‑148 (MQ=255)
                                                                                                                   tCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTCTGGCTGTAACATCCGCGCGAGGTACCATGCGGCGTAAGCGCTGGATGACACGCGACTGGt                             >  2:2389358/1‑148 (MQ=255)
                                                                                                                    cTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTg                            >  1:558837/1‑148 (MQ=255)
                                                                                                                      aCGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCa                          <  1:1181155/148‑1 (MQ=255)
                                                                                                                      aCGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCa                          <  2:2266239/148‑1 (MQ=255)
                                                                                                                                  ttGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAagaaga              <  1:2516477/148‑1 (MQ=255)
                                                                                                                                  ttGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAagaaga              <  1:1616954/148‑1 (MQ=255)
                                                                                                                                     ttCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcg           >  2:519800/1‑148 (MQ=255)
                                                                                                                                     ttCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcg           >  1:3491390/1‑148 (MQ=255)
                                                                                                                                     ttCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcg           <  1:481558/148‑1 (MQ=255)
                                                                                                                                     ttCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcg           >  1:1707266/1‑148 (MQ=255)
                                                                                                                                       cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg         >  1:945153/1‑148 (MQ=255)
                                                                                                                                        gCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGa        >  1:2024980/1‑148 (MQ=255)
                                                                                                                                        gCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGa        <  2:157145/148‑1 (MQ=255)
                                                                                                                                        gCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGa        <  2:167687/148‑1 (MQ=255)
                                                                                                                                         ccATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGAt       >  1:3002148/1‑148 (MQ=255)
                                                                                                                                          cATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGATCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATg      <  1:1710273/148‑1 (MQ=255)
                                                                                                                                           atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     <  1:3119170/148‑1 (MQ=255)
                                                                                                                                           atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     <  1:15008/148‑1 (MQ=255)
                                                                                                                                           atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     <  1:1277303/148‑1 (MQ=255)
                                                                                                                                           atatCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     <  2:1927378/148‑1 (MQ=255)
                                                                                                                                            tatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCt    >  1:2869895/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGA  >  CP009273/3217144‑3217433

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