Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,175,178 C→A 8.9% intergenic (‑293/‑94) katG ← / → PP_3669 catalase‑peroxidase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,175,1780CA 8.9% 747.7 / 45.3 235intergenic (‑293/‑94)katG/PP_3669catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (109/105);  new base A (8/13);  total (117/118)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 238 total aligned reads displayed.

CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAG  >  NC_002947/4175029‑4175326
                                                                                                                                                     |                                                                                                                                                    
      cTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                                <  1:2713381/150‑1 (MQ=255)
      cTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                                <  2:4823813/150‑1 (MQ=255)
       tGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTACCGGTCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                               <  1:5667884/150‑1 (MQ=255)
        gACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATg                                                                                                                                              >  1:5037014/1‑150 (MQ=255)
         ccgcgCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTCTCAGGCGCAGGCCTGTGCCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                             <  1:5877429/149‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                             <  2:1292290/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                             <  2:1420818/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                             <  2:3789310/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                             <  2:240196/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                             <  2:386975/150‑1 (MQ=255)
          cgcgCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAg                                                                                                                                            <  2:2239492/150‑1 (MQ=255)
           gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAgg                                                                                                                                           >  2:2829649/1‑150 (MQ=255)
           gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAgg                                                                                                                                           >  2:1082421/1‑150 (MQ=255)
            cgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGc                                                                                                                                          <  1:5411787/150‑1 (MQ=255)
               tGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCgg                                                                                                                                       >  2:321526/1‑150 (MQ=255)
                 cccTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGca                                                                                                                                     >  2:3534739/1‑150 (MQ=255)
                    tGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACAt                                                                                                                                  >  1:5539599/1‑150 (MQ=255)
                      tCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATcc                                                                                                                                >  2:4148953/1‑150 (MQ=255)
                      tCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATcc                                                                                                                                >  2:6760108/1‑150 (MQ=255)
                        gcgGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCaa                                                                                                                              >  1:4837800/1‑149 (MQ=255)
                            aaTTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCGTATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaa                                                                                                                          <  2:5798059/150‑1 (MQ=255)
                              ttgggggggCCACGATAAAGACTAGCCGAGGTGGGGCTGAGGTCTAATTGCGTGGCTGTTGGAGCCTGATTGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaagg                                                                                                                        <  2:1687647/145‑1 (MQ=255)
                              ttCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATTAGGCCGGCACATACATacaatg                                                                                                                        >  1:467754/1‑148 (MQ=255)
                              ttCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACCTCCATacaagg                                                                                                                        >  1:1714891/1‑150 (MQ=255)
                              ttCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaagg                                                                                                                        >  1:2577237/1‑150 (MQ=255)
                               tCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggc                                                                                                                       <  2:5039898/150‑1 (MQ=255)
                               tCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggc                                                                                                                       <  2:2168460/150‑1 (MQ=255)
                                cgcgGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggca                                                                                                                      >  2:154831/1‑150 (MQ=255)
                                cgcgGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggca                                                                                                                      <  1:6760108/150‑1 (MQ=255)
                                   ggggCCACGAGAACGACTAGCCGAGCTTGGGCAGAGGTCTAATAGCGTGGCGGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggcacaa                                                                                                                   <  2:3116836/150‑1 (MQ=255)
                                   ggggCCACGAGAAAGACTAGCCGTGCTTGGGCAGAGGTATAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTTTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggcacaa                                                                                                                   <  2:269959/150‑1 (MQ=255)
                                       ccACGAGAAAGACTAGCCGAGCTTGGGCTTTGGTATAATAGTGTGGCTGTTGGAGCGTGATAGGCTGGCTGTGTCTGGCGCAGGCCTTTATCGGCCTAATCGCGGTCAAGCCCGCGTTGAGGCCGGCACATCCATACAAGGCACAAGGcc                                                                                                               <  1:4866629/150‑1 (MQ=255)
                                         aCGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGcccc                                                                                                             <  2:6841767/150‑1 (MQ=255)
                                         aCGAGAAAGACTAGCCGAGCTGGGGCATAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCTGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGcccc                                                                                                             <  2:6435516/150‑1 (MQ=255)
                                             gaAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCAACGCGGTCAAGCCCGCGATGAGGCCGGCCCATCCATACAAGGCACAAGGCCCCgcgc                                                                                                         >  1:75217/1‑150 (MQ=255)
                                                aGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCAt                                                                                                      >  1:2123434/1‑150 (MQ=255)
                                                 gACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATg                                                                                                     >  2:2845304/1‑150 (MQ=255)
                                                  aCTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCTATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGc                                                                                                    <  1:3262292/150‑1 (MQ=255)
                                                    tAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtt                                                                                                  >  1:4346740/1‑150 (MQ=255)
                                                    tAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAGGCCga                                                                                                                                       >  1:3162556/1‑112 (MQ=38)
                                                    tAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAGGCCga                                                                                                                                       <  2:3162556/113‑2 (MQ=38)
                                                     aGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTc                                                                                                 >  1:4392551/1‑150 (MQ=255)
                                                      gCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATTAGGCCGGCAAATCCCTACAAGGCACAAGGCCCCGCGCCATGCTTCa                                                                                                >  1:5343606/1‑150 (MQ=255)
                                                       ccgcggttGGGCAGAGGTCTAATCGCGTGGCTGTTGGCGCCTGATCGGGGGGCTCTGTCAGGCGCCGGCCTGTAGCGCCCCCATCGCGGACACGCCCGCGATGCGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCGTCaa                                                                                               <  1:2354550/144‑1 (MQ=255)
                                                       ccGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCCTGCTTCaa                                                                                               >  1:4124314/1‑150 (MQ=255)
                                                       ccGAGCTTGGGCAGAGATCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCAAACAAGGCACAAGGCCCCGCGCCATGCTTCaa                                                                                               >  1:6581409/1‑150 (MQ=255)
                                                           gCTTGGGCAGAGGTCTATTTGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCTGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCTATGTGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATaca                                                                                           <  1:1108295/150‑1 (MQ=255)
                                                            cTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATacac                                                                                          >  2:1172230/1‑150 (MQ=255)
                                                            cTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGCATGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCTCACAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATacac                                                                                          <  2:6283498/150‑1 (MQ=255)
                                                             ttGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATAcaca                                                                                         <  2:3822805/150‑1 (MQ=255)
                                                               gggCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGc                                                                                       >  2:2816272/1‑150 (MQ=255)
                                                                ggCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGcc                                                                                      >  2:932689/1‑150 (MQ=255)
                                                                         atAATAGTGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCttt                                                                             <  1:2668767/149‑1 (MQ=255)
                                                                          tAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTg                                                                            >  2:1000927/1‑150 (MQ=255)
                                                                               gtgTGGCTGTTGGTGCCTGTTTGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGGGGTCACGCCCGCGATGTGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAg                                                                       <  2:3028069/148‑1 (MQ=255)
                                                                                 gTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGa                                                                     <  2:5015379/150‑1 (MQ=255)
                                                                                    cctgacagAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCACATCGCGGTCAAGCACGCGATGAGGCCGGCCCATCCAAACAACGCACAAGGCCCCGCGCCCAGCTTCAATACCCCGCCATACCCATTGCCCAGGACGc                                                                  >  1:3518888/8‑150 (MQ=255)
                                                                                      tGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCg                                                                >  1:5625955/1‑150 (MQ=255)
                                                                                          ggAGCCTGAGAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGcc                                                            <  2:6179261/150‑1 (MQ=255)
                                                                                           gAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGccc                                                           <  2:3930517/150‑1 (MQ=255)
                                                                                               cTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCGTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCAt                                                       <  1:5972400/150‑1 (MQ=255)
                                                                                                tGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACCAGGCCCCGCGCCATGCGTCAATACACAGCCATACCCGTTGCCCAGGACGCCGAGCCCTCATg                                                      >  2:686433/1‑150 (MQ=255)
                                                                                                 gATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGa                                                     <  2:6050422/150‑1 (MQ=255)
                                                                                                   ttggctggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCGGCGTTGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAAc                                                   <  1:6430939/148‑1 (MQ=255)
                                                                                                      gctggctCTGTCTGGCGCAGGCCTTTATCGCCCTCATCGCGGTCAAGCCCGCTATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACCGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGc                                                <  1:1888950/150‑1 (MQ=255)
                                                                                                         ggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAAc                                             <  2:876193/150‑1 (MQ=255)
                                                                                                            tctgtctgGCGCCGGCCTGTATCGGCCTCATCGCGGACAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGAc                                          <  1:1466241/150‑1 (MQ=255)
                                                                                                            tctGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGAc                                          <  1:955755/150‑1 (MQ=255)
                                                                                                            tctGTCAGGCGCAGGCCTGTATCGGCCTCATCGCCGGCAAGCCCGCGATGAGGCCGGCCCAACCAAACAAGGCACAAGGCCCCGCGCCAAGCTTCAATACACACCCAAACCCGTTGCCCAGGACGCCGCGCCCACAAGAACCCCAACGAc                                          >  1:6514561/1‑150 (MQ=255)
                                                                                                                tCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTgc                                      >  1:2817206/1‑150 (MQ=255)
                                                                                                                    gcgcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCCGGACGCCGAGCCCTCATGAACCGCAACGACCTGcgccg                                  >  2:6501797/1‑150 (MQ=255)
                                                                                                                          gCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCGTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                            >  2:465012/1‑150 (MQ=255)
                                                                                                                            cTGTATCGGCCTCATCGCGGTCAAGCCCGCTATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACa                          <  1:377340/150‑1 (MQ=255)
                                                                                                                            cTGTATCGGCCTCATCGCGGACACGCCCGCGTTGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCACAACGACCTGCGCCGCGTCGACa                          <  1:3072898/150‑1 (MQ=255)
                                                                                                                             tGTATCGGCCTCATCGCGTTCAAGCCCGCTATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACAt                         <  1:3610580/150‑1 (MQ=255)
                                                                                                                               tATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGa                       <  1:5088985/150‑1 (MQ=255)
                                                                                                                               tATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGa                       <  1:1610832/150‑1 (MQ=255)
                                                                                                                                    gCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACCTGAACctg                  >  2:1110649/1‑150 (MQ=255)
                                                                                                                                    gCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctg                  >  2:5371729/1‑150 (MQ=255)
                                                                                                                                      cTCATCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgct                >  1:1077456/1‑150 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg               <  1:3267627/150‑1 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg               >  1:1296527/1‑150 (MQ=255)
                                                                                                                                       tCATCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg               >  2:1137001/1‑150 (MQ=255)
                                                                                                                                       tCATCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg               >  1:6610323/1‑150 (MQ=255)
                                                                                                                                         aTCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGt             >  1:5277535/1‑150 (MQ=255)
                                                                                                                                          tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg            >  1:2914710/1‑150 (MQ=255)
                                                                                                                                          tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg            >  1:5786888/1‑150 (MQ=255)
                                                                                                                                           cgcgGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGc           <  2:469907/150‑1 (MQ=255)
                                                                                                                                              ggTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGt        >  1:3926432/1‑150 (MQ=255)
                                                                                                                                              ggTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGt        >  1:1792307/1‑150 (MQ=255)
                                                                                                                                               gTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGtt       >  1:2501627/1‑150 (MQ=255)
                                                                                                                                               gTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGtt       >  2:5127489/1‑150 (MQ=255)
                                                                                                                                                tCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTc      <  1:2794362/150‑1 (MQ=255)
                                                                                                                                                tCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTc      >  2:181283/1‑150 (MQ=255)
                                                                                                                                                tCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTc      >  2:5711584/1‑150 (MQ=255)
                                                                                                                                                  aaGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGa    >  1:3121218/1‑150 (MQ=255)
                                                                                                                                                   aGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGaa   <  2:4017143/150‑1 (MQ=255)
                                                                                                                                                   aGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGaa   <  2:5467451/150‑1 (MQ=255)
                                                                                                                                                   aGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGaa   <  2:3279699/150‑1 (MQ=255)
                                                                                                                                                    gCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAg  >  1:6400098/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAG  >  NC_002947/4175029‑4175326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: