Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,278,846 C→A 5.8% I385I (ATC→ATA nrfA → nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,278,8460CA 5.8% 1193.1 / 30.4 308I385I (ATC→ATAnrfAnitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base C (151/139);  new base A (11/7);  total (162/146)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 322 total aligned reads displayed.

ATTAACGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTC  >  CP009273/4278697‑4278995
                                                                                                                                                     |                                                                                                                                                     
aTTAACGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATc                                                                                                                                                       <  1:2329728/150‑1 (MQ=255)
aTTAACGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATc                                                                                                                                                       >  1:1761863/1‑150 (MQ=255)
aTTAACGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATc                                                                                                                                                       >  1:629909/1‑150 (MQ=255)
aTTAACGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATc                                                                                                                                                       <  2:1521031/150‑1 (MQ=255)
  tAACGACCTGAAAATCAAGGTTGAAGATAAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCGGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATAACCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGc                                                                                                                                                     >  2:3216330/1‑150 (MQ=255)
    aCGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCtt                                                                                                                                                   <  1:4081952/150‑1 (MQ=255)
    aCGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCtt                                                                                                                                                   <  2:3575426/150‑1 (MQ=255)
      gACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTcc                                                                                                                                                 >  2:3945369/1‑150 (MQ=255)
       aCCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTccc                                                                                                                                                <  2:4263159/150‑1 (MQ=255)
        ccTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCa                                                                                                                                               >  2:3481466/1‑150 (MQ=255)
                tCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCTCTTCCCACGGCATTc                                                                                                                                       >  1:483431/1‑150 (MQ=255)
                tCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTc                                                                                                                                       >  1:1915050/1‑150 (MQ=255)
                  aaGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCat                                                                                                                                     >  2:547815/1‑150 (MQ=255)
                  aaGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCat                                                                                                                                     >  1:2310014/1‑150 (MQ=255)
                       ttaaGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATAGCTTCCCACGGCATTCATATgca                                                                                                                                >  2:2888890/3‑150 (MQ=255)
                        gAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATgcac                                                                                                                               >  1:2797537/1‑150 (MQ=255)
                          aGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATgcacgc                                                                                                                             >  2:826799/1‑150 (MQ=255)
                             tCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCAcc                                                                                                                          >  1:1302475/1‑150 (MQ=255)
                               aaCTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCgg                                                                                                                        >  2:2716808/1‑150 (MQ=255)
                                aCTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGga                                                                                                                       <  2:3613054/150‑1 (MQ=255)
                                 cTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGgaa                                                                                                                      <  2:813584/150‑1 (MQ=255)
                                     ttCACGCTCACATCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGaagaag                                                                                                                  >  2:3315875/1‑150 (MQ=255)
                                      tCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAgg                                                                                                                 <  2:3332806/150‑1 (MQ=255)
                                       cACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGt                                                                                                                <  1:1397169/150‑1 (MQ=255)
                                         cGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGttt                                                                                                              <  1:4177050/150‑1 (MQ=255)
                                          gCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAATCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTa                                                                                                             <  1:983913/150‑1 (MQ=255)
                                           cTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTAc                                                                                                            >  1:202250/1‑150 (MQ=255)
                                           cTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTAc                                                                                                            <  2:3969609/150‑1 (MQ=255)
                                            tCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACg                                                                                                           >  1:4006880/1‑150 (MQ=255)
                                             cACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACgg                                                                                                          <  2:1206222/150‑1 (MQ=255)
                                               cTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGAt                                                                                                        <  1:86477/150‑1 (MQ=255)
                                                ttCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATg                                                                                                       <  2:2057775/150‑1 (MQ=255)
                                                ttCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATg                                                                                                       <  2:3858745/150‑1 (MQ=255)
                                                 tCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGc                                                                                                      >  2:2412003/1‑150 (MQ=255)
                                                 tCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGc                                                                                                      >  2:1406052/1‑150 (MQ=255)
                                                  cGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCt                                                                                                     <  1:2945027/150‑1 (MQ=255)
                                                  cGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCt                                                                                                     <  2:3810883/150‑1 (MQ=255)
                                                   gAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTc                                                                                                    <  2:2099462/150‑1 (MQ=255)
                                                   gAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTc                                                                                                    <  2:3426979/150‑1 (MQ=255)
                                                   gAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATCCACGCACCGGAAGAAGGTTTACGGATGCTc                                                                                                    <  1:4278551/150‑1 (MQ=255)
                                                   gAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATAGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTc                                                                                                    <  1:4401709/150‑1 (MQ=255)
                                                     aGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCTCGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCgg                                                                                                  <  1:120980/150‑1 (MQ=255)
                                                     aGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCgg                                                                                                  >  1:2575365/1‑150 (MQ=255)
                                                      gCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGACGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGTTCGCTTGCCACGGCATTCATATGCACGCACCGGTAGAAGGTTTACGGATGCTCGGt                                                                                                 <  2:2363206/150‑1 (MQ=255)
                                                       cGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTa                                                                                                <  2:2680115/150‑1 (MQ=255)
                                                         aaaGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACg                                                                                              <  1:4604260/150‑1 (MQ=255)
                                                          aaGCAGCGCTGGATGCAGGCGCGACGGAAGGTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACgg                                                                                             >  1:3895319/1‑150 (MQ=255)
                                                           agcagcGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGc                                                                                            >  1:2975315/1‑150 (MQ=255)
                                                                 gcTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGa                                                                                      >  2:695180/1‑150 (MQ=255)
                                                                      aTGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAg                                                                                 <  1:2663862/150‑1 (MQ=255)
                                                                       tGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAgc                                                                                <  2:11025/150‑1 (MQ=255)
                                                                         cAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAgcgg                                                                              >  2:1061492/1‑150 (MQ=255)
                                                                         cAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAgcgg                                                                              <  1:1314954/150‑1 (MQ=255)
                                                                              gcgACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGAt                                                                         <  2:1862614/150‑1 (MQ=255)
                                                                                gACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATgc                                                                       <  1:3769034/150‑1 (MQ=255)
                                                                                  cGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATgcac                                                                     <  1:1432930/150‑1 (MQ=255)
                                                                                    gAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATgcacgc                                                                   >  1:1533174/1‑150 (MQ=255)
                                                                                    gAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATgcacgc                                                                   >  1:3728713/1‑150 (MQ=255)
                                                                                     aaGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATgcacgca                                                                  >  1:2466495/1‑150 (MQ=255)
                                                                                         tGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCaa                                                              >  2:3040296/1‑150 (MQ=255)
                                                                                           aaaTGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAAc                                                            <  2:4163472/150‑1 (MQ=255)
                                                                                            aaTGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACt                                                           >  1:1000999/1‑150 (MQ=255)
                                                                                            aaTGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACt                                                           >  1:287983/1‑150 (MQ=255)
                                                                                                 aGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGgcgc                                                      <  1:704546/150‑1 (MQ=255)
                                                                                                 aGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGgcgc                                                      <  2:4035704/150‑1 (MQ=255)
                                                                                                     taTTCAGGACGATATCCGTCTTGCCGAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGTAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCct                                                  <  1:4550807/149‑1 (MQ=255)
                                                                                                     aaTTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCct                                                  <  2:408900/150‑1 (MQ=255)
                                                                                                      aTTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCctg                                                 >  1:3030860/1‑150 (MQ=255)
                                                                                                         cAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCctgctg                                              >  1:3585506/1‑150 (MQ=255)
                                                                                                         cAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCctgctg                                              >  1:1148039/1‑150 (MQ=255)
                                                                                                         cAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCctgctg                                              <  2:2234345/150‑1 (MQ=255)
                                                                                                         cAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTACCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCcggctg                                              <  1:2514595/150‑1 (MQ=255)
                                                                                                             aCGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATAGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGa                                          >  1:265485/1‑150 (MQ=255)
                                                                                                                atatCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCa                                       <  2:29498/150‑1 (MQ=255)
                                                                                                                 tatCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCaa                                      <  1:3529527/150‑1 (MQ=255)
                                                                                                                 tatCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCaa                                      >  1:1770553/1‑150 (MQ=255)
                                                                                                                    ccGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAgg                                   <  1:1055613/150‑1 (MQ=255)
                                                                                                                    ccGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAgg                                   <  2:636726/150‑1 (MQ=255)
                                                                                                                    ccGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAgg                                   >  1:1271224/1‑150 (MQ=255)
                                                                                                                     cGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGc                                  >  1:2727570/1‑150 (MQ=255)
                                                                                                                      gTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCa                                 <  2:722649/150‑1 (MQ=255)
                                                                                                                      gTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCa                                 >  1:369518/1‑150 (MQ=255)
                                                                                                                       tCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGTAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCAt                                <  2:2531334/150‑1 (MQ=255)
                                                                                                                       tCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCAt                                >  2:454222/1‑150 (MQ=255)
                                                                                                                         aTGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGTGACCAAAGGCATCa                              >  2:2818534/1‑150 (MQ=255)
                                                                                                                         aTGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGTGACCAAAGGCATCa                              >  2:1775341/1‑150 (MQ=255)
                                                                                                                          tGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCAc                             >  2:1796825/1‑150 (MQ=255)
                                                                                                                           gCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCAcc                            <  1:3343401/150‑1 (MQ=255)
                                                                                                                               aGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATg                        >  2:4384722/1‑150 (MQ=255)
                                                                                                                                     gcTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATcc                  <  1:4338393/150‑1 (MQ=255)
                                                                                                                                      cTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCa                 >  2:2182371/1‑150 (MQ=255)
                                                                                                                                       tGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAg                >  2:3755061/1‑150 (MQ=255)
                                                                                                                                       tGGGATCTGGCGATCGCTTCCCAAGGCTTTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCTGCGGATGCTCTCACCAAACTGGCTCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAg                <  2:2271676/150‑1 (MQ=255)
                                                                                                                                       tGGGATCTGGCGATAGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAg                <  1:3365641/150‑1 (MQ=255)
                                                                                                                                        gggATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGa               >  2:3671017/1‑150 (MQ=255)
                                                                                                                                         ggATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGAt              >  1:571512/1‑150 (MQ=255)
                                                                                                                                             cttgCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAATGCATCACCCATGAAATCCAGATCCCg          <  1:3322641/147‑1 (MQ=255)
                                                                                                                                                 cGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGata      <  1:646430/150‑1 (MQ=255)
                                                                                                                                                  gATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGatat     >  2:479450/1‑150 (MQ=255)
                                                                                                                                                     cGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATAtctc  <  1:2211368/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTAACGACCTGAAAATCAAGGTTGAAGATCAACTGGTTCACGCTCACTTCGAAGCGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTC  >  CP009273/4278697‑4278995

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: