Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,596,973 C→A 11.9% intergenic (‑89/‑114) PP_1395 ← / → PP_1396 AraC family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,596,9730CA11.9% 802.7 / 80.2 268intergenic (‑89/‑114)PP_1395/PP_1396AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (113/123);  new base A (16/16);  total (129/139)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 276 total aligned reads displayed.

GCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTC  >  NC_002947/1596824‑1597122
                                                                                                                                                     |                                                                                                                                                     
gCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGcc                                                                                                                                                       >  1:1094042/1‑150 (MQ=255)
gCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGcc                                                                                                                                                       >  2:501916/1‑150 (MQ=255)
          cTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAgg                                                                                                                                             <  1:16081/150‑1 (MQ=255)
           tCCAGGCGCTGCAGGTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGa                                                                                                                                            <  2:910481/150‑1 (MQ=255)
             ccggCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGCTGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATc                                                                                                                                          <  1:4943328/148‑1 (MQ=255)
             cAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATc                                                                                                                                          >  1:4069215/1‑150 (MQ=255)
             cAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATc                                                                                                                                          >  1:2530908/1‑150 (MQ=255)
             cAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATc                                                                                                                                          >  1:6469036/1‑150 (MQ=255)
             cAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATc                                                                                                                                          <  1:1932028/150‑1 (MQ=255)
              aGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                         <  2:2537783/150‑1 (MQ=255)
              aGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCa                                                                                                                                         <  2:6203983/150‑1 (MQ=255)
               ggCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAg                                                                                                                                        <  1:1965674/150‑1 (MQ=255)
               ggCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAg                                                                                                                                        <  2:2133446/150‑1 (MQ=255)
               ggCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAg                                                                                                                                        <  1:3172707/150‑1 (MQ=255)
               ggCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAACTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAg                                                                                                                                        <  1:6744968/150‑1 (MQ=255)
                gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                       <  1:2384606/150‑1 (MQ=255)
                gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                       <  2:5685311/150‑1 (MQ=255)
                gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                       <  1:6034291/150‑1 (MQ=255)
                gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAgg                                                                                                                                       <  1:4210845/150‑1 (MQ=255)
                gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCAgg                                                                                                                                       <  1:5915023/150‑1 (MQ=255)
                    tGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAAt                                                                                                                                   >  1:4347426/1‑150 (MQ=255)
                        gTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGAc                                                                                                                               >  2:2595184/1‑150 (MQ=255)
                          tCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGAACCTGCAGCGCGGCAGGATCAGGCAATCCACAt                                                                                                                             <  2:6469036/150‑1 (MQ=255)
                           ccAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATg                                                                                                                            >  2:4438254/1‑150 (MQ=255)
                                 tcggttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTaa                                                                                                                      >  2:5334270/1‑150 (MQ=255)
                                          tGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAc                                                                                                             <  2:5475060/150‑1 (MQ=255)
                                             aCGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAcgcg                                                                                                          >  1:6843878/1‑150 (MQ=255)
                                              cGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAcgcgc                                                                                                         <  1:2668374/150‑1 (MQ=255)
                                                       ggtCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGt                                                                                                <  2:2825079/150‑1 (MQ=255)
                                                        gtCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCCATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTg                                                                                               >  2:1901098/1‑150 (MQ=255)
                                                           atgatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAATGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGTGTGAgc                                                                                            <  2:3937224/150‑1 (MQ=255)
                                                           atgatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGCTTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAgc                                                                                            <  2:5899908/150‑1 (MQ=255)
                                                             gatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAgcgc                                                                                          >  1:3357311/1‑150 (MQ=255)
                                                                  cTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAAc                                                                                     <  1:3713694/150‑1 (MQ=255)
                                                                   tttCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAAt                                                                                                                     <  2:2764130/117‑1 (MQ=255)
                                                                     tCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGTTTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCtt                                                                                  <  1:3074045/150‑1 (MQ=255)
                                                                       ctctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTgg                                                                                >  2:3856471/1‑150 (MQ=255)
                                                                       ctctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTgg                                                                                <  2:6087079/150‑1 (MQ=255)
                                                                        tctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGc                                                                               >  1:557088/1‑150 (MQ=255)
                                                                           gctgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGc                                                                            <  2:1751287/150‑1 (MQ=255)
                                                                            ctgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGcc                                                                           >  1:4569179/1‑150 (MQ=255)
                                                                             tgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCt                                                                          <  2:3212816/150‑1 (MQ=255)
                                                                             tgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCGAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGccg                                                                          >  1:178153/1‑149 (MQ=255)
                                                                              gACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTg                                                                         >  2:1629379/1‑150 (MQ=255)
                                                                                ccGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGc                                                                       <  2:4483606/150‑1 (MQ=255)
                                                                                 cGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCt                                                                      <  2:5170521/150‑1 (MQ=255)
                                                                                       agcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccg                                                                >  2:3404226/1‑150 (MQ=255)
                                                                                              gTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCt                                                         <  1:1167863/150‑1 (MQ=255)
                                                                                              gTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGAGGATTGCCTGATCCTGACGAGCTGCAGGATCAGGCAAACGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGATCGCATAACCTTGGCAGCCTGGATTTTCAGTCCGCCt                                                         >  2:4292249/1‑150 (MQ=255)
                                                                                                  gggCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCtt                                                     >  1:1290541/1‑150 (MQ=255)
                                                                                                   ggCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTc                                                    >  2:3681420/1‑150 (MQ=255)
                                                                                                    gCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCa                                                   >  2:2237685/1‑150 (MQ=255)
                                                                                                      aaaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                                 <  2:4183404/150‑1 (MQ=255)
                                                                                                      aaaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                                 <  1:2001287/150‑1 (MQ=255)
                                                                                                      aaaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                                 >  1:3889102/1‑150 (MQ=255)
                                                                                                      aaaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                                 <  2:3806392/150‑1 (MQ=255)
                                                                                                        aGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtccg                                                            <  1:2781428/137‑1 (MQ=255)
                                                                                                        aGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGga                                               <  2:4528796/150‑1 (MQ=255)
                                                                                                            gggCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCATGGAGAtt                                           <  1:6292571/150‑1 (MQ=255)
                                                                                                            gggCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGAtt                                           >  2:6016933/1‑150 (MQ=255)
                                                                                                                gcTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGc                                       <  2:5661382/150‑1 (MQ=255)
                                                                                                                gcTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGc                                       >  1:6636602/1‑150 (MQ=255)
                                                                                                                  tAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCaa                                     <  2:6707608/150‑1 (MQ=255)
                                                                                                                     gccattacTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATAGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCAAACCTTGGCAGCCAGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAAAGa                                  >  1:4929664/8‑150 (MQ=25)
                                                                                                                        aaaGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGAccc                               <  1:3468280/150‑1 (MQ=255)
                                                                                                                         aaGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCt                              <  1:536141/150‑1 (MQ=255)
                                                                                                                         aaGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCt                              <  2:5775048/150‑1 (MQ=255)
                                                                                                                          aGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCtt                             <  2:6851792/150‑1 (MQ=255)
                                                                                                                           gCTGCATGTGGATTGCCTGATCCTGCCGCCCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTa                            >  1:4173124/1‑150 (MQ=255)
                                                                                                                            cTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaa                           <  1:6488226/150‑1 (MQ=255)
                                                                                                                            cTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaa                           <  2:2121177/150‑1 (MQ=255)
                                                                                                                                aTGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaacaa                       >  2:4842927/1‑150 (MQ=255)
                                                                                                                                aTGTGGATTGCCTGATCCGGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCGAACCGTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCCATGACCCTTAaacaa                       >  2:6485746/1‑150 (MQ=255)
                                                                                                                                 tgtgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAAGCGCCGCGCGGGTGAGCGCCTAAACATGGCAGCCTGGCATTTCCGTACGCCTGCTTCAAGGAGATTGTGCAAATACCCTTAaacaac                      >  1:170704/1‑150 (MQ=255)
                                                                                                                                  gtgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAAcc                     >  2:6618407/1‑150 (MQ=255)
                                                                                                                                  gtgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAAcc                     <  2:474472/150‑1 (MQ=255)
                                                                                                                                  gtgGATTGCCTGATACTGCCGCGCTGCAGGATCAGGCAATAGAAATGCAGTAATGGCCAAACGCGCCGCGCGGGTGAGCGCCTAAACTTGGCATCCTGGCTTTTACGTCCGCCTTCTTAAAGGAGATTGTGAAATGACCCATAAACAAcc                     >  1:2205522/1‑150 (MQ=255)
                                                                                                                                    ggATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCgg                   <  1:5664969/150‑1 (MQ=255)
                                                                                                                                           cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTACACAACCGGTcaccca            >  2:4951340/1‑150 (MQ=255)
                                                                                                                                           cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:1252197/1‑150 (MQ=255)
                                                                                                                                           cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:2475785/1‑150 (MQ=255)
                                                                                                                                           cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  1:5829255/1‑150 (MQ=255)
                                                                                                                                           cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  1:5461314/1‑150 (MQ=255)
                                                                                                                                           cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:2933796/1‑150 (MQ=255)
                                                                                                                                           cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccca            >  2:3033288/1‑150 (MQ=255)
                                                                                                                                            tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccac           >  1:782616/1‑150 (MQ=255)
                                                                                                                                            tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccac           >  1:3605305/1‑150 (MQ=255)
                                                                                                                                            tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccac           >  2:3075714/1‑150 (MQ=255)
                                                                                                                                            tGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccac           <  1:81607/150‑1 (MQ=255)
                                                                                                                                             gATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccacc          <  2:3149385/150‑1 (MQ=255)
                                                                                                                                             gATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccacc          >  1:2589171/1‑150 (MQ=255)
                                                                                                                                               tCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaaa                          >  2:1449986/1‑132 (MQ=255)
                                                                                                                                               tCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaaa                          <  1:1449986/132‑1 (MQ=255)
                                                                                                                                               tCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCtt        <  1:4861056/150‑1 (MQ=255)
                                                                                                                                                ccTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTg                                                                                 <  1:958008/76‑1 (MQ=255)
                                                                                                                                                ccTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGATCACCCACCTTg       >  2:3263285/1‑150 (MQ=255)
                                                                                                                                                  tGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCt     >  1:3084508/1‑150 (MQ=255)
                                                                                                                                                    ccGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCttt   >  2:5253852/1‑150 (MQ=255)
                                                                                                                                                     cgcgCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTc  >  1:4407049/1‑150 (MQ=255)
                                                                                                                                                     cgcgCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTc  >  2:2297583/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCGGCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTTTC  >  NC_002947/1596824‑1597122

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: