Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,875,627 C→A 5.4% G19C (GGT→TGT)  dnaN ← DNA polymerase III, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,875,6270CA 5.4% 1408.2 / 32.9 350G19C (GGT→TGT) dnaNDNA polymerase III, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (146/185);  new base A (6/13);  total (152/198)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 364 total aligned reads displayed.

CGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGA  >  CP009273/3875478‑3875776
                                                                                                                                                     |                                                                                                                                                     
cgtcgCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGAcc                                                                                                                                                       >  1:4125803/1‑150 (MQ=255)
cgtcgCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGAcc                                                                                                                                                       <  1:3143122/150‑1 (MQ=255)
   cgCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCt                                                                                                                                                    >  2:3427516/1‑150 (MQ=255)
       ccTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCAcc                                                                                                                                                >  1:2429691/1‑150 (MQ=255)
        cTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCt                                                                                                                                               >  2:3777292/1‑150 (MQ=255)
         tGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCtg                                                                                                                                              >  1:4921103/1‑150 (MQ=255)
         tGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACAGCTCACCtg                                                                                                                                              >  2:893105/1‑150 (MQ=255)
            cTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCtgttg                                                                                                                                           >  2:1151599/1‑150 (MQ=255)
              cGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTa                                                                                                                                         >  1:2493426/1‑150 (MQ=255)
                tgtgGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGc                                                                                                                                       >  2:3336486/1‑150 (MQ=255)
                     cTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGttt                                                                                                                                  <  2:1003381/150‑1 (MQ=255)
                      tGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGtttt                                                                                                                                 <  2:3154002/150‑1 (MQ=255)
                       gAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTta                                                                                                                                <  2:71464/150‑1 (MQ=255)
                         aCCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGTACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTGTtaag                                                                                                                              >  1:3360926/1‑149 (MQ=255)
                         aCCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTtaat                                                                                                                              >  2:4359370/1‑150 (MQ=255)
                          ccAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTtaata                                                                                                                             <  2:4084501/150‑1 (MQ=255)
                          ccAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTtaata                                                                                                                             >  1:4277747/1‑150 (MQ=255)
                            aGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATaaa                                                                                                                           >  1:2136257/1‑150 (MQ=255)
                             gcgcAACACGTGCCACCATTTCCATCGCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAAt                                                                                                                          <  1:2960113/150‑1 (MQ=255)
                               gcAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGc                                                                                                                        <  1:3592589/150‑1 (MQ=255)
                               gcAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGc                                                                                                                        >  1:4824498/1‑150 (MQ=255)
                                cAACACGTGCCACCATTTCCCTCTCGAGATCAGTCCCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCt                                                                                                                       <  2:3216853/150‑1 (MQ=255)
                                 aaCACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTc                                                                                                                      <  1:766189/150‑1 (MQ=255)
                                  acacGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACAGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCa                                                                                                                     >  1:4313086/1‑150 (MQ=255)
                                      gTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGtt                                                                                                                 <  2:3710118/150‑1 (MQ=255)
                                        gCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCt                                                                                                               >  1:1407684/1‑150 (MQ=255)
                                         ccaccaTTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTa                                                                                                              >  1:2050362/1‑150 (MQ=255)
                                         ccaccaTTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTa                                                                                                              <  2:173125/150‑1 (MQ=255)
                                          caccaTTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTAc                                                                                                             <  1:304241/150‑1 (MQ=255)
                                           accaTTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACg                                                                                                            <  2:2856308/150‑1 (MQ=255)
                                              aTTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTa                                                                                                         <  1:5134/150‑1 (MQ=255)
                                                ttCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGGAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTaaa                                                                                                       <  2:3077718/150‑1 (MQ=255)
                                                 tCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAAt                                                                                                      <  1:2655566/150‑1 (MQ=255)
                                                  ccATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAAtt                                                                                                     <  1:2447956/150‑1 (MQ=255)
                                                    aTCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTTGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTc                                                                                                   <  1:4466740/150‑1 (MQ=255)
                                                       tcGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATa                                                                                                >  2:1421054/1‑150 (MQ=255)
                                                         gagaTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAgg                                                                                              >  2:44533/1‑150 (MQ=255)
                                                            aTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCAGAGGttt                                                                                           >  2:2686474/1‑150 (MQ=255)
                                                              cAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTAc                                                                                         <  2:857333/150‑1 (MQ=255)
                                                                 tACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGAt                                                                                      <  2:4632262/150‑1 (MQ=255)
                                                                 gaCCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGAt                                                                                      <  2:2331363/149‑1 (MQ=255)
                                                                  aCCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATg                                                                                     >  2:2175417/1‑150 (MQ=255)
                                                                     ggTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACa                                                                                  >  2:368089/1‑150 (MQ=255)
                                                                      gTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACaa                                                                                 <  2:2244897/150‑1 (MQ=255)
                                                                          gCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGtt                                                                             <  1:1003086/150‑1 (MQ=255)
                                                                          gCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGtt                                                                             <  1:473701/150‑1 (MQ=255)
                                                                           cGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTc                                                                            <  1:1154383/150‑1 (MQ=255)
                                                                            gACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCt                                                                           <  1:5058581/150‑1 (MQ=255)
                                                                             aCAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTg                                                                          <  1:2855024/150‑1 (MQ=255)
                                                                             aCAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTg                                                                          >  1:1957721/1‑150 (MQ=255)
                                                                             aCAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTg                                                                          <  2:2042255/150‑1 (MQ=255)
                                                                              cAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGa                                                                         <  2:1222452/150‑1 (MQ=255)
                                                                               aaCGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGAt                                                                        <  1:3525356/150‑1 (MQ=255)
                                                                                aCGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGAtt                                                                       <  1:956042/150‑1 (MQ=255)
                                                                                aCGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGAtt                                                                       <  1:4788029/150‑1 (MQ=255)
                                                                                 cGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTa                                                                      >  1:2218213/1‑150 (MQ=255)
                                                                                 cGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTa                                                                      <  2:4204860/150‑1 (MQ=255)
                                                                                    aCCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAAt                                                                   <  2:3527137/150‑1 (MQ=255)
                                                                                      cGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATCTTCTGATTAAAttt                                                                 >  2:4223571/1‑150 (MQ=255)
                                                                                           gCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGaaaa                                                            <  1:3312377/150‑1 (MQ=255)
                                                                                           gCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGaaaa                                                            >  1:39307/1‑150 (MQ=255)
                                                                                            cAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGaaaaa                                                           >  2:1423476/1‑150 (MQ=255)
                                                                                             aaCCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAAt                                                          <  2:2955719/150‑1 (MQ=255)
                                                                                               ccTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAAtct                                                        <  1:4211265/150‑1 (MQ=255)
                                                                                               ccTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACAGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAAtct                                                        <  1:5108694/150‑1 (MQ=255)
                                                                                                cTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAAtctt                                                       >  2:4861157/1‑150 (MQ=255)
                                                                                                cTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAAtctt                                                       >  1:3379755/1‑150 (MQ=255)
                                                                                                   tAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATcttctt                                                    <  1:2982506/150‑1 (MQ=255)
                                                                                                   tAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATcttctt                                                    >  1:3242681/1‑150 (MQ=255)
                                                                                                   tAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATcttctt                                                    <  2:4659446/150‑1 (MQ=255)
                                                                                                      cagcagATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGa                                                 >  2:1199200/1‑150 (MQ=255)
                                                                                                        gcagATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGata                                               <  2:1267599/150‑1 (MQ=255)
                                                                                                        gcagATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGata                                               <  2:2702919/150‑1 (MQ=255)
                                                                                                           gATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCg                                            <  1:3025741/150‑1 (MQ=255)
                                                                                                           gATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACAGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCg                                            <  2:3646020/150‑1 (MQ=255)
                                                                                                                    gaATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGctctct                                   <  1:4605651/150‑1 (MQ=255)
                                                                                                                     aaTCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCtt                                  >  2:4006770/1‑150 (MQ=255)
                                                                                                                        cGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCAc                               <  2:1966185/150‑1 (MQ=255)
                                                                                                                        cGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCGGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCAc                               >  1:4094709/1‑150 (MQ=255)
                                                                                                                         ggTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTATTTGATATCGTGGCTCTCTTCACg                              >  2:4134250/1‑150 (MQ=255)
                                                                                                                              cGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACt                         >  1:1441051/1‑150 (MQ=255)
                                                                                                                               gTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTg                        >  2:1634156/1‑150 (MQ=255)
                                                                                                                               gTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTg                        >  1:1836042/1‑150 (MQ=255)
                                                                                                                               gTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTg                        >  1:2507945/1‑150 (MQ=255)
                                                                                                                               gTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTg                        <  2:4580799/150‑1 (MQ=255)
                                                                                                                               gTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTg                        >  2:4752299/1‑150 (MQ=255)
                                                                                                                                     cgacCACCTAACGGACAGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGAt                  >  1:1684176/1‑150 (MQ=255)
                                                                                                                                       accaccTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCt                >  1:2687275/1‑150 (MQ=255)
                                                                                                                                          accTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTAc             <  2:809249/150‑1 (MQ=255)
                                                                                                                                          accTAACGGACAGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTAc             <  1:111243/150‑1 (MQ=255)
                                                                                                                                             tAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggc          >  2:605391/1‑150 (MQ=255)
                                                                                                                                              aaCGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggca         >  1:3832644/1‑150 (MQ=255)
                                                                                                                                              aaCGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggca         >  1:4049229/1‑150 (MQ=255)
                                                                                                                                                cGGACAGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggcagg       <  1:2895617/150‑1 (MQ=255)
                                                                                                                                                 ggACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggcaggc      >  2:2545983/1‑150 (MQ=255)
                                                                                                                                                    ccGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATg   >  2:1351137/1‑150 (MQ=255)
                                                                                                                                                    ccGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATg   >  1:1179947/1‑150 (MQ=255)
                                                                                                                                                     cGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGa  >  1:4275384/1‑150 (MQ=255)
                                                                                                                                                     cGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGa  >  2:1117142/1‑150 (MQ=255)
                                                                                                                                                     cGCTCAACTGTTGTAGCTGTTTTAATAAATGCTCACGTTCTACTGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCAtta  <  2:1827340/150‑3 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTCGCTCCTGGCTCGTGTGGCTGAACCAGCGCAACACGTGCCACCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGA  >  CP009273/3875478‑3875776

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: