Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,176,823 G→T 88.8% D516Y (GAC→TAC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,176,8230GT88.8% 863.4 / 100.4 305D516Y (GAC→TAC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base T (124/147);  total (141/164)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 305 total aligned reads displayed.

CTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTA  >  CP009273/4176677‑4176969
                                                                                                                                                  |                                                                                                                                                  
cTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACt                                                                                                                                                 <  1:2859510/150‑1 (MQ=255)
 tGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACtg                                                                                                                                                >  1:1021071/1‑150 (MQ=255)
 tGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACtg                                                                                                                                                <  2:2676591/150‑1 (MQ=255)
      aCGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTa                                                                                                                                           <  2:4529620/150‑1 (MQ=255)
        gTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAAc                                                                                                                                         >  1:1157981/1‑150 (MQ=255)
          aCCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGc                                                                                                                                       >  1:3943654/1‑150 (MQ=255)
             gCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCaa                                                                                                                                    <  2:2775385/150‑1 (MQ=255)
               tGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAg                                                                                                                                  <  1:2439994/150‑1 (MQ=255)
                    cTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGaa                                                                                                                             <  2:132435/150‑1 (MQ=255)
                    cTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGaa                                                                                                                             <  1:193657/150‑1 (MQ=255)
                          gCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATg                                                                                                                       <  1:2532174/150‑1 (MQ=255)
                           cccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATgg                                                                                                                      <  1:2987328/150‑1 (MQ=255)
                            ccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGt                                                                                                                     <  2:1746304/150‑1 (MQ=255)
                            ccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGt                                                                                                                     <  2:3977720/150‑1 (MQ=255)
                            ccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGt                                                                                                                     <  1:3649828/150‑1 (MQ=255)
                                cTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTg                                                                                                                 <  2:1647979/150‑1 (MQ=255)
                                  ggAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGAc                                                                                                               <  2:1659010/150‑1 (MQ=255)
                                  ggAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGAc                                                                                                               <  2:1855682/150‑1 (MQ=255)
                                   gAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACt                                                                                                              >  1:3389033/1‑150 (MQ=255)
                                     aGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTgg                                                                                                            >  1:999372/1‑150 (MQ=255)
                                      gcgcGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGc                                                                                                           <  2:2670815/150‑1 (MQ=255)
                                       cgcgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGcc                                                                                                          <  1:4460046/150‑1 (MQ=255)
                                         cgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCg                                                                                                        <  1:1480824/150‑1 (MQ=255)
                                          gTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGa                                                                                                       >  2:3411482/1‑150 (MQ=255)
                                             cgcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGGTAACGCCAAAGGCGAAGGCATGGTACTGACTGGCCCGAaag                                                                                                    >  2:1455010/1‑150 (MQ=255)
                                              gcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaaga                                                                                                   <  2:3200652/150‑1 (MQ=255)
                                                 gatgatGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGc                                                                                                >  1:4330811/1‑150 (MQ=255)
                                                  atgatgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCa                                                                                               <  2:2440139/150‑1 (MQ=255)
                                                     atgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGaa                                                                                            <  2:2020280/150‑1 (MQ=255)
                                                      tgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAAc                                                                                           >  2:3752855/1‑150 (MQ=255)
                                                         cTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTc                                                                                        >  1:1382844/1‑150 (MQ=255)
                                                          tACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCt                                                                                       <  1:1547795/150‑1 (MQ=255)
                                                           aCCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTg                                                                                      <  1:889783/150‑1 (MQ=255)
                                                                caacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCg                                                                                 <  2:873470/150‑1 (MQ=255)
                                                                caacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCg                                                                                 <  2:223837/150‑1 (MQ=255)
                                                                  acaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGct                                                                               <  1:3290613/150‑1 (MQ=255)
                                                                  acaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGct                                                                               <  1:985871/150‑1 (MQ=255)
                                                                  acaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGct                                                                               <  2:4535207/150‑1 (MQ=255)
                                                                  acaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGct                                                                               <  1:3579757/150‑1 (MQ=255)
                                                                   caTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGctc                                                                              >  1:1238924/1‑150 (MQ=255)
                                                                      ccTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctgg                                                                           >  1:2795743/1‑150 (MQ=255)
                                                                         gTCCCCGGCGAACGGCTAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGTTATGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtct                                                                        <  1:5748/150‑1 (MQ=255)
                                                                           ccccGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtctgg                                                                      >  2:1523949/1‑150 (MQ=255)
                                                                            cccGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGc                                                                     >  2:206682/1‑150 (MQ=255)
                                                                            cccGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGc                                                                     >  1:3437546/1‑150 (MQ=255)
                                                                              cggcgaacggcgaaCCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCtt                                                                   <  2:4383711/150‑1 (MQ=255)
                                                                                 cgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtctc                                                                <  2:3599478/150‑1 (MQ=255)
                                                                                 cgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtctc                                                                >  2:353645/1‑150 (MQ=255)
                                                                                 cgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtctc                                                                <  1:2247240/150‑1 (MQ=255)
                                                                                  gaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTctct                                                               <  1:1770471/150‑1 (MQ=255)
                                                                                  gaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTctct                                                               <  1:1659162/150‑1 (MQ=255)
                                                                                    acggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGc                                                             >  1:21190/1‑150 (MQ=255)
                                                                                     cggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCa                                                            <  2:1923042/150‑1 (MQ=255)
                                                                                      ggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCAt                                                           <  1:317703/150‑1 (MQ=255)
                                                                                      ggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCAt                                                           <  2:1414481/150‑1 (MQ=255)
                                                                                      ggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCAt                                                           <  2:2140911/150‑1 (MQ=255)
                                                                                         gaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcg                                                        <  2:490500/150‑1 (MQ=255)
                                                                                             cAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGcgcgcg                                                    >  1:461393/1‑150 (MQ=255)
                                                                                               atcatcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGtt                                                  <  1:635259/150‑1 (MQ=255)
                                                                                               atcatcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGtt                                                  <  1:4109185/150‑1 (MQ=255)
                                                                                                 catcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaa                                                >  2:2203652/1‑150 (MQ=255)
                                                                                                  atcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaaa                                               >  2:3709807/1‑150 (MQ=255)
                                                                                                  atcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaaa                                               >  1:1453083/1‑150 (MQ=255)
                                                                                                   tcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAg                                              <  1:1957565/150‑1 (MQ=255)
                                                                                                   tcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAg                                              <  1:854239/150‑1 (MQ=255)
                                                                                                     gTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTg                                            >  2:1203768/1‑150 (MQ=255)
                                                                                                         cGTCTCAGGACGTTGTACTGGGTCTGTACTACATGCCCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGACCGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTa                                        <  2:2833337/150‑1 (MQ=255)
                                                                                                          gTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTAt                                       <  2:3646268/150‑1 (MQ=255)
                                                                                                             tcAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCAc                                    >  1:3569813/1‑150 (MQ=255)
                                                                                                             tcAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCAc                                    <  1:360125/150‑1 (MQ=255)
                                                                                                              cAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCAcc                                   <  2:1656924/150‑1 (MQ=255)
                                                                                                                ggACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGa                                 >  2:750561/1‑150 (MQ=255)
                                                                                                                ggACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGa                                 >  2:4093113/1‑150 (MQ=255)
                                                                                                                 gACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAg                                >  2:2515022/1‑150 (MQ=255)
                                                                                                                    gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTCAAGTGCGTATCACCGAGTAt                             >  1:3123852/1‑150 (MQ=255)
                                                                                                                    gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTAt                             <  2:420768/150‑1 (MQ=255)
                                                                                                                    gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTAt                             <  2:4025408/150‑1 (MQ=255)
                                                                                                                    gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTAt                             >  1:1089474/1‑150 (MQ=255)
                                                                                                                      tGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGa                           >  2:2848114/1‑150 (MQ=255)
                                                                                                                       gTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGaa                          >  1:2958423/1‑150 (MQ=255)
                                                                                                                        tACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGaaa                         <  2:2101979/150‑1 (MQ=255)
                                                                                                                             ggTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGAt                    >  2:2919075/1‑150 (MQ=255)
                                                                                                                                cTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCt                 >  1:4535265/1‑150 (MQ=255)
                                                                                                                                 tGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTa                <  1:3911027/150‑1 (MQ=255)
                                                                                                                                   tactacATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAAc              >  1:1726347/1‑150 (MQ=255)
                                                                                                                                   tactacATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAAc              <  1:3689202/150‑1 (MQ=255)
                                                                                                                                    actacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAATAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACg             <  2:29878/150‑1 (MQ=255)
                                                                                                                                      tacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGt           >  2:3865179/1‑150 (MQ=255)
                                                                                                                                         aTGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGaa        >  1:2594391/1‑150 (MQ=255)
                                                                                                                                         aTGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGaa        <  2:4126249/150‑1 (MQ=255)
                                                                                                                                          tGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAAt       <  2:2301035/150‑1 (MQ=255)
                                                                                                                                           gACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACTGTGAAtt      >  2:4482439/1‑150 (MQ=255)
                                                                                                                                           gACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAAtt      >  2:923835/1‑150 (MQ=255)
                                                                                                                                           gACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAAtt      >  1:3202548/1‑150 (MQ=255)
                                                                                                                                            aCCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATta     >  1:4036015/1‑150 (MQ=255)
                                                                                                                                            aCCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATta     <  2:316195/150‑1 (MQ=255)
                                                                                                                                             cccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtag    <  1:3836535/150‑1 (MQ=255)
                                                                                                                                              ccGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagt   <  1:2080123/150‑1 (MQ=255)
                                                                                                                                               cGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagta  >  1:3776745/1‑150 (MQ=255)
                                                                                                                                               cGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagta  >  1:1865167/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTA  >  CP009273/4176677‑4176969

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: