Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,124,384 C→A 100% H208N (CAT→AAT)  katG → catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,124,3840CA98.7% 1012.6 / 2.3 318H208N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (161/153);  total (163/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 319 total aligned reads displayed.

CCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGAC  >  CP009273/4124235‑4124530
                                                                                                                                                     |                                                                                                                                                  
ccTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGACGACGGCTGGGAACCGGATCTGGATGTTAACGGGGGTGATGACAAAGCCTGGCTGACTCACCGTc                                                                                                                                                    >  2:1669118/1‑150 (MQ=255)
  tGGGCCGACCTGTTTATCCTCTCGGGTAAAGTGGCGCAAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAAt                                                                                                                                                  >  2:3656245/1‑150 (MQ=255)
   gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                 <  2:1341937/150‑1 (MQ=255)
    ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATcc                                                                                                                                                <  1:4580057/150‑1 (MQ=255)
     gCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCg                                                                                                                                               >  1:1196708/1‑150 (MQ=255)
      ccGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCgg                                                                                                                                              >  2:261414/1‑150 (MQ=255)
       cGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGTCTGACTCACCGTAATCCGGa                                                                                                                                             >  1:5241588/1‑150 (MQ=255)
       cGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGa                                                                                                                                             >  2:1453551/1‑150 (MQ=255)
        gACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGaa                                                                                                                                            <  2:3347197/150‑1 (MQ=255)
         aCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTACTCCGGAAg                                                                                                                                           >  2:1189664/1‑150 (MQ=255)
          ccTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAgc                                                                                                                                          <  1:3646331/150‑1 (MQ=255)
           cTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGCTGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcg                                                                                                                                         >  1:2019845/1‑150 (MQ=255)
            tGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcgc                                                                                                                                        <  1:1306639/150‑1 (MQ=255)
             gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCt                                                                                                                                       >  1:2079066/1‑150 (MQ=255)
             gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCt                                                                                                                                       <  1:2704805/150‑1 (MQ=255)
                 aTCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCTCCGTCATCCGGAAGCGCTGGCg                                                                                                                                   >  2:3029310/1‑150 (MQ=255)
                  tCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGGCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGa                                                                                                                                  >  1:4793301/1‑150 (MQ=255)
                    cTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGaaa                                                                                                                                >  1:3314904/1‑150 (MQ=255)
                         gggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCAcc                                                                                                                           >  2:2403615/1‑150 (MQ=255)
                                  ggCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCAGGAAGCGCTGGCGAAAGCACCGCTGGGTGc                                                                                                                  >  2:4810256/1‑150 (MQ=255)
                                   gcgcTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCa                                                                                                                 <  1:4755092/150‑1 (MQ=255)
                                      cTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAcc                                                                                                              <  2:4355206/150‑1 (MQ=255)
                                       tAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCg                                                                                                             >  1:2601781/1‑150 (MQ=255)
                                        aGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCga                                                                                                            <  1:4308167/150‑1 (MQ=255)
                                          aaaaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgaga                                                                                                          >  2:2938134/1‑150 (MQ=255)
                                               tCCGGCTTCCGTACCTTCGGGTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGt                                                                                                     <  2:255440/150‑1 (MQ=255)
                                                 cGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCt                                                                                                   >  2:4153943/1‑150 (MQ=255)
                                                  ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                                                                                  >  1:816591/1‑150 (MQ=255)
                                                  ggCTTCCGTACCTTCGGTTTTGGTGCCGGACGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTTGCTGACTCACCGTAATCCTGAAGCGCTGTCGAAAGCACCGCTGGGTGCAACCTAGATGGGTATg                                                                                                  >  1:741706/1‑150 (MQ=255)
                                                      tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAttt                                                                                              <  1:2337755/150‑1 (MQ=255)
                                                       ccGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTa                                                                                             >  2:4971530/1‑150 (MQ=255)
                                                        cGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTAc                                                                                            <  1:4536955/150‑1 (MQ=255)
                                                         gTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACg                                                                                           >  2:5279580/1‑150 (MQ=255)
                                                         gTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACg                                                                                           >  1:1858259/1‑150 (MQ=255)
                                                            ccTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTa                                                                                        >  2:3842539/1‑150 (MQ=255)
                                                             cTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTaa                                                                                       >  1:41691/1‑150 (MQ=255)
                                                                cGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAccc                                                                                    <  1:1079767/150‑1 (MQ=255)
                                                                 ggTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCg                                                                                   >  2:1653864/1‑150 (MQ=255)
                                                                     ttGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAg                                                                               <  1:4648133/150‑1 (MQ=255)
                                                                      tGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAgg                                                                              <  2:1434282/150‑1 (MQ=255)
                                                                        gTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAATAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGcc                                                                            <  1:1345833/150‑1 (MQ=255)
                                                                         tGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGccc                                                                           >  1:3442105/1‑150 (MQ=255)
                                                                         tGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGccc                                                                           <  1:4863556/150‑1 (MQ=255)
                                                                             ggTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt                                                                       >  1:3950492/1‑150 (MQ=255)
                                                                              gTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATc                                                                      >  2:1807112/1‑150 (MQ=255)
                                                                                 gTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                   <  1:1644944/150‑1 (MQ=255)
                                                                                 gTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                   >  1:2415207/1‑150 (MQ=255)
                                                                                 gTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                   <  1:2740381/150‑1 (MQ=255)
                                                                                   gAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgc                                                                 >  2:4227870/1‑150 (MQ=255)
                                                                                   gAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgc                                                                 <  1:80339/150‑1 (MQ=255)
                                                                                      gACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggc                                                              <  2:3709911/150‑1 (MQ=255)
                                                                                         gTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGaa                                                           >  1:607293/1‑150 (MQ=255)
                                                                                          tCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAAc                                                          <  1:3837752/150‑1 (MQ=255)
                                                                                           cTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAAcc                                                         <  2:2191371/150‑1 (MQ=255)
                                                                                           cTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAAcc                                                         >  1:4233607/1‑150 (MQ=255)
                                                                                              ggAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCt                                                      <  1:3174797/150‑1 (MQ=255)
                                                                                                aaCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCttt                                                    <  1:44783/150‑1 (MQ=255)
                                                                                                aaCCGGATCTGGATGTTAACTGGGGTGAAGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCttt                                                    <  2:3517775/150‑1 (MQ=255)
                                                                                                      aTCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCgg                                              <  1:5119134/150‑1 (MQ=255)
                                                                                                        cTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgca                                            <  1:1017528/150‑1 (MQ=255)
                                                                                                         tGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                           >  1:867349/1‑150 (MQ=255)
                                                                                                          ggATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcagc                                          >  1:2976348/1‑150 (MQ=255)
                                                                                                          ggATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcagc                                          >  2:1772555/1‑150 (MQ=255)
                                                                                                             tGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCagcagc                                       >  1:2565839/1‑150 (MQ=255)
                                                                                                             tGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTATGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCagcagc                                       <  2:1542295/150‑1 (MQ=255)
                                                                                                             tGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAACGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCagcagc                                       >  2:3946910/1‑150 (MQ=255)
                                                                                                              gTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCt                                      <  2:2165545/150‑1 (MQ=255)
                                                                                                               ttAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTa                                     <  1:1490749/150‑1 (MQ=255)
                                                                                                               ttAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTa                                     <  1:2947407/150‑1 (MQ=255)
                                                                                                                tAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                    >  2:768464/1‑150 (MQ=255)
                                                                                                                  aCTGGGGTGCTGAAAAAGCCTGGCTGACTCACCGGAATCCGGTAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTTCGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATcc                                  <  1:2701274/150‑1 (MQ=255)
                                                                                                                  aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATcc                                  <  2:4614014/150‑1 (MQ=255)
                                                                                                                   cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcg                                 <  2:735633/150‑1 (MQ=255)
                                                                                                                   cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcg                                 >  1:3717009/1‑150 (MQ=255)
                                                                                                                   cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcg                                 <  1:996720/150‑1 (MQ=255)
                                                                                                                    tGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgc                                >  1:771414/1‑150 (MQ=255)
                                                                                                                     ggggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgcg                               >  2:3186530/1‑150 (MQ=255)
                                                                                                                     ggggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGACTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgcg                               <  1:1510284/150‑1 (MQ=255)
                                                                                                                       ggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgcgcg                             <  2:3297860/150‑1 (MQ=255)
                                                                                                                        gTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGa                            >  1:2505659/1‑150 (MQ=255)
                                                                                                                        gTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGa                            >  1:3792480/1‑150 (MQ=255)
                                                                                                                           atgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                         <  2:1789413/150‑1 (MQ=255)
                                                                                                                           atgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                         >  1:655674/1‑150 (MQ=255)
                                                                                                                            tgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCtt                        >  2:1991978/1‑150 (MQ=255)
                                                                                                                             gaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTc                       <  2:2771019/150‑1 (MQ=255)
                                                                                                                               aaaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCgg                     >  1:2639859/1‑150 (MQ=255)
                                                                                                                                  aGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCaa                  >  1:3740544/1‑150 (MQ=255)
                                                                                                                                        gCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATggg            <  2:3094657/150‑1 (MQ=255)
                                                                                                                                        gCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATggg            >  1:1430347/1‑150 (MQ=255)
                                                                                                                                         cTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGc           <  2:2079066/150‑1 (MQ=255)
                                                                                                                                          tGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCa          >  1:3051114/1‑150 (MQ=255)
                                                                                                                                          tGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCa          >  2:2507273/1‑150 (MQ=255)
                                                                                                                                           gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt         <  1:1453551/150‑1 (MQ=255)
                                                                                                                                            aCTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATg        <  1:303102/150‑1 (MQ=255)
                                                                                                                                              tCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa      >  1:2971801/1‑150 (MQ=255)
                                                                                                                                              tCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa      <  2:816591/150‑1 (MQ=255)
                                                                                                                                               cACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAac     <  1:497645/150‑1 (MQ=255)
                                                                                                                                               cACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAac     <  1:5171753/150‑1 (MQ=255)
                                                                                                                                                aCCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacg    >  2:3204256/1‑150 (MQ=255)
                                                                                                                                                  cGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacgac  <  1:3490780/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGAC  >  CP009273/4124235‑4124530

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: