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,3840CA99.1% 692.2 / ‑2.0 213H208N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (109/102);  total (109/104)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 213 total aligned reads displayed.

GGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGA  >  CP009273/4124238‑4124529
                                                                                                                                                  |                                                                                                                                                 
gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGCTGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAAt                                                                                                                                                 >  2:2452981/1‑149 (MQ=255)
gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAAt                                                                                                                                                 <  1:1428776/149‑1 (MQ=255)
gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAAt                                                                                                                                                 <  2:1072352/149‑1 (MQ=255)
gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAAt                                                                                                                                                 >  1:3978858/1‑149 (MQ=255)
 ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCTTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                <  2:1447243/149‑1 (MQ=255)
 ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                <  1:3609295/149‑1 (MQ=255)
 ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                >  2:2963562/1‑149 (MQ=255)
 ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                <  2:3495159/149‑1 (MQ=255)
 ggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATc                                                                                                                                                <  1:958938/149‑1 (MQ=255)
   ccGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCg                                                                                                                                              <  2:1020027/149‑1 (MQ=255)
    cGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCgg                                                                                                                                             >  2:2245632/1‑149 (MQ=255)
      aCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGaa                                                                                                                                           >  1:1433078/1‑149 (MQ=255)
        cTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgc                                                                                                                                         <  2:63349/149‑1 (MQ=255)
          gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcgc                                                                                                                                       >  1:748090/1‑149 (MQ=255)
          gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcgc                                                                                                                                       >  1:3666809/1‑149 (MQ=255)
          gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAgcgc                                                                                                                                       <  2:704074/149‑1 (MQ=255)
           tttATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCt                                                                                                                                      >  1:3442620/1‑149 (MQ=255)
           tttATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCt                                                                                                                                      >  1:1125217/1‑149 (MQ=255)
              aTCCTCGCGGGTACCGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGc                                                                                                                                   >  2:2522461/1‑149 (MQ=255)
               tCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCg                                                                                                                                  >  2:3895254/1‑149 (MQ=255)
                  tCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGaaa                                                                                                                               <  2:225118/149‑1 (MQ=255)
                  tCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGaaa                                                                                                                               <  2:3786034/149‑1 (MQ=255)
                      gggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCAc                                                                                                                           <  1:3935969/149‑1 (MQ=255)
                         tAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGc                                                                                                                        <  1:1702394/149‑1 (MQ=255)
                         tAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGc                                                                                                                        <  1:445075/149‑1 (MQ=255)
                              tGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGt                                                                                                                   <  2:2413414/149‑1 (MQ=255)
                                  gcTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCaa                                                                                                               >  2:2571710/1‑149 (MQ=255)
                                  gcTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCaa                                                                                                               >  1:464802/1‑149 (MQ=255)
                                   cTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAc                                                                                                              <  2:3548289/149‑1 (MQ=255)
                                         caCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGAt                                                                                                        <  2:3853499/148‑1 (MQ=255)
                                         aaCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGAt                                                                                                        >  2:798744/1‑149 (MQ=255)
                                           cTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATgg                                                                                                      >  1:1537883/1‑149 (MQ=255)
                                                gCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                                                                                 >  2:2015535/1‑149 (MQ=255)
                                                  ttCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                                                                               >  1:87934/1‑149 (MQ=255)
                                                   tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAtt                                                                                              <  1:3673225/149‑1 (MQ=255)
                                                   tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGAGAAAGCACCGCTGTGTGCAACCGAGATGGGTCTGAtt                                                                                              >  2:2420306/1‑149 (MQ=255)
                                                     cGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTa                                                                                            <  2:3213392/149‑1 (MQ=255)
                                                          cTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTa                                                                                       >  2:2104667/1‑149 (MQ=255)
                                                           ttCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTaa                                                                                      >  2:219094/1‑149 (MQ=255)
                                                           ttCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTaa                                                                                      >  2:1627911/1‑149 (MQ=255)
                                                            tCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAc                                                                                     <  1:661755/149‑1 (MQ=255)
                                                              ggTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAccc                                                                                   <  1:3774863/149‑1 (MQ=255)
                                                              ggTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAccc                                                                                   <  2:1536937/149‑1 (MQ=255)
                                                              ggTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAccc                                                                                   <  2:3585384/149‑1 (MQ=255)
                                                               gTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCg                                                                                  >  1:1188557/1‑149 (MQ=255)
                                                                    ggTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAgg                                                                             <  1:3183922/149‑1 (MQ=255)
                                                                    ggTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAgg                                                                             <  1:1850216/149‑1 (MQ=255)
                                                                         cGGTCGTGAATACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCgg                                                                        >  1:3570334/1‑149 (MQ=255)
                                                                         cGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCgg                                                                        >  1:7545/1‑149 (MQ=255)
                                                                           gTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt                                                                      <  2:2325996/149‑1 (MQ=255)
                                                                               tGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                  >  2:3620942/1‑149 (MQ=255)
                                                                               tGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcaca                                                                  >  2:1207469/1‑149 (MQ=255)
                                                                                 aaGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgc                                                                >  1:3408956/1‑149 (MQ=255)
                                                                                     cGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggcg                                                            <  2:3258778/149‑1 (MQ=255)
                                                                                     cGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggcg                                                            <  2:3258887/149‑1 (MQ=255)
                                                                                           ggAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGc                                                      >  1:1827436/1‑149 (MQ=255)
                                                                                            gAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCt                                                     <  1:2828957/149‑1 (MQ=255)
                                                                                              aCCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCttt                                                   >  1:3998073/1‑149 (MQ=255)
                                                                                               ccGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTc                                                  <  2:1228626/149‑1 (MQ=255)
                                                                                               ccGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTc                                                  <  1:3083568/149‑1 (MQ=255)
                                                                                                  gATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGc                                               >  1:1450858/1‑149 (MQ=255)
                                                                                                   aTCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCg                                              >  2:2063580/1‑149 (MQ=255)
                                                                                                       ggATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                          >  2:2389670/1‑149 (MQ=255)
                                                                                                        gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcagc                                         >  1:2835810/1‑149 (MQ=255)
                                                                                                         aTGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcagca                                        >  1:1605633/1‑149 (MQ=255)
                                                                                                          tGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGcagcag                                       >  2:421183/1‑149 (MQ=255)
                                                                                                           gTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCTAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCagcagc                                      <  2:105690/149‑1 (MQ=255)
                                                                                                              aaCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                   <  1:1088658/149‑1 (MQ=255)
                                                                                                               aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                                  <  1:2122029/149‑1 (MQ=255)
                                                                                                               aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                                  <  2:1603216/149‑1 (MQ=255)
                                                                                                               aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                                  <  1:43335/149‑1 (MQ=255)
                                                                                                               aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATc                                  <  1:1705072/149‑1 (MQ=255)
                                                                                                                 tGGGGTGATGAAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcg                                >  1:2357798/1‑149 (MQ=255)
                                                                                                                        atgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGAcc                         >  2:1587904/1‑149 (MQ=255)
                                                                                                                         tgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                        <  1:3674794/149‑1 (MQ=255)
                                                                                                                         tgaAAAAGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                        <  2:1461936/149‑1 (MQ=255)
                                                                                                                           aaaaaGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTc                      <  1:1125933/149‑1 (MQ=255)
                                                                                                                           aaaaaGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTc                      <  1:1812397/149‑1 (MQ=255)
                                                                                                                             aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCgg                    >  2:2689253/1‑149 (MQ=255)
                                                                                                                             aaaGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCgg                    >  1:2323011/1‑149 (MQ=255)
                                                                                                                               aGCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCa                  >  1:2699266/1‑149 (MQ=255)
                                                                                                                                gCCTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCaa                 <  2:2872753/149‑1 (MQ=255)
                                                                                                                                  cTGGCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACa               >  2:1034651/1‑149 (MQ=255)
                                                                                                                                     gCTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg            >  1:2398348/1‑149 (MQ=255)
                                                                                                                                      cTGACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAAATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATggg           <  1:2522461/149‑1 (MQ=255)
                                                                                                                                        gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCa         <  2:2792306/149‑1 (MQ=255)
                                                                                                                                        gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCa         >  2:3102197/1‑149 (MQ=255)
                                                                                                                                        gACTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCa         >  2:545031/1‑149 (MQ=255)
                                                                                                                                         aCTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt        <  1:1884079/149‑1 (MQ=255)
                                                                                                                                         aCTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt        <  1:1929794/149‑1 (MQ=255)
                                                                                                                                         aCTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt        >  1:3279462/1‑149 (MQ=255)
                                                                                                                                         aCTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCAt        <  1:218450/149‑1 (MQ=255)
                                                                                                                                          cTCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATg       >  1:1470024/1‑149 (MQ=255)
                                                                                                                                           tCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGa      >  1:3440166/1‑149 (MQ=255)
                                                                                                                                           tCACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGa      >  2:3959997/1‑149 (MQ=255)
                                                                                                                                            cACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa     <  2:4004704/149‑1 (MQ=255)
                                                                                                                                            cACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa     <  2:4004615/149‑1 (MQ=255)
                                                                                                                                            cACCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGaa     <  1:2249494/149‑1 (MQ=255)
                                                                                                                                             aCCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAac    <  1:798446/149‑1 (MQ=255)
                                                                                                                                             aCCGTAATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAac    >  1:1959928/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGA  >  CP009273/4124238‑4124529

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