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,3840CA96.2% 621.5 / 14.3 209H208N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (122/79);  total (126/83)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 209 total aligned reads displayed.

GGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAA  >  CP009273/4124239‑4124526
                                                                                                                                                 |                                                                                                                                              
  ccGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATcc                                                                                                                                            <  2:211863/148‑1 (MQ=255)
  ccGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATcc                                                                                                                                            >  1:1407523/1‑148 (MQ=255)
        tGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgc                                                                                                                                      >  1:3305065/1‑148 (MQ=255)
         gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcg                                                                                                                                     <  1:1161499/148‑1 (MQ=255)
              tCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGc                                                                                                                                >  1:1778802/1‑148 (MQ=255)
                cTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTATTGGTGACGGTCGTGAATACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCAGGCTGACTCACCGTAATCCGGAAGCGCTGGCGa                                                                                                                              >  2:1692139/1‑148 (MQ=255)
                atCGAGGGTAACGTGGCGCTAGAAAACTCCGGCTTACGTAACTTCGGTTTTGGTTCCGGTCCTGAAGACGTCTGGGAAACGGATCTGGATGTTAACTGGGGTGAGGAAAAAGACTAGCTGACTCACAGTAATCCGGAAGACCTGGCGa                                                                                                                              >  2:1692072/2‑148 (MQ=255)
                  cgcgGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGaaa                                                                                                                            <  2:2124198/148‑1 (MQ=255)
                     gggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCa                                                                                                                         <  2:1325610/148‑1 (MQ=255)
                      ggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCAc                                                                                                                        <  1:814449/148‑1 (MQ=255)
                           cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAATCCTTGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                   >  1:3040296/1‑148 (MQ=255)
                           cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                   <  1:447517/148‑1 (MQ=255)
                           cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                   >  1:3041018/1‑148 (MQ=255)
                             tGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTggg                                                                                                                 >  2:930191/1‑148 (MQ=255)
                              ggCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGt                                                                                                                >  1:1461556/1‑148 (MQ=255)
                               gcgcTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTg                                                                                                               >  1:1520601/1‑148 (MQ=255)
                                 gcTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCa                                                                                                             <  1:1117270/148‑1 (MQ=255)
                                     gAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCg                                                                                                         <  2:1778802/148‑1 (MQ=255)
                                     gAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCg                                                                                                         <  2:606259/148‑1 (MQ=255)
                                      aaaaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCga                                                                                                        >  2:1707175/1‑148 (MQ=255)
                                        aaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgaga                                                                                                      <  2:1065031/148‑1 (MQ=255)
                                        aaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgaga                                                                                                      <  2:2427676/148‑1 (MQ=255)
                                          cTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATg                                                                                                    <  1:970351/148‑1 (MQ=255)
                                              ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTc                                                                                                <  2:376959/148‑1 (MQ=255)
                                                 ttCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGa                                                                                             <  1:1017441/148‑1 (MQ=255)
                                                  tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                                                                            >  2:1971317/1‑148 (MQ=255)
                                                  tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                                                                            >  2:671257/1‑148 (MQ=255)
                                                         cTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGtt                                                                                     >  1:2443030/1‑148 (MQ=255)
                                                             ggTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAcc                                                                                 <  2:813803/148‑1 (MQ=255)
                                                               ttttGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCg                                                                               >  2:3229275/1‑148 (MQ=255)
                                                               ttttGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCg                                                                               >  2:2584744/1‑148 (MQ=255)
                                                                    gTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAgg                                                                          <  1:1291949/148‑1 (MQ=255)
                                                                          gTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGa                                                                    <  1:2513123/148‑1 (MQ=255)
                                                                           tCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt                                                                   <  1:2113833/148‑1 (MQ=255)
                                                                           tCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt                                                                   >  2:1660011/1‑148 (MQ=255)
                                                                            cGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATc                                                                  >  1:1497511/1‑148 (MQ=255)
                                                                              tGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcac                                                                >  2:1281832/1‑148 (MQ=255)
                                                                              tGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcac                                                                >  2:1875262/1‑148 (MQ=255)
                                                                               gAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                               <  2:2081048/148‑1 (MQ=255)
                                                                                aaGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAg                                                              <  2:1472951/148‑1 (MQ=255)
                                                                                aaGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAg                                                              <  2:3312409/148‑1 (MQ=255)
                                                                                   aCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcgg                                                           <  2:84603/148‑1 (MQ=255)
                                                                                   aCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcgg                                                           >  1:2101263/1‑148 (MQ=255)
                                                                                       cTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGaa                                                       <  2:249600/148‑1 (MQ=255)
                                                                                           gAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGc                                                   <  2:1796805/148‑1 (MQ=255)
                                                                                            aaCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCt                                                  <  1:850313/148‑1 (MQ=255)
                                                                                                   tCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCg                                           >  2:493281/1‑148 (MQ=255)
                                                                                                   tCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCg                                           >  2:492759/1‑148 (MQ=255)
                                                                                                     tGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgc                                         <  1:345515/148‑1 (MQ=255)
                                                                                                      ggATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgca                                        >  2:1823894/1‑148 (MQ=255)
                                                                                                      ggATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgca                                        >  1:1225268/1‑148 (MQ=255)
                                                                                                      ggATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgca                                        >  1:2886028/1‑148 (MQ=255)
                                                                                                       gATTTTTACTGGTGTGATTAATAAGCCTGGCTGACTCACCGTAATCCGGTAGCGCTGGCGAAAGCACCGCTGTGTGCAAACGAGATGGGTCTTATTTACTTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGATTTCTGCGgcag                                       <  1:3133506/148‑1 (MQ=255)
                                                                                                       gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                       >  1:1953214/1‑148 (MQ=255)
                                                                                                       gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                       >  1:3509447/1‑148 (MQ=255)
                                                                                                       gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                       <  1:1365145/148‑1 (MQ=255)
                                                                                                          gTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGcagcag                                    >  2:2543100/1‑148 (MQ=255)
                                                                                                            tAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCt                                  >  2:1589440/1‑148 (MQ=255)
                                                                                                            tAACTGGGGTGATGAAAAAGCCGGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCt                                  <  2:3158992/148‑1 (MQ=255)
                                                                                                             aaCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTa                                 >  1:2414125/1‑148 (MQ=255)
                                                                                                              aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                >  1:3003616/1‑148 (MQ=255)
                                                                                                               cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                               >  2:2155274/1‑148 (MQ=255)
                                                                                                               cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                               >  1:2962885/1‑148 (MQ=255)
                                                                                                               cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                               >  1:3391062/1‑148 (MQ=255)
                                                                                                               cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                               >  2:1245540/1‑148 (MQ=255)
                                                                                                                  gggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgc                            >  1:2251169/1‑148 (MQ=255)
                                                                                                                  gggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcgc                            >  1:2251302/1‑148 (MQ=255)
                                                                                                                       atgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGAc                       >  2:1059200/1‑148 (MQ=255)
                                                                                                                       atgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGAc                       <  1:3070700/148‑1 (MQ=255)
                                                                                                                        tgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGAcc                      <  2:1032459/148‑1 (MQ=255)
                                                                                                                         gaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                     >  1:2927820/1‑148 (MQ=255)
                                                                                                                         gaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                     <  2:1972450/148‑1 (MQ=255)
                                                                                                                          aaaaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCtt                    >  1:2096519/1‑148 (MQ=255)
                                                                                                                           aaaaGCCTGGCTGACTCACCGGAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTc                   <  2:120358/148‑1 (MQ=255)
                                                                                                                            aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTGTCCGCGCGACCTTCg                  >  1:1812316/1‑148 (MQ=255)
                                                                                                                            aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCg                  >  1:600531/1‑148 (MQ=255)
                                                                                                                            aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAGTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCg                  >  2:1622328/1‑148 (MQ=255)
                                                                                                                             aaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCgg                 >  1:888794/1‑148 (MQ=255)
                                                                                                                              aGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGc                <  1:563255/148‑1 (MQ=255)
                                                                                                                               gCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCa               >  1:3380541/1‑148 (MQ=255)
                                                                                                                               gCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCa               <  2:3504308/148‑1 (MQ=255)
                                                                                                                               gCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCa               >  1:1334245/1‑148 (MQ=255)
                                                                                                                               gCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCa               >  1:3553106/1‑148 (MQ=255)
                                                                                                                                ccTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCaa              >  1:1251009/1‑148 (MQ=255)
                                                                                                                                ccTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCaa              >  1:668668/1‑148 (MQ=255)
                                                                                                                                  tGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACa            <  2:1820790/148‑1 (MQ=255)
                                                                                                                                   ggCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACAt           >  2:2101079/1‑148 (MQ=255)
                                                                                                                                    gCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATg          >  2:2865121/1‑148 (MQ=255)
                                                                                                                                     cTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg         <  2:3041018/148‑1 (MQ=255)
                                                                                                                                     cTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg         <  2:3040296/148‑1 (MQ=255)
                                                                                                                                       gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGc       >  2:741940/1‑148 (MQ=255)
                                                                                                                                        aCTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCa      >  1:99476/1‑148 (MQ=255)
                                                                                                                                         cTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt     <  2:3096822/148‑1 (MQ=255)
                                                                                                                                         cTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt     >  2:3513073/1‑148 (MQ=255)
                                                                                                                                         cTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt     >  2:3513322/1‑148 (MQ=255)
                                                                                                                                         cTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt     >  2:903192/1‑148 (MQ=255)
                                                                                                                                         cTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGGCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt     >  2:2956598/1‑148 (MQ=255)
                                                                                                                                           cACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGa   <  2:575092/148‑1 (MQ=255)
                                                                                                                                            aCCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa  >  1:2814467/1‑148 (MQ=255)
                                                                                                                                            aCCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa  <  1:830674/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAA  >  CP009273/4124239‑4124526

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: