Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,172,146 T→G 10.7% intergenic (‑288/‑99) 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,172,1460TG10.7% 252.6 / 21.9 105intergenic (‑288/‑99)katG/PP_3669catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (36/55);  new base G (3/8);  total (39/66)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 105 total aligned reads displayed.

ACGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTT  >  NC_002947/4171999‑4172294
                                                                                                                                                   |                                                                                                                                                    
aCGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCa                                                                                                                                                    <  1:617953/150‑1 (MQ=255)
aCGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCa                                                                                                                                                    >  2:2956636/1‑150 (MQ=255)
   cAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGc                                                                                                                                                 <  2:803734/150‑1 (MQ=255)
   cAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGc                                                                                                                                                 <  1:1259743/150‑1 (MQ=255)
         cTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                           <  2:3078895/150‑1 (MQ=255)
         cTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                           <  2:1143876/150‑1 (MQ=255)
          tGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                          <  2:4128760/150‑1 (MQ=255)
            aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                        >  1:2152730/1‑150 (MQ=255)
            aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                        <  2:2183319/150‑1 (MQ=255)
            aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                        <  1:1433229/150‑1 (MQ=255)
            aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                        <  2:3595764/150‑1 (MQ=255)
            aCGCGCGGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                        <  1:1948490/150‑1 (MQ=255)
             cgcgCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                       <  1:1262813/150‑1 (MQ=255)
             cgcgCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                       <  2:3104829/150‑1 (MQ=255)
              gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAgg                                                                                                                                      >  2:4755163/1‑150 (MQ=255)
               cgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAGGc                                                                                                                                     <  2:3767249/150‑1 (MQ=255)
                gcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGcc                                                                                                                                    >  1:101553/1‑150 (MQ=255)
                 cTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGTGGTCTAATAGAGTGGCTGTTGGAGCCTGATTGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCg                                                                                                                                   <  1:610352/150‑1 (MQ=255)
                       tGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGGTGAGGCCGGCACAt                                                                                                                             >  1:3316703/1‑150 (MQ=255)
                           gcgGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCAt                                                                                                                         <  2:4029645/150‑1 (MQ=255)
                           gcgGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCAACGCGGTCAAGCCCGCCATGAGGCCGGCACATCCAt                                                                                                                         >  1:2445863/1‑150 (MQ=255)
                                  tCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCTTAGGTCTAATAGTGTGGCTGTTGGAGCCTGTTAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggc                                                                                                                  <  2:4748830/150‑1 (MQ=255)
                                  tCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGTGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCTGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCTATGAGGCCGGCACATCCATacaaggc                                                                                                                  <  1:1268479/150‑1 (MQ=255)
                                  tCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATCGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggc                                                                                                                  <  2:4748798/150‑1 (MQ=255)
                                  tCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggc                                                                                                                  >  2:4942656/1‑150 (MQ=255)
                                      ggggCCAGGAGTAAGTCTTGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGTTCGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTAATCGCGGTCAAGCCGGGGATGTGGCCGGCACATCCATacaaggcacaa                                                                                                              <  2:2672197/150‑1 (MQ=255)
                                      ggggCCACGAGAAAGACTAGCCTTGATTGGGCAGAGGTCTAATTGAGGGGGTGTTGGAGCCTGATAGGCTGGATCTGTCTGGCGCTGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggcacaa                                                                                                              <  1:4868757/150‑1 (MQ=255)
                                        ggCCACGAGAAAGACTAGCGGAGCTGGGGCAGTGGTATAATAGTGTGGCTGTTGGAGCCTGATTGGCTGGCTCTGTCAGGCGCTGGCCTGTATCGGCCTCATCGCGGACAAGCCCGCGATGAGGCCGGCACATCCATacaaggcacaagg                                                                                                            <  2:461480/150‑1 (MQ=255)
                                         gCCACGAGAAAGACTAGCCGTGGTTGGGCTGTGTTATAATTGCGTGGGTGTTGGCGCCTGATTGGGTGGCTCTGTCAGGCGCCGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGTGGCCGGCACATCCATacaaggcacaaggc                                                                                                           <  1:318162/150‑1 (MQ=255)
                                         gCCACGACAACGCCTAGCCGCGCTGGGGCCGAGGTCTAATAGAGGGGCTGTTGGAGCCTGATAGGCTGGCTCTTACAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCACGCCCGCGATGAGGCCGGCACATCCATacaaggcacaaggc                                                                                                           <  1:5246956/150‑1 (MQ=255)
                                           cACGATAATGACTAGCCGAGCTTGGGCAGAGGTATAATAGTGTGGCTGTTGGAGCCTGATAGGGTGGCTCTTTCAGGCGCAGGCCTTTATGGGCCTCATCGCGGTAACGCCCGCGTTGAGGCCGGCACATCCATACAAGGCACAAGGccc                                                                                                         <  1:3840402/150‑1 (MQ=255)
                                           cACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                        >  1:2483886/1‑119 (MQ=39)
                                           cACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                        <  2:2483886/119‑1 (MQ=39)
                                            aCGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGcccc                                                                                                        >  1:4031394/1‑150 (MQ=255)
                                                gaACGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGGGGGTGTTGGAGCCTGATCGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCgcgc                                                                                                    <  2:3785351/150‑1 (MQ=255)
                                                gaAAGACTCGCCGAGCTTGGGCTGAGGTCTAATAGCGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTAGCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCgcgc                                                                                                    <  2:3785373/150‑1 (MQ=255)
                                                gaAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCgcgc                                                                                                    >  2:639564/1‑150 (MQ=255)
                                                  acgACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATTAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCa                                                                                                  <  2:2975471/148‑1 (MQ=255)
                                                  aaGACTAGCCGTGCTTGGGCAGAGGTATAATAGTGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCa                                                                                                  <  2:1604672/150‑1 (MQ=255)
                                                    gACTATCCGAGCATGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCACATCGAGGTAAATCCCGCGATGAAGAAGGCACATCCATACAAGGCACAAGGCCCCGCGCCATg                                                                                                >  2:1594690/1‑150 (MQ=255)
                                                     aCTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGc                                                                                               <  1:4988055/150‑1 (MQ=255)
                                                       tAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtc                                                                                             >  1:5205995/1‑149 (MQ=255)
                                                         gCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCa                                                                                           >  1:2105828/1‑150 (MQ=255)
                                                            gtgCTTGGGCAGAGGTATATTAGAGTGGCTGTTGGAGCCTGATTGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCAGCGCGGTCAAGCCCGCGTTGTGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATa                                                                                        <  1:743013/148‑1 (MQ=255)
                                                            gAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtcggcta                                                                                        <  2:5205995/150‑7 (MQ=255)
                                                              gCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATTAGGCCGGCACATCCAAACAAGGCACAATGCCCCGCGCCATGCTTCAATaca                                                                                      >  1:4489126/1‑150 (MQ=255)
                                                              gCTGGGGCAGAGCTCTAATCGCGTGGCTGTTGGAGCCTGATTGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGCCACGCCCGCCATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATaca                                                                                      <  1:2838581/150‑1 (MQ=255)
                                                                   ggCAGTGGACTAATCGCGGGGCTGTTGGAGCGTGTTTGGCTGGCTCTGTCCGGCGCTGGCCTGTATCGGCCTCATCGGGGTCACGCCCGCTATTAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGcc                                                                                 <  1:48478/150‑1 (MQ=255)
                                                                       gagGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGTTGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATAc                                                                             <  2:674545/150‑1 (MQ=255)
                                                                        agGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATAcc                                                                            <  2:3035254/150‑1 (MQ=255)
                                                                             tAATAGAGTGGCTGTTGGAGCCTGTTAGGGTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCTGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTg                                                                       <  2:1533376/150‑1 (MQ=255)
                                                                             tAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTg                                                                       <  2:3816460/150‑1 (MQ=255)
                                                                             tAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCAAACCGGTCAAGCCCCCGAAGATGCCCGCACCAACAAACAATGCACAATGCCCCGCGCCAAGCTTAAAAACACACCCAAACCCattt                                                                       >  1:3810689/1‑149 (MQ=255)
                                                                               aTAGTGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGcc                                                                     <  2:1892123/150‑1 (MQ=255)
                                                                               aTAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATTCAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGcc                                                                     >  1:4306322/1‑150 (MQ=255)
                                                                                tAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGccc                                                                    >  2:530991/1‑150 (MQ=255)
                                                                                 agagTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCa                                                                   >  2:4548840/1‑150 (MQ=255)
                                                                                   cgTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCACGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAgg                                                                 <  2:1064643/149‑1 (MQ=255)
                                                                                       gCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGc                                                             >  2:3471332/1‑150 (MQ=255)
                                                                                          gTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGa                                                          <  1:3471332/150‑1 (MQ=255)
                                                                                             ggAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGcc                                                       <  1:4942656/150‑1 (MQ=255)
                                                                                              gAGCCTGCTCGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCACGCCCGCGATGCGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGccc                                                      <  2:1308422/150‑1 (MQ=255)
                                                                                               aGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCt                                                     >  2:4454276/1‑150 (MQ=255)
                                                                                                gcgttaTAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTATTCGCGGACAACCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTc                                                    <  1:2486743/145‑1 (MQ=255)
                                                                                                gCCTGATCGGCGGGCTCGGTCAGGCGCCGGCCTGTATCGGCCTCATCGCGGTCACGCCCGCGTTGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTc                                                    <  1:1110495/150‑1 (MQ=255)
                                                                                                 ccTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCa                                                   <  2:5122647/150‑1 (MQ=255)
                                                                                                    gATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCACATGa                                                >  1:1250735/1‑150 (MQ=255)
                                                                                                          ctggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCa                                          <  1:4272192/150‑1 (MQ=255)
                                                                                                           tggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAATCCCGCGATGAGGGCGGCACATCCATACAAGGCAGAAGGCCCCGCGCCTTGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGTGCCCTCATGAACCGCaa                                         >  2:2844838/1‑150 (MQ=255)
                                                                                                            ggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAAc                                        >  2:2789878/1‑150 (MQ=255)
                                                                                                            ggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAAc                                        <  2:375862/150‑1 (MQ=255)
                                                                                                                  gTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCTTTGTGGCCGGCACATCCATACAAGGCACATGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTg                                  <  1:2956636/150‑1 (MQ=255)
                                                                                                                   tCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTTCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTgc                                 >  1:2575196/1‑150 (MQ=255)
                                                                                                                       gcgcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGcgccg                             >  1:1771810/1‑150 (MQ=255)
                                                                                                                        cgcTGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGTGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGcgccgc                            <  1:4279486/150‑1 (MQ=255)
                                                                                                                         gcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCAGGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCcgcg                           >  2:4172307/1‑150 (MQ=255)
                                                                                                                              ccTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGAc                      <  2:4480947/150‑1 (MQ=255)
                                                                                                                               cTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACa                     <  2:1776579/150‑1 (MQ=255)
                                                                                                                               cTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACa                     <  2:4517396/150‑1 (MQ=255)
                                                                                                                                tGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACAt                    <  2:2418371/150‑1 (MQ=255)
                                                                                                                                  tATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGa                  <  2:2230888/150‑1 (MQ=255)
                                                                                                                                  tATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCGGCGCCGCGTCGACAAGa                  >  2:278012/1‑150 (MQ=255)
                                                                                                                                      ggCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACct              >  2:3613001/1‑150 (MQ=255)
                                                                                                                                        ccTCTTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgc            <  1:4551903/150‑1 (MQ=255)
                                                                                                                                        ccTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgc            <  1:1743343/150‑1 (MQ=255)
                                                                                                                                          tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCa                                                                   <  1:4874030/93‑1 (MQ=255)
                                                                                                                                          tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCa                                                                   >  2:4874030/1‑93 (MQ=255)
                                                                                                                                          tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg          >  2:2571582/1‑150 (MQ=255)
                                                                                                                                          tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg          >  2:546135/1‑150 (MQ=255)
                                                                                                                                          tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg          >  2:724203/1‑150 (MQ=255)
                                                                                                                                          tCATCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg          >  2:1625619/1‑150 (MQ=255)
                                                                                                                                          tCATCGCGGGCTTGACCGCGATGAGGCCGATACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg          <  1:925070/150‑1 (MQ=255)
                                                                                                                                            aTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGt        <  2:2055852/150‑1 (MQ=255)
                                                                                                                                             tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg       >  2:3459008/1‑150 (MQ=255)
                                                                                                                                             tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg       >  1:2371585/1‑150 (MQ=255)
                                                                                                                                             tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg       <  2:741919/150‑1 (MQ=255)
                                                                                                                                             tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg       >  1:4685412/1‑150 (MQ=255)
                                                                                                                                              cgcgGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGc      >  1:2306239/1‑150 (MQ=255)
                                                                                                                                               gcgGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACAt                    <  1:4550210/135‑1 (MQ=255)
                                                                                                                                               gcgGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTGTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACAt                    >  2:4550210/1‑135 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTT  >  NC_002947/4171999‑4172294

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: