Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,172,151 C→A 9.4% 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,172,1510CA 9.4% 320.1 / 17.7 104intergenic (‑293/‑94)katG/PP_3669catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (38/56);  new base A (2/8);  total (40/64)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAGC  >  NC_002947/4172002‑4172300
                                                                                                                                                     |                                                                                                                                                     
cAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGc                                                                                                                                                       <  2:803734/150‑1 (MQ=255)
cAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGc                                                                                                                                                       <  1:1259743/150‑1 (MQ=255)
      cTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                                 <  2:1143876/150‑1 (MQ=255)
      cTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                                 <  2:3078895/150‑1 (MQ=255)
       tGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                                <  2:4128760/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                              <  2:2183319/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                              <  2:2701627/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                              >  1:2152730/1‑150 (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)
               tGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCgg                                                                                                                                        <  2:4448675/150‑1 (MQ=255)
                    tGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGGTGAGGCCGGCACAt                                                                                                                                   >  1:3316703/1‑150 (MQ=255)
                        gcgGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCAt                                                                                                                               <  2:4029645/150‑1 (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)
                                      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)
                                                                                                        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                          <  1:3110887/150‑1 (MQ=255)
                                                                                                                               tATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGa                        <  2:2230888/150‑1 (MQ=255)
                                                                                                                               tATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCGGCGCCGCGTCGACAAGa                        >  2:278012/1‑150 (MQ=255)
                                                                                                                                     ccTCTTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgc                  <  1:4551903/150‑1 (MQ=255)
                                                                                                                                     ccTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgc                  <  1:1743343/150‑1 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCa                                                                         >  2:4874030/1‑93 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCa                                                                         <  1:4874030/93‑1 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg                >  2:724203/1‑150 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg                >  2:2571582/1‑150 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg                >  2:546135/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:741919/150‑1 (MQ=255)
                                                                                                                                          tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg             >  1:2371585/1‑150 (MQ=255)
                                                                                                                                          tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg             >  2:3459008/1‑150 (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)
                                                                                                                                                  aaGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGa     >  1:375354/1‑150 (MQ=255)
                                                                                                                                                   aGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGaa    <  1:4755163/150‑1 (MQ=255)
                                                                                                                                                    gCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAg   >  1:1198307/1‑150 (MQ=255)
                                                                                                                                                     cccGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAGc  >  1:2891064/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAGC  >  NC_002947/4172002‑4172300

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: