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,3840CA97.7% 538.3 / 3.7 172H208N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (92/76);  total (94/78)
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.67e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 172 total aligned reads displayed.

GCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACG  >  CP009273/4124240‑4124528
                                                                                                                                                |                                                                                                                                                
gCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                               <  2:1944608/148‑1 (MQ=255)
gCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                               <  2:1141865/148‑1 (MQ=255)
gCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                               <  1:3158523/148‑1 (MQ=255)
    aCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGa                                                                                                                                           >  2:1986935/1‑148 (MQ=255)
    aCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGa                                                                                                                                           >  2:2313813/1‑148 (MQ=255)
    aCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGa                                                                                                                                           <  1:2123735/148‑1 (MQ=255)
    aCCTGTTTATCCTCGCGGGGAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGa                                                                                                                                           <  2:1172740/148‑1 (MQ=255)
      cTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAg                                                                                                                                         >  1:2503302/1‑148 (MQ=255)
       tGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgc                                                                                                                                        >  1:2663063/1‑148 (MQ=255)
        gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcg                                                                                                                                       >  2:1422883/1‑148 (MQ=255)
        gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcg                                                                                                                                       >  2:1425109/1‑148 (MQ=255)
        gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcg                                                                                                                                       >  1:1714872/1‑148 (MQ=255)
        gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcg                                                                                                                                       >  2:3195536/1‑148 (MQ=255)
        gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTACGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcg                                                                                                                                       >  1:3026242/1‑148 (MQ=255)
                 cgcgGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGaaa                                                                                                                              <  1:499462/148‑1 (MQ=255)
                 cgcgGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGaaa                                                                                                                              <  1:2421773/148‑1 (MQ=255)
                     ggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCAc                                                                                                                          >  2:562380/1‑148 (MQ=255)
                          cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                     >  1:2078088/1‑148 (MQ=255)
                          cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                     >  2:2981086/1‑148 (MQ=255)
                          cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                     >  1:3578310/1‑148 (MQ=255)
                          cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                     <  1:1590217/148‑1 (MQ=255)
                          cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGACTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTg                                                                                                                     >  2:1215407/1‑148 (MQ=255)
                                gcTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCa                                                                                                               <  2:1248445/148‑1 (MQ=255)
                                 cTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCaa                                                                                                              <  1:2560017/148‑1 (MQ=255)
                                   aGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAcc                                                                                                            >  1:1333153/1‑148 (MQ=255)
                                   aGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAcc                                                                                                            >  1:3651751/1‑148 (MQ=255)
                                   aGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAcc                                                                                                            >  1:779108/1‑148 (MQ=255)
                                     aaaaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCga                                                                                                          >  1:1203698/1‑148 (MQ=255)
                                      aaaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgag                                                                                                         <  2:105719/148‑1 (MQ=255)
                                       aaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgaga                                                                                                        >  1:934515/1‑148 (MQ=255)
                                       aaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgaga                                                                                                        <  2:2010699/148‑1 (MQ=255)
                                         cTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATg                                                                                                      <  2:1883745/148‑1 (MQ=255)
                                           ccGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGCTGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATggg                                                                                                    >  1:1613394/1‑148 (MQ=255)
                                             ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTc                                                                                                  <  1:1018965/148‑1 (MQ=255)
                                             ggCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTc                                                                                                  <  1:1019011/148‑1 (MQ=255)
                                              gCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCt                                                                                                 <  1:853518/148‑1 (MQ=255)
                                               cTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                                                                                >  2:2212853/1‑148 (MQ=255)
                                                ttCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGa                                                                                               <  1:2977630/148‑1 (MQ=255)
                                                 tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                                                                              >  1:497700/1‑148 (MQ=255)
                                                 tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                                                                              <  2:2040547/148‑1 (MQ=255)
                                                  ccGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAtt                                                                                             >  2:2648554/1‑148 (MQ=255)
                                                      aCCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACg                                                                                         <  1:2726886/148‑1 (MQ=255)
                                                        cTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGtt                                                                                       >  1:3802008/1‑148 (MQ=255)
                                                         ttCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTa                                                                                      <  1:2212853/148‑1 (MQ=255)
                                                          tCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTaa                                                                                     <  1:3599632/148‑1 (MQ=255)
                                                          tCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTaa                                                                                     <  2:1340063/148‑1 (MQ=255)
                                                                ttGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGa                                                                               >  1:3354138/1‑148 (MQ=255)
                                                                ttGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGa                                                                               <  1:191171/148‑1 (MQ=255)
                                                                    tGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGACGGc                                                                           <  1:276561/148‑1 (MQ=255)
                                                                     gCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGcc                                                                          <  1:292462/148‑1 (MQ=255)
                                                                     gCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGcc                                                                          <  2:3037803/148‑1 (MQ=255)
                                                                     gCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGcc                                                                          <  2:1009903/148‑1 (MQ=255)
                                                                          tCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt                                                                     >  1:1527144/1‑148 (MQ=255)
                                                                              gAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                 >  1:242260/1‑148 (MQ=255)
                                                                              gAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                 <  2:1113586/148‑1 (MQ=255)
                                                                                aGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgc                                                               >  1:3242800/1‑148 (MQ=255)
                                                                                  aCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcgg                                                             >  1:630644/1‑148 (MQ=255)
                                                                                   cGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggc                                                            >  2:1078204/1‑148 (MQ=255)
                                                                                    gTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggcg                                                           <  1:204750/148‑1 (MQ=255)
                                                                                         ggAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCg                                                      >  2:3037259/1‑148 (MQ=255)
                                                                                         ggAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCg                                                      >  2:3455350/1‑148 (MQ=255)
                                                                                            aCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCtt                                                   >  1:395683/1‑148 (MQ=255)
                                                                                             ccGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCttt                                                  >  2:2236572/1‑148 (MQ=255)
                                                                                              cGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTc                                                 <  2:3488154/148‑1 (MQ=255)
                                                                                                  tCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCg                                             >  1:1496426/1‑148 (MQ=255)
                                                                                                      gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                         >  2:1532294/1‑148 (MQ=255)
                                                                                                         gTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGTGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGcagcag                                      >  2:2913915/1‑148 (MQ=255)
                                                                                                         gTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGcagcag                                      >  2:2913870/1‑148 (MQ=255)
                                                                                                         gTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGcagcag                                      >  1:2389765/1‑148 (MQ=255)
                                                                                                          ttAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCTGCagcagc                                     >  1:3269410/1‑148 (MQ=255)
                                                                                                           tAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCt                                    <  2:850442/148‑1 (MQ=255)
                                                                                                           tAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCt                                    >  1:2018500/1‑148 (MQ=255)
                                                                                                             aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                  <  1:1422883/148‑1 (MQ=255)
                                                                                                             aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                  <  1:783962/148‑1 (MQ=255)
                                                                                                             aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                  <  1:783853/148‑1 (MQ=255)
                                                                                                              cttgtgTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                                 <  1:2105008/143‑1 (MQ=255)
                                                                                                                ggggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcg                               <  1:1804456/148‑1 (MQ=255)
                                                                                                                 gggTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAGTTACGTTAACCCGGAAGGCCCGGATCACAGGGGCGAACCGCTGTCTGCGGCAGCAGCTCTCcgc                              >  2:3520256/1‑148 (MQ=255)
                                                                                                                       tgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGAcc                        <  2:2031210/148‑1 (MQ=255)
                                                                                                                       tgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGAcc                        <  2:3578310/148‑1 (MQ=255)
                                                                                                                       tgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGAcc                        <  2:2078088/148‑1 (MQ=255)
                                                                                                                        gaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                       <  2:2669352/148‑1 (MQ=255)
                                                                                                                         aaaaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCtt                      <  1:1292600/148‑1 (MQ=255)
                                                                                                                         aaaaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCtt                      <  1:1292645/148‑1 (MQ=255)
                                                                                                                           aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCg                    >  2:514186/1‑148 (MQ=255)
                                                                                                                            aaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCgg                   >  2:3268567/1‑148 (MQ=255)
                                                                                                                             aGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGc                  >  2:3501679/1‑148 (MQ=255)
                                                                                                                             aGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGc                  >  2:1163795/1‑148 (MQ=255)
                                                                                                                             aGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGc                  >  2:2794419/1‑148 (MQ=255)
                                                                                                                              gCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCa                 >  1:3081663/1‑148 (MQ=255)
                                                                                                                                cTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAAc               <  2:3549952/148‑1 (MQ=255)
                                                                                                                                 tGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACa              <  1:3195536/148‑1 (MQ=255)
                                                                                                                                 tGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACa              >  1:1836559/1‑148 (MQ=255)
                                                                                                                                  ggCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCACCAt             >  1:3206773/1‑148 (MQ=255)
                                                                                                                                    cTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg           <  1:2637828/148‑1 (MQ=255)
                                                                                                                                    cTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg           <  2:2846900/148‑1 (MQ=255)
                                                                                                                                      gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGc         >  1:2114983/1‑148 (MQ=255)
                                                                                                                                         tCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATg      >  2:1827070/1‑148 (MQ=255)
                                                                                                                                             cGTAATCCGGAAGCGCTGGCGAAAGTACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacg  <  2:3676350/148‑1 (MQ=255)
                                                                                                                                             cGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacg  <  2:3802008/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACG  >  CP009273/4124240‑4124528

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: ‑