Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,125,407 C→A 100% H549N (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,125,4070CA95.2% 908.4 / 27.3 304H549N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (4/11);  new base A (143/144);  total (147/157)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 304 total aligned reads displayed.

TGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCG  >  CP009273/4125260‑4125551
                                                                                                                                                   |                                                                                                                                                
 gggATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATg                                                                                                                                               <  2:1953344/150‑1 (MQ=255)
   gATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTa                                                                                                                                             >  2:2952571/1‑150 (MQ=255)
    aTGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGTATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTAc                                                                                                                                            <  2:4221083/150‑1 (MQ=255)
    aTGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTAc                                                                                                                                            >  2:3122287/1‑150 (MQ=255)
        gAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGtt                                                                                                                                        >  1:3184579/1‑150 (MQ=255)
         aaCGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGttt                                                                                                                                       <  1:4800307/150‑1 (MQ=255)
          aCGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTg                                                                                                                                      <  1:208451/150‑1 (MQ=255)
           cgccgcAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTgc                                                                                                                                     <  1:2492880/150‑1 (MQ=255)
           cgccgcAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTgc                                                                                                                                     <  1:3690356/150‑1 (MQ=255)
             ccgcAGCCGTTCGTGCTCTGCCTTTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTgcgc                                                                                                                                   <  2:2639477/150‑1 (MQ=255)
              cgcAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGcc                                                                                                                                  >  1:2263881/1‑150 (MQ=255)
               gcATCCGTTAGTGCTCTGCCTGTTCATGAGAAAATACAGAAAGAGTCTTGTAAAGCCTCGCTGTCGGATATAATAGTGCTTGCTGGTGTTGTTTGAGTTGATAAATCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGAGCCg                                                                                                                                 >  2:2222509/1‑150 (MQ=255)
               gcAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCg                                                                                                                                 <  2:1845075/150‑1 (MQ=255)
                cAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCgg                                                                                                                                >  1:4611576/1‑150 (MQ=255)
                   ccGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTc                                                                                                                             <  1:2361779/150‑1 (MQ=255)
                    cGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTcg                                                                                                                            >  2:1313185/1‑150 (MQ=255)
                        cGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGtt                                                                                                                        <  1:4447835/150‑1 (MQ=255)
                          tGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGa                                                                                                                      <  1:5138600/150‑1 (MQ=255)
                           gCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGAt                                                                                                                     <  1:1507337/150‑1 (MQ=255)
                           gCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGGTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGAt                                                                                                                     >  1:4403243/1‑150 (MQ=255)
                            ctctGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATg                                                                                                                    <  1:4196803/150‑1 (MQ=255)
                             tctGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATgc                                                                                                                   <  2:4107511/150‑1 (MQ=255)
                                gCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGcgcg                                                                                                                >  2:1383251/1‑150 (MQ=255)
                                  cTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTc                                                                                                              <  2:1140426/150‑1 (MQ=255)
                                   tGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCa                                                                                                             <  1:221993/150‑1 (MQ=255)
                                     ttCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAgg                                                                                                           <  1:2832416/150‑1 (MQ=255)
                                     ttCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAgg                                                                                                           <  1:90423/150‑1 (MQ=255)
                                         ggAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCa                                                                                                       >  2:4053760/1‑150 (MQ=255)
                                           agaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGa                                                                                                     >  2:4181393/1‑150 (MQ=255)
                                              aaaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTg                                                                                                  >  1:770421/1‑150 (MQ=255)
                                                aTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                                                >  2:1606504/1‑150 (MQ=255)
                                                aTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                                                <  1:527589/150‑1 (MQ=255)
                                                 tCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACa                                                                                               <  2:1770634/150‑1 (MQ=255)
                                                  ccAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACAt                                                                                              <  1:2041327/150‑1 (MQ=255)
                                                  ccAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACAt                                                                                              >  1:1367731/1‑150 (MQ=255)
                                                   cAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACAtt                                                                                             <  2:2107298/150‑1 (MQ=255)
                                                    aGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTg                                                                                            >  1:3908426/1‑150 (MQ=255)
                                                    aGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTg                                                                                            >  2:1810283/1‑150 (MQ=255)
                                                    aGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTg                                                                                            >  2:4614493/1‑150 (MQ=255)
                                                     gAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGATTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTga                                                                                           <  2:1658022/150‑1 (MQ=255)
                                                      aaaGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTgag                                                                                          >  1:3577067/1‑150 (MQ=255)
                                                      aaaGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTgag                                                                                          <  2:3228421/150‑1 (MQ=255)
                                                       aaGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTgaga                                                                                         >  2:5017462/1‑150 (MQ=255)
                                                          agTCTGGTAAAGCCTCGCTGGCGGATATCCTAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGt                                                                                      >  2:977543/1‑150 (MQ=255)
                                                           gTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGtt                                                                                     >  1:1116125/1‑150 (MQ=255)
                                                           gTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGtt                                                                                     <  2:914078/150‑1 (MQ=255)
                                                            tCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGttt                                                                                    <  1:2152251/150‑1 (MQ=255)
                                                              tGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGa                                                                                  >  2:3140431/1‑150 (MQ=255)
                                                                gTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgc                                                                                >  2:107327/1‑150 (MQ=255)
                                                                  aaaGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCGGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                              >  2:2705515/1‑150 (MQ=255)
                                                                  aaaGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                              >  1:1775953/1‑150 (MQ=255)
                                                                  aaaGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGAGGTTTGAgctg                                                                              >  2:3899065/1‑150 (MQ=255)
                                                                      ccTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTgg                                                                          >  2:4734216/1‑150 (MQ=255)
                                                                          gCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGcc                                                                      >  1:185554/1‑150 (MQ=255)
                                                                               cGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTg                                                                 <  1:4591182/150‑1 (MQ=255)
                                                                                ggATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGc                                                                <  2:757062/150‑1 (MQ=255)
                                                                                 gATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAACGCCGCAAGCGCCTCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCt                                                               <  2:4805159/150‑1 (MQ=255)
                                                                                  atatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTg                                                              >  2:4496279/1‑150 (MQ=255)
                                                                                   tatCATAGTGCTGGGTGGTGTGGTTGGTGGTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGa                                                             >  2:2823977/1‑150 (MQ=255)
                                                                                    atCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATCCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGAc                                                            <  1:3872679/150‑1 (MQ=255)
                                                                                     tCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACg                                                           <  2:4898383/150‑1 (MQ=255)
                                                                                       aTAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGt                                                         <  2:3068286/150‑1 (MQ=255)
                                                                                        tAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGtt                                                        <  1:3744730/150‑1 (MQ=255)
                                                                                         aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGTAGCCAATTGCTGACGGttt                                                       >  2:146744/1‑150 (MQ=255)
                                                                                         aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGttt                                                       >  1:3457575/1‑150 (MQ=255)
                                                                                         aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGttt                                                       >  2:233221/1‑150 (MQ=255)
                                                                                          gTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTc                                                      >  2:3447392/1‑150 (MQ=255)
                                                                                            gctggctgGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCg                                                    <  1:3850748/150‑1 (MQ=255)
                                                                                             ctggctggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGt                                                   <  2:3387802/150‑1 (MQ=255)
                                                                                             ctggctggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGt                                                   >  2:2369142/1‑150 (MQ=255)
                                                                                             ctggctggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGt                                                   <  2:2876368/150‑1 (MQ=255)
                                                                                                gctggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAAc                                                <  1:3721155/150‑1 (MQ=255)
                                                                                                   ggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTAt                                             <  2:4939843/150‑1 (MQ=255)
                                                                                                   ggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTAt                                             <  1:3814017/150‑1 (MQ=255)
                                                                                                    gtgtGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATc                                            <  1:1226485/150‑1 (MQ=255)
                                                                                                       tggttggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATcgcg                                         >  1:3964595/1‑150 (MQ=255)
                                                                                                       tggttggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATcgcg                                         <  1:4318079/150‑1 (MQ=255)
                                                                                                       tggttggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATcgcg                                         <  2:105559/150‑1 (MQ=255)
                                                                                                        ggttggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCgcgc                                        >  1:3279814/1‑150 (MQ=255)
                                                                                                          ttggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTc                                      >  1:2161493/1‑150 (MQ=255)
                                                                                                           tggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCg                                     >  2:820267/1‑150 (MQ=255)
                                                                                                           tggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCg                                     >  1:2917376/1‑150 (MQ=255)
                                                                                                              tgtTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCt                                  <  2:550276/150‑1 (MQ=255)
                                                                                                              tgtTGAGAAAGCCGCAAGCGCCGCAGGTGTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCt                                  <  1:23078/150‑1 (MQ=255)
                                                                                                                ttGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCACTTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTgg                                >  2:1902922/1‑150 (MQ=255)
                                                                                                                    gaAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGt                            >  2:621109/1‑150 (MQ=255)
                                                                                                                     aaaGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGtt                           >  2:3956085/1‑150 (MQ=255)
                                                                                                                      aaGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGttt                          >  2:2276670/1‑150 (MQ=255)
                                                                                                                      aaGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGttt                          <  1:3122287/150‑1 (MQ=255)
                                                                                                                              aGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCaccacc                  >  1:1933312/1‑150 (MQ=255)
                                                                                                                                 gccgcAGGAGTGAGCATTCCTGTCCCGTGTGCGCCGGGTCGCGTGGATGCGCGTCAGGATCAGACGGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAg               <  2:4374403/150‑1 (MQ=255)
                                                                                                                                  ccgcAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGANt              <  1:2703040/150‑1 (MQ=255)
                                                                                                                                  ccgcAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGt              <  1:4970455/150‑1 (MQ=255)
                                                                                                                                   cgcAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTc             >  1:2786920/1‑150 (MQ=255)
                                                                                                                                      aGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCAct          >  2:2101829/1‑150 (MQ=255)
                                                                                                                                         tttGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCActgct       <  1:859216/150‑1 (MQ=255)
                                                                                                                                           tGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGa     >  2:3294169/1‑150 (MQ=255)
                                                                                                                                           agAGCATTAATGTACCGATTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGATGGAGCCAAATGCTGACGGTTTCCGTAAATATCGCGCTCGTCTGGACGTTTCCACCACCTAGTCACTGCTGa     >  2:3858269/2‑150 (MQ=255)
                                                                                                                                              gCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCg  >  2:3655320/1‑150 (MQ=255)
                                                                                                                                              gCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCg  <  2:1558032/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCG  >  CP009273/4125260‑4125551

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: