Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,361,360 G→C 5.4% G94A (GGC→GCC)  PP_1181 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,361,3600GC 5.4% 1803.4 / 63.8 534G94A (GGC→GCC) PP_1181two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (266/237);  new base C (7/22);  total (273/261)
Fisher's exact test for biased strand distribution p-value = 3.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 543 total aligned reads displayed.

CATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAACCGCC  >  NC_002947/1361211‑1361509
                                                                                                                                                     |                                                                                                                                                     
   catcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCt                                                                                                                                                    <  2:8394157/150‑1 (MQ=255)
   catcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCt                                                                                                                                                    <  1:879958/150‑1 (MQ=255)
     tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCCGTACCAGGGGCCGATCCTCATGCTCCCCGCGCGCAGTGATGAACTGGACCAGGGCCAGGGCCTgg                                                                                                                                                  >  2:7459483/1‑150 (MQ=255)
     tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGACCTGGACCAGGTCCAGGGCCTgg                                                                                                                                                  >  2:5647882/1‑150 (MQ=255)
     tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTgg                                                                                                                                                  >  2:5868858/1‑150 (MQ=255)
        aCAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGAcc                                                                                                                                               >  2:3085639/1‑150 (MQ=255)
        aCAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGAcc                                                                                                                                               >  2:221537/1‑150 (MQ=255)
          agccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCTg                                                                                                                                             <  2:6300302/150‑1 (MQ=255)
            ccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTggg                                                                                                                                           <  1:1141510/150‑1 (MQ=255)
             cagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgc                                                                                                                                          >  1:7715251/1‑150 (MQ=255)
               gccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgcgc                                                                                                                                        >  2:6088729/1‑150 (MQ=255)
               gccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgcgc                                                                                                                                        >  1:1762452/1‑150 (MQ=255)
               gccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgcgc                                                                                                                                        >  2:2402836/1‑150 (MQ=255)
                   gACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcgac                                                                                                                                    <  1:6514896/150‑1 (MQ=255)
                     ccTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCTGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcgacga                                                                                                                                  <  2:3376335/150‑1 (MQ=255)
                     ccTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcgacga                                                                                                                                  <  2:1278979/150‑1 (MQ=255)
                        ggTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCTGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGGGATGAACTGGACCAGGTCCCGGGCCTGGACCTGGGCGCCGACGACTa                                                                                                                               <  1:3717023/150‑1 (MQ=255)
                        ggTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCTATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTa                                                                                                                               <  1:8418773/150‑1 (MQ=255)
                           gATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGt                                                                                                                            <  1:6773456/150‑1 (MQ=255)
                               cTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGc                                                                                                                        >  1:2132463/1‑150 (MQ=255)
                               atCGACCTGGTGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGc                                                                                                                        >  2:4526566/2‑150 (MQ=255)
                                  gACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCaaa                                                                                                                     >  1:6639223/1‑150 (MQ=255)
                                        aTGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTg                                                                                                               <  2:5492198/150‑1 (MQ=255)
                                            tCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCGGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGcc                                                                                                           >  1:318853/1‑150 (MQ=255)
                                            tCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGcc                                                                                                           >  2:3280205/1‑150 (MQ=255)
                                             ccccGGCGAGGACGGCCTGAGCATCTGCCGCGGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGccc                                                                                                          <  1:2429181/150‑1 (MQ=255)
                                                    gAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtg                                                                                                   <  2:3170096/150‑1 (MQ=255)
                                                      ggACGGCCTGAGCATCTGCCGCCGGGTGCGCAGACAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGAACAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTTCGCCCACGCTtgct                                                                                                 >  1:987717/1‑150 (MQ=255)
                                                       gACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtgctg                                                                                                >  2:5220338/1‑150 (MQ=255)
                                                         cGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTgctgct                                                                                              <  2:5015838/150‑1 (MQ=255)
                                                         cGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTgctgct                                                                                              <  1:8553048/150‑1 (MQ=255)
                                                          ggCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGctgctg                                                                                             >  2:8493674/1‑150 (MQ=255)
                                                            ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGc                                                                                           <  1:2024601/150‑1 (MQ=255)
                                                             cTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGGTGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCa                                                                                          >  1:1486075/1‑150 (MQ=255)
                                                              tGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCAc                                                                                         >  2:692477/1‑150 (MQ=255)
                                                              tGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCAc                                                                                         <  2:171438/150‑1 (MQ=255)
                                                                  cATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTAt                                                                                     >  1:156375/1‑150 (MQ=255)
                                                                  cATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTAt                                                                                     >  2:823371/1‑150 (MQ=255)
                                                                     cTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATccc                                                                                  >  2:7714805/1‑149 (MQ=255)
                                                                     cTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCa                                                                                  >  1:1903982/1‑150 (MQ=255)
                                                                      tGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGTTGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAg                                                                                 <  2:1486075/150‑1 (MQ=255)
                                                                        ccgcggGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCTGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGc                                                                               <  1:4767755/150‑1 (MQ=255)
                                                                        ccgccgGGTGCGCAGCCAGTCCCTGGGGCCGATCCTCATTGTCACGGCGCGCAGTGATTATATGGACCAGGTCCAGGGCCTGGACGTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCCCGCTTGCTGCTGGCACGTATCCAGGc                                                                               <  1:2756322/150‑1 (MQ=255)
                                                                         cgccgGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGcc                                                                              <  1:4185593/150‑1 (MQ=255)
                                                                            cgGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCctg                                                                           <  2:2548722/150‑1 (MQ=255)
                                                                             ggtggcgcTGCCTGTACCAGGGGCATATCCTCATGCTCCCCGCGCGCTGTTATGAAGTGGCCCTGGGCCAGGGCGTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCctgc                                                                          <  1:2352926/146‑1 (MQ=255)
                                                                                 gcgcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTgcg                                                                      <  1:7289798/150‑1 (MQ=255)
                                                                                  cgcAGCCAGTACCCGGGGCCGATCCTCAGGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACATCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                                     <  1:7872077/150‑1 (MQ=255)
                                                                                  cgcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                                     <  2:2877666/150‑1 (MQ=255)
                                                                                  cgcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                                     <  1:1279817/150‑1 (MQ=255)
                                                                                    cAGCCAGTACCTGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGcgccg                                                                   <  1:375062/150‑1 (MQ=255)
                                                                                    cAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTCTCCAGGCCCTGCTGcgccg                                                                   >  2:8641683/1‑150 (MQ=255)
                                                                                    cAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGcgccg                                                                   >  2:5751746/1‑150 (MQ=255)
                                                                                     aGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGcgccgc                                                                  >  2:4030388/1‑150 (MQ=255)
                                                                                       ccAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                                >  1:3338497/1‑150 (MQ=255)
                                                                                       ccAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                                <  1:5871626/150‑1 (MQ=255)
                                                                                       ccAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                                >  1:8333367/1‑150 (MQ=255)
                                                                                        cAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCTGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGc                                                               <  1:7528328/150‑1 (MQ=255)
                                                                                          gTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGa                                                             >  1:5403295/1‑150 (MQ=255)
                                                                                            aCCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGaaa                                                           >  1:511933/1‑150 (MQ=255)
                                                                                             ccAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAAc                                                          <  1:611579/150‑1 (MQ=255)
                                                                                              cAGGGGCCTATCCTCATGCTCACCGCGCGCAGTGATGAACGGGCCCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAAcc                                                         <  1:2963322/150‑1 (MQ=255)
                                                                                                ggggCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                                       <  2:4754456/150‑1 (MQ=255)
                                                                                                 gggCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTg                                                      <  2:6820314/150‑1 (MQ=255)
                                                                                                  ggCCGATCCTCATGCTCACCGCGCGCAGTGATGACCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTgg                                                     <  1:2610094/150‑1 (MQ=255)
                                                                                                  ggCCGATCCTCATGCTCACCGCGCGCAGGGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTgg                                                     <  1:7996419/150‑1 (MQ=255)
                                                                                                   gCCGACCCTCAGGCTCCCCGCGCGCAGGGATGACCGGGCCCCGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGa                                                    <  1:4428998/150‑1 (MQ=255)
                                                                                                     cGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACa                                                  <  1:494861/150‑1 (MQ=255)
                                                                                                     cGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACa                                                  >  2:2056321/1‑150 (MQ=255)
                                                                                                     cGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGGCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACa                                                  <  1:6939243/150‑1 (MQ=255)
                                                                                                      gATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAg                                                 <  2:4335210/150‑1 (MQ=255)
                                                                                                            cATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGcgcc                                                                    >  1:4530726/1‑125 (MQ=255)
                                                                                                            cATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgcg                                           >  1:5766687/1‑150 (MQ=255)
                                                                                                               gCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCa                                        >  2:3857227/1‑150 (MQ=255)
                                                                                                               gCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCa                                        >  2:7008633/1‑150 (MQ=255)
                                                                                                                   aCCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGAc                                    >  2:8018203/1‑150 (MQ=255)
                                                                                                                      gcgcgcTGTGATGAACTGGTCCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTg                                 <  1:8115591/150‑1 (MQ=255)
                                                                                                                        gcgcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGc                               <  2:4573832/150‑1 (MQ=255)
                                                                                                                          gcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCt                             >  1:1349241/1‑150 (MQ=255)
                                                                                                                              tgctgaCCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGGCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAACCCGTGGACAGCAAGCGCCAGGACCTGGCCTTCgg                         <  1:6253201/147‑1 (MQ=255)
                                                                                                                                 tgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTgc                      <  2:4230387/150‑1 (MQ=255)
                                                                                                                                  gaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTgcg                     >  1:5624471/1‑150 (MQ=255)
                                                                                                                                   aaCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTgcgc                    >  2:6489293/1‑150 (MQ=255)
                                                                                                                                      tGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGc                 >  2:1836911/1‑150 (MQ=255)
                                                                                                                                         aCCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTa              <  2:2733078/150‑1 (MQ=255)
                                                                                                                                         aCCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTa              >  1:6300830/1‑150 (MQ=255)
                                                                                                                                          ccAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTAt             >  1:5287980/1‑150 (MQ=255)
                                                                                                                                          ccAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTAt             >  2:7575348/1‑150 (MQ=255)
                                                                                                                                          ccAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTAt             >  2:3691661/1‑150 (MQ=255)
                                                                                                                                          ccAGGTCCAGGCACTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTAt             >  2:1620814/1‑150 (MQ=255)
                                                                                                                                           cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGtctc            >  2:7186401/1‑150 (MQ=255)
                                                                                                                                           cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc            >  1:1010909/1‑150 (MQ=255)
                                                                                                                                           cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc            >  1:7359091/1‑150 (MQ=255)
                                                                                                                                            aGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCg           >  2:1370309/1‑150 (MQ=255)
                                                                                                                                             ggTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGa          >  1:4661761/1‑150 (MQ=255)
                                                                                                                                              gTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGAc         >  1:4034592/1‑150 (MQ=255)
                                                                                                                                              gTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGATACGCAGCGCACCGaa         <  1:3506676/150‑2 (MQ=255)
                                                                                                                                              gTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCGTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAACCCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGAc         >  2:6653028/1‑150 (MQ=255)
                                                                                                                                                 cAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAAc      >  1:3112860/1‑150 (MQ=255)
                                                                                                                                                    ggCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAACCGc   <  1:387955/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAACCGCC  >  NC_002947/1361211‑1361509

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: