Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,172,149 A→T 9.9% intergenic (‑291/‑96) 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,1490AT 9.9% 248.7 / 18.9 103intergenic (‑291/‑96)katG/PP_3669catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (36/50);  new base T (2/8);  total (38/65)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

Only 100 of 105 total aligned reads displayed.

CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAA  >  NC_002947/4172002‑4172298
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                                                            <  2:2701627/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                            <  2:2183319/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                            >  1:2152730/1‑150 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                            <  2:3595764/150‑1 (MQ=255)
         aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                            <  1:1433229/150‑1 (MQ=255)
         aCGCGCGGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                            <  1:1948490/150‑1 (MQ=255)
          cgcgCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                           <  2:3104829/150‑1 (MQ=255)
          cgcgCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                           <  1:1262813/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)
                                   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                                                                                                                                            <  2:2483886/119‑1 (MQ=39)
                                        cACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                            >  1:2483886/1‑119 (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)
                                                      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:724203/1‑150 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg              >  2:546135/1‑150 (MQ=255)
                                                                                                                                       tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg              >  2:2571582/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:4685412/1‑150 (MQ=255)
                                                                                                                                          tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg           >  2:3459008/1‑150 (MQ=255)
                                                                                                                                          tCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg           >  1:2371585/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)
                                                                                                                                               gTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGtt      <  1:546135/150‑1 (MQ=255)
                                                                                                                                                  aaGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGa   >  1:375354/1‑150 (MQ=255)
                                                                                                                                                   aGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGaa  <  1:4755163/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAA  >  NC_002947/4172002‑4172298

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: