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% 1029.1 / 1.9 312H208N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (166/141);  total (168/144)
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.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 312 total aligned reads displayed.

TGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGAC  >  CP009273/4124237‑4124530
                                                                                                                                                   |                                                                                                                                                  
tGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAAt                                                                                                                                                  >  1:4386105/1‑150 (MQ=255)
 gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                 <  2:2345023/150‑1 (MQ=255)
 gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTAGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                 >  1:3225811/1‑150 (MQ=255)
  ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGCTGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATcc                                                                                                                                                >  2:3336347/1‑150 (MQ=255)
  ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATcc                                                                                                                                                >  2:1304611/1‑150 (MQ=255)
  ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATcc                                                                                                                                                <  2:638757/150‑1 (MQ=255)
   gCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCg                                                                                                                                               <  1:103815/150‑1 (MQ=255)
     cGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGa                                                                                                                                             >  1:747162/1‑150 (MQ=255)
     cGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGa                                                                                                                                             >  1:3462018/1‑150 (MQ=255)
      gACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGaa                                                                                                                                            >  1:15221/1‑150 (MQ=255)
        ccTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgc                                                                                                                                          >  2:3003444/1‑150 (MQ=255)
          tGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAgcgc                                                                                                                                        <  2:243358/150‑1 (MQ=255)
            tttATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTg                                                                                                                                      >  1:1808638/1‑150 (MQ=255)
              tATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGc                                                                                                                                    <  2:3274772/150‑1 (MQ=255)
               aTCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCg                                                                                                                                   >  1:512735/1‑150 (MQ=255)
               aTCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGGACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCg                                                                                                                                   <  1:1427522/150‑1 (MQ=255)
                      cgGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCAc                                                                                                                            >  1:4857929/1‑150 (MQ=255)
                      cgGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCAc                                                                                                                            >  2:842956/1‑150 (MQ=255)
                        ggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCg                                                                                                                          <  2:4487483/150‑1 (MQ=255)
                         gTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGc                                                                                                                         >  2:4567257/1‑150 (MQ=255)
                          tAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCt                                                                                                                        >  2:3831969/1‑150 (MQ=255)
                           aaCGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                       >  1:134081/1‑150 (MQ=255)
                                ggCGCTAGAAAACTCCGGCTTCAGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACAGGATCTGGATGTTATCTGGTGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGc                                                                                                                  <  1:81276/150‑1 (MQ=255)
                                    cTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAcc                                                                                                              >  1:2686446/1‑150 (MQ=255)
                                    cTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAcc                                                                                                              <  1:1248009/150‑1 (MQ=255)
                                        aaaaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgaga                                                                                                          >  1:3792328/1‑150 (MQ=255)
                                          aaCTCCGGCTTCCGTACCTTCGGTGTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGGTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGACGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATg                                                                                                        >  2:2147859/1‑150 (MQ=255)
                                              ccGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTc                                                                                                    <  1:1749096/150‑1 (MQ=255)
                                                ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                                                                                  >  1:2068138/1‑150 (MQ=255)
                                                ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                                                                                  >  1:1897447/1‑150 (MQ=255)
                                                ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                                                                                  <  1:1870021/150‑1 (MQ=255)
                                                ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                                                                                  <  1:4428988/150‑1 (MQ=255)
                                                   ttCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAtt                                                                                               >  1:5047560/1‑150 (MQ=255)
                                                   ttCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAtt                                                                                               >  1:4335427/1‑150 (MQ=255)
                                                    tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAttt                                                                                              <  1:899618/150‑1 (MQ=255)
                                                     ccGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTa                                                                                             >  2:941543/1‑150 (MQ=255)
                                                       gTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACg                                                                                           >  1:1652484/1‑150 (MQ=255)
                                                        tACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGt                                                                                          >  2:4669344/1‑150 (MQ=255)
                                                          ccTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTa                                                                                        <  1:852288/150‑1 (MQ=255)
                                                           cTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTaa                                                                                       <  2:4976952/150‑1 (MQ=255)
                                                             tCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAcc                                                                                     <  2:1583293/150‑1 (MQ=255)
                                                              cGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAccc                                                                                    <  2:1393174/150‑1 (MQ=255)
                                                               ggTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCg                                                                                   >  1:3335726/1‑150 (MQ=255)
                                                                  tttGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGNGCAACCGAGATGGGTCTGATTTACGTTAACCCGGaa                                                                                <  1:2673945/150‑1 (MQ=255)
                                                                   ttGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAg                                                                               >  2:4586053/1‑150 (MQ=255)
                                                                    tGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAgg                                                                              >  2:1696879/1‑150 (MQ=255)
                                                                        gCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCg                                                                          <  1:1703705/150‑1 (MQ=255)
                                                                        gCCGGGCGGGAAGACGTCGGGGAACCGGCTCTGGATGTTACCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGCGATGGGTCGGATTTACGTTAACCCGGAAGGCCCg                                                                          <  2:3923592/150‑1 (MQ=255)
                                                                          cGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGa                                                                        >  1:4183933/1‑150 (MQ=255)
                                                                          cGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGa                                                                        <  2:4151872/150‑1 (MQ=255)
                                                                               gTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                   <  1:1494762/150‑1 (MQ=255)
                                                                                 gAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgc                                                                 >  2:3145690/1‑150 (MQ=255)
                                                                                 gAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgc                                                                 <  2:4929840/150‑1 (MQ=255)
                                                                                    gACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggc                                                              <  2:3920599/150‑1 (MQ=255)
                                                                                    gACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggc                                                              <  1:3125495/150‑1 (MQ=255)
                                                                                       gTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGaa                                                           >  1:1259774/1‑150 (MQ=255)
                                                                                        tCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAAc                                                          >  2:4369396/1‑150 (MQ=255)
                                                                                          tGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCg                                                        <  2:2041079/150‑1 (MQ=255)
                                                                                             gAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCtt                                                     >  1:213751/1‑150 (MQ=255)
                                                                                              aaCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCttt                                                    >  2:2500119/1‑150 (MQ=255)
                                                                                               aCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTc                                                   <  2:2660036/150‑1 (MQ=255)
                                                                                               aCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTc                                                   <  1:4161927/150‑1 (MQ=255)
                                                                                               aCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTc                                                   <  1:3991433/150‑1 (MQ=255)
                                                                                                ccGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCtttat                                                  <  2:671247/150‑3 (MQ=255)
                                                                                                ccGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCt                                                  <  2:1183960/150‑1 (MQ=255)
                                                                                                 cGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTg                                                 >  1:2388473/1‑150 (MQ=255)
                                                                                                   gATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCg                                               <  2:2721358/150‑1 (MQ=255)
                                                                                                   gATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCg                                               >  2:178861/1‑150 (MQ=255)
                                                                                                    aTCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCgg                                              >  2:3540355/1‑150 (MQ=255)
                                                                                                       tGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                           >  1:2214738/1‑150 (MQ=255)
                                                                                                              tAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                    >  1:2252911/1‑150 (MQ=255)
                                                                                                               aaCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                                   <  1:665177/150‑1 (MQ=255)
                                                                                                                aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcagca                                         >  2:1458635/1‑143 (MQ=255)
                                                                                                                 cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcg                                 <  1:998344/150‑1 (MQ=255)
                                                                                                                  tGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgc                                <  1:2108683/150‑1 (MQ=255)
                                                                                                                   ggggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgcg                               >  2:5101270/1‑150 (MQ=255)
                                                                                                                   ggggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgcg                               <  1:3651966/150‑1 (MQ=255)
                                                                                                                   ggggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgcg                               <  2:769394/150‑1 (MQ=255)
                                                                                                                    gggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgcgc                              >  1:2948968/1‑150 (MQ=255)
                                                                                                                           gaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTc                       <  2:1376621/150‑1 (MQ=255)
                                                                                                                             aaaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCgg                     >  2:2213809/1‑150 (MQ=255)
                                                                                                                             aaaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCgg                     <  2:4504402/150‑1 (MQ=255)
                                                                                                                              aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGc                    >  2:4071634/1‑150 (MQ=255)
                                                                                                                              aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGc                    >  1:2795545/1‑150 (MQ=255)
                                                                                                                               aaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCa                   <  2:225255/150‑1 (MQ=255)
                                                                                                                                  ccTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACa                >  2:2468395/1‑150 (MQ=255)
                                                                                                                                   cTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACAt               >  2:3744928/1‑150 (MQ=255)
                                                                                                                                   cTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACAt               <  2:2893188/150‑1 (MQ=255)
                                                                                                                                   cTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACAt               <  2:5157930/150‑1 (MQ=255)
                                                                                                                                    tgtcTGACTCACCTTAATCCGGAAGCGCTGGAGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTAAGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATg              <  2:3506323/147‑1 (MQ=255)
                                                                                                                                     ggCTGACTCACCGTACTCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg             >  2:35832/1‑150 (MQ=255)
                                                                                                                                     ggCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg             >  1:3163371/1‑150 (MQ=255)
                                                                                                                                     ggCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg             <  1:247277/150‑1 (MQ=255)
                                                                                                                                         gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt         <  2:3225721/150‑1 (MQ=255)
                                                                                                                                         gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt         <  2:701270/150‑1 (MQ=255)
                                                                                                                                            tCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa      >  1:426673/1‑150 (MQ=255)
                                                                                                                                              aCCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacg    >  1:1983195/1‑150 (MQ=255)
                                                                                                                                                cGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacgac  >  1:1611741/1‑150 (MQ=255)
                                                                                                                                                cGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacgac  >  2:2828505/1‑150 (MQ=255)
                                                                                                                                                cGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacgac  >  1:2887761/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGAC  >  CP009273/4124237‑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: