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,2890AT99.5% 624.7 / ‑4.6 191H433L (CAC→CTC) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (84/106);  total (85/106)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 191 total aligned reads displayed.

AACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGA  >  CP009273/3217145‑3217433
                                                                                                                                                |                                                                                                                                                
aaCTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTcgc                                                                                                                                               >  1:3203829/1‑148 (MQ=255)
       tgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaa                                                                                                                                        <  1:451744/148‑1 (MQ=255)
        gATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaaa                                                                                                                                       <  1:222259/148‑1 (MQ=255)
         atatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaaaa                                                                                                                                      >  2:1026211/1‑148 (MQ=255)
         atatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaaaa                                                                                                                                      >  2:2308010/1‑148 (MQ=255)
            tCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCa                                                                                                                                   >  1:1117362/1‑148 (MQ=255)
            tCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCa                                                                                                                                   >  1:909368/1‑148 (MQ=255)
             cTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCAt                                                                                                                                  <  2:2690099/148‑1 (MQ=255)
              tACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATc                                                                                                                                 <  2:1334469/148‑1 (MQ=255)
                   ccTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCg                                                                                                                            <  1:21169/148‑1 (MQ=255)
                       tttGTGATGGAGGAAACAGACGACGAATAGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCATTGGGAAAAATTCTACGTAGAGCGCAATGTACGCCATATCCTCGCGCAGAAAAACCATCAGACGaaca                                                                                                                        >  2:1843641/1‑148 (MQ=255)
                        ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcat                                                                                                                       >  1:65887/1‑148 (MQ=255)
                         tgtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatc                                                                                                                      <  1:2454422/148‑1 (MQ=255)
                           tgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcat                                                                                                                    <  2:225708/148‑1 (MQ=255)
                            gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcatc                                                                                                                   <  1:2120976/148‑1 (MQ=255)
                               ggCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTgg                                                                                                                >  2:2058296/1‑148 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                           >  2:1248156/1‑148 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                           <  2:2661742/148‑1 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                           >  2:1983944/1‑148 (MQ=255)
                                     aaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTg                                                                                                          <  1:402548/148‑1 (MQ=255)
                                      aCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTgg                                                                                                         >  2:3003972/1‑148 (MQ=255)
                                           cGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAATTCTACTTAGAGCGCATTGTTCGACATATACTAGCGCAGAAAAACCATCCGACGATCATCATCTGGTCGCTGGGCAAt                                                                                                    >  2:3302935/1‑148 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAAt                                                                                                <  2:1927291/148‑1 (MQ=255)
                                                      aCGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATCTCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTa                                                                                         <  2:535281/148‑1 (MQ=255)
                                                        ggCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATg                                                                                       >  1:3341672/1‑148 (MQ=255)
                                                         gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATgg                                                                                      <  2:825993/148‑1 (MQ=255)
                                                         gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATgg                                                                                      >  1:96977/1‑148 (MQ=255)
                                                           tttGCTTATGTCGGCGATTTTAGCCGTATTACCGACGATCCTCAGTGGGAAAAAGTCTACGTCGAGCGCTTTGTTCGCCATATCCTCGCGCAGTAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCt                                                                                    <  1:1470103/148‑1 (MQ=255)
                                                             tGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGt                                                                                  >  1:2268046/1‑148 (MQ=255)
                                                                tAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAAc                                                                               <  1:1659513/148‑1 (MQ=255)
                                                                 aaTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACa                                                                              >  2:3217025/1‑148 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  2:3372133/148‑1 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  1:2925182/148‑1 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  1:2899011/148‑1 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             >  2:207309/1‑148 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             >  2:3331567/1‑148 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  1:1147286/148‑1 (MQ=255)
                                                                   tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCCGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATc                                                                            >  1:2400915/1‑148 (MQ=255)
                                                                   tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATc                                                                            >  2:3333882/1‑148 (MQ=255)
                                                                   tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATc                                                                            >  1:461555/1‑148 (MQ=255)
                                                                       ggCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcgcg                                                                        <  1:2722654/148‑1 (MQ=255)
                                                                       ggCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcgcg                                                                        >  2:3257714/1‑148 (MQ=255)
                                                                          gATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGa                                                                     <  1:1002445/148‑1 (MQ=255)
                                                                           atatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGAt                                                                    >  1:3499278/1‑148 (MQ=255)
                                                                            tatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATg                                                                   <  2:1141078/148‑1 (MQ=255)
                                                                            tatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATg                                                                   >  1:1614805/1‑148 (MQ=255)
                                                                             atTAGCCGTATTCCCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGt                                                                  <  2:2016578/148‑1 (MQ=255)
                                                                             atTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGt                                                                  >  1:1838623/1‑148 (MQ=255)
                                                                                aGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTAcc                                                               <  1:3065082/148‑1 (MQ=255)
                                                                                  ccGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCAt                                                             >  2:1750476/1‑148 (MQ=255)
                                                                                    gTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgc                                                           >  2:1606789/1‑148 (MQ=255)
                                                                                          ccGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaa                                                     <  2:3137470/148‑1 (MQ=255)
                                                                                          ccGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaa                                                     <  2:1503295/148‑1 (MQ=255)
                                                                                           cgacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaaa                                                    <  1:2857201/148‑1 (MQ=255)
                                                                                           cgacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaaa                                                    <  1:2532958/148‑1 (MQ=255)
                                                                                           cgacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaaa                                                    >  2:2523646/1‑148 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  2:2911099/148‑1 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  2:1161441/148‑1 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  2:1787491/148‑1 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  1:936313/148‑1 (MQ=255)
                                                                                              cgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAgcg                                                 <  1:1361664/148‑1 (MQ=255)
                                                                                                aTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCt                                               >  2:1837806/1‑148 (MQ=255)
                                                                                                  ccGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTgg                                             >  1:1499488/1‑148 (MQ=255)
                                                                                                  ccGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTgg                                             <  2:2464675/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  2:1747945/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  2:753592/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  2:1041028/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  1:1983944/148‑1 (MQ=255)
                                                                                                      aGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGa                                         <  2:2010169/148‑1 (MQ=255)
                                                                                                       gTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  1:626536/137‑1 (MQ=255)
                                                                                                       gTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   >  2:626736/1‑137 (MQ=255)
                                                                                                        tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                       >  1:2241452/1‑148 (MQ=255)
                                                                                                        tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                       >  2:2777340/1‑148 (MQ=255)
                                                                                                        tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                       <  2:2653080/148‑1 (MQ=255)
                                                                                                         gggAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGacac                                      >  1:1863368/1‑148 (MQ=255)
                                                                                                          ggAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcg                                     <  2:1840434/148‑1 (MQ=255)
                                                                                                           gAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcgc                                    >  1:2152149/1‑148 (MQ=255)
                                                                                                               aaGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACt                                >  1:2632253/1‑148 (MQ=255)
                                                                                                                aGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTg                               <  1:1671287/148‑1 (MQ=255)
                                                                                                                 gTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTgg                              <  2:3499278/148‑1 (MQ=255)
                                                                                                                      cGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCAt                         <  1:359868/148‑1 (MQ=255)
                                                                                                                        tCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGGGGCGAAAGCGCTGGTTGACACGCGACTGGTGCATTa                       <  1:1811304/148‑1 (MQ=255)
                                                                                                                        tCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTCGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTa                       >  2:828867/1‑148 (MQ=255)
                                                                                                                            gcgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaa                   >  2:463447/1‑148 (MQ=255)
                                                                                                                            gcgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaa                   <  1:3305855/148‑1 (MQ=255)
                                                                                                                             cgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaag                  <  2:912458/148‑1 (MQ=255)
                                                                                                                             cgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaag                  <  1:2201260/148‑1 (MQ=255)
                                                                                                                                aTTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGaagaag               <  2:362625/148‑1 (MQ=255)
                                                                                                                                 ttGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAagaaga              <  2:350053/148‑1 (MQ=255)
                                                                                                                                     tCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgc          >  2:2470667/1‑148 (MQ=255)
                                                                                                                                      cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg         <  2:928856/148‑1 (MQ=255)
                                                                                                                                      cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg         <  2:3198349/148‑1 (MQ=255)
                                                                                                                                      cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg         <  1:3390641/148‑1 (MQ=255)
                                                                                                                                        ccATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGAt       >  1:2617742/1‑148 (MQ=255)
                                                                                                                                         cATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATg      >  2:1656499/1‑148 (MQ=255)
                                                                                                                                         cATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATg      >  1:3237347/1‑148 (MQ=255)
                                                                                                                                         cATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATg      <  1:1026211/148‑1 (MQ=255)
                                                                                                                                          atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGACAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     >  2:2867240/1‑148 (MQ=255)
                                                                                                                                          atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     >  2:824129/1‑148 (MQ=255)
                                                                                                                                             tCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGa  <  2:2752279/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGA  >  CP009273/3217145‑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: