Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,162 A→T 10.3% intergenic (+47/+76) PP_4295 → / ← PP_4296 TetR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,1620AT10.3% 441.4 / 28.8 148intergenic (+47/+76)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (57/71);  new base T (4/11);  total (63/85)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

GCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGAT  >  NC_002947/4891014‑4891309
                                                                                                                                                    |                                                                                                                                                   
gCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcac                                                                                                                                                    <  1:3002295/150‑1 (MQ=255)
 cGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcacc                                                                                                                                                   <  1:4786995/150‑1 (MQ=255)
 cGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcacc                                                                                                                                                   >  2:1893379/1‑150 (MQ=255)
   gTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcaccac                                                                                                                                                 >  1:4674747/1‑150 (MQ=255)
   gTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcaccac                                                                                                                                                 >  1:3737987/1‑150 (MQ=255)
     gACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAt                                                                                                                                               <  1:319013/150‑1 (MQ=255)
         ggTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCt                                                                                                                                           <  2:434770/150‑1 (MQ=255)
            aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                        <  2:4654590/150‑1 (MQ=255)
            aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGACTGAACCATCTCGAAGGGCCATTCGCCGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                        <  2:4019132/150‑1 (MQ=255)
             aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAg                                                                                                                                       <  2:2679239/150‑1 (MQ=255)
             aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAg                                                                                                                                       <  1:4431761/150‑1 (MQ=255)
              ggCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAgg                                                                                                                                      <  2:3341969/150‑1 (MQ=255)
               gCCAAGCTGGCAGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGGATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACAATCTCGAAGGGACATTCGCGGGCAAGCCCGCTCCACCAGCATCACCAAATGCCTCACGGt                                                                                                                                     >  2:2938353/1‑150 (MQ=255)
                ccAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGtt                                                                                                                                    <  1:1558258/150‑1 (MQ=255)
                ccAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGtt                                                                                                                                    <  2:3872950/150‑1 (MQ=255)
                  caGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGCATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCAACTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                  <  2:1601790/149‑1 (MQ=255)
                  aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                  <  2:4320770/150‑1 (MQ=255)
                  aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                  >  2:3929935/1‑150 (MQ=255)
                   aGCTGGCCGACAGTGTTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGATTGGTTGTGAGCCTGAGGCAGCCTGAAGCATCTCGTAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCACATGCGTGAAGGTtgtg                                                                                                                                 <  2:3824035/150‑1 (MQ=255)
                   aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:3328018/150‑1 (MQ=255)
                   aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:4593385/150‑1 (MQ=255)
                   aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:1666174/150‑1 (MQ=255)
                   aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:4657739/150‑1 (MQ=255)
                   aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:1514414/150‑1 (MQ=255)
                    gCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCGGAAGGTTGTgg                                                                                                                                >  1:390561/1‑150 (MQ=255)
                     cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                               <  1:5441252/150‑1 (MQ=255)
                        gCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAATCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTc                                                                                                                            >  2:574924/1‑150 (MQ=255)
                         ccGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                           <  2:1394576/150‑1 (MQ=255)
                          cGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                           >  1:3584271/1‑149 (MQ=255)
                              aGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                      <  1:221975/150‑1 (MQ=255)
                               gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  1:1317116/150‑1 (MQ=255)
                               gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  2:5104069/150‑1 (MQ=255)
                               gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCACTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  2:918417/150‑1 (MQ=255)
                               gTGATTATGACGCGGCGGCCGAGACCATCAGTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  2:4321479/150‑1 (MQ=255)
                                tGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACCTGCCTGAAGGTTGTGGAGTCCCTGTggg                                                                                                                    >  2:5007229/1‑150 (MQ=255)
                                  aTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAg                                                                                                                  <  2:586635/150‑1 (MQ=18)
                                   ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGc                                                                                                                 <  1:1494836/150‑1 (MQ=18)
                                   ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 >  1:4227205/1‑134 (MQ=255)
                                   ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:4227205/134‑1 (MQ=255)
                                   ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                 <  2:2320382/150‑1 (MQ=255)
                                     aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               <  1:3763251/150‑1 (MQ=18)
                                     aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               <  1:2308331/150‑1 (MQ=18)
                                     aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCgg                                                                                                               <  1:1662750/150‑1 (MQ=255)
                                         cgcgGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCtt                                                                                                           <  1:4297999/150‑1 (MQ=17)
                                            ggcggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGcc                                                                                                        <  2:3813785/150‑1 (MQ=255)
                                               ggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCcgc                                                                                                     <  2:133574/150‑1 (MQ=255)
                                                   gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               >  2:4161471/1‑136 (MQ=17)
                                                   gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               <  1:4162016/136‑1 (MQ=17)
                                                       ccATCATTCGTCTGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGc                                                                                             <  2:4157397/150‑1 (MQ=255)
                                                        catcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGcc                                                                                            >  1:4487245/1‑150 (MQ=255)
                                                         atcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATGTCGATGGGCCATTCGCGGGCAAGCCCGGTCCCACAGGATCAACACATGCCTGAAGGTTGTGGAGTCCGTGTGGGAGCGGGCTTGCCCGCGAATGGGGCa                                                                                           >  2:2899392/1‑150 (MQ=39)
                                                               cGTATGGTCTTGAAAGGTGGTGAGCCTGAGGCAGCCTGAACCACCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAg                                                                                     <  2:3538131/150‑1 (MQ=255)
                                                                gtttGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATCACCACCTGCCTGAAGGTTGTGGTGTCCCTGTGGGAGGGGGCTTGCCCGCGAATGGGCCAGAACAgg                                                                                    <  2:5162648/147‑1 (MQ=39)
                                                                   tGGTCTTGAAAGGTTGTGTGCCTGAGGCAGCCTGAACCATCTCTAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCCCTTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGGTTGCCCGCGAATGGGCCAGAACAGGCaa                                                                                 <  2:1965184/150‑1 (MQ=39)
                                                                     gTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATa                                                                               <  2:3148952/150‑1 (MQ=255)
                                                                      tCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATag                                                                              >  2:3679941/1‑150 (MQ=255)
                                                                             aGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                                       <  1:5601818/150‑1 (MQ=255)
                                                                                ttGTGAGCCTGAGGCCGCCTGAACAATCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCCCCACATGCCTGAAGGTTGTGGAGCCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGc                                                                    <  1:2707262/150‑1 (MQ=38)
                                                                                tggtgCGCCTGCGGCAGCCTGAACCATCTCGACGGGCCATGCGCGGGCACGCCCGCGCCCCCCGGATCCCCACCTGCCTGACGGGTGTGGCGTCCCTGTGGGAGCGGGCTGGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGc                                                                    <  2:1135695/148‑1 (MQ=255)
                                                                                             gCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCa                                                       >  2:2105995/1‑150 (MQ=255)
                                                                                                   tGACCCCTCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGcc                                                 <  2:5027479/150‑1 (MQ=255)
                                                                                                          tctcGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCCCCACATGCCTGTAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACCGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcc                                          >  2:3589163/1‑150 (MQ=255)
                                                                                                            tcGACGGGCCATACGCGGGCACGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGCCCCTGTGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcccc                                        <  1:3469601/150‑1 (MQ=39)
                                                                                                            ccGACGGGCCTTGCGCGGGCACGCCCGCCCCCACAGGACCCCCACAGCCCTGACGGTGGTGGCGCCCCGGTGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcccc                                        <  1:5260095/149‑1 (MQ=255)
                                                                                                              gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                      >  1:4127105/1‑150 (MQ=255)
                                                                                                                  ggCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAcctcc                                  >  2:5248050/1‑150 (MQ=255)
                                                                                                                  ggCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCAGTGGGAGCGGGCTTGCCCGCGAATGGGCCCGAACAGGCACTAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCaccacc                                  >  2:5247612/1‑150 (MQ=255)
                                                                                                                    ccATTCGCGGGCAAGCCCGCTCCCACAGGGACGACACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc                                >  2:3785388/1‑150 (MQ=17)
                                                                                                                            gggCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAAAAGAGAATTATGCCGCTACCCCGGCATCTGACCGCAGCCCCACCTCCACGATCGCAc                        >  1:4504149/1‑150 (MQ=39)
                                                                                                                                 aGCCCGCTCCCCCTGGATCACCACAGGCCTGATGGTTGTGGTGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATc                   <  2:4698736/150‑1 (MQ=255)
                                                                                                                                 aGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATc                   <  2:4674747/150‑1 (MQ=255)
                                                                                                                                   cccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGc                 >  1:1833223/1‑150 (MQ=255)
                                                                                                                                   cccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGc                 >  1:1313007/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                >  1:2766609/1‑150 (MQ=18)
                                                                                                                                    ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGTCAATAGAGATTTAGGCCGCTACCCCGGCATCTTCCCGCAGCCCCACATCCTCGATCTCACTGATCGca                >  1:4429635/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                >  1:331340/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                >  1:5556559/1‑150 (MQ=255)
                                                                                                                                      gCTCCCACCGTATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCgagga              <  2:44697/150‑5 (MQ=255)
                                                                                                                                      gCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTCCCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGcaca              >  2:3003929/1‑150 (MQ=255)
                                                                                                                                      gCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGcaca              >  1:4930921/1‑150 (MQ=255)
                                                                                                                                        tCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACt            >  1:4852984/1‑150 (MQ=255)
                                                                                                                                         cccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTg           <  2:4127105/150‑1 (MQ=255)
                                                                                                                                          ccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGc          >  2:3727829/1‑150 (MQ=255)
                                                                                                                                          ccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGc          >  1:951113/1‑150 (MQ=255)
                                                                                                                                           cacaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt         >  1:5317852/1‑150 (MQ=255)
                                                                                                                                           cacaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt         >  1:1217515/1‑150 (MQ=255)
                                                                                                                                            acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  1:1035668/1‑150 (MQ=255)
                                                                                                                                            acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  2:3273153/1‑150 (MQ=255)
                                                                                                                                            acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  2:3272633/1‑150 (MQ=255)
                                                                                                                                            acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  2:3123316/1‑150 (MQ=255)
                                                                                                                                            acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  2:5327437/1‑150 (MQ=255)
                                                                                                                                             ccggATCACCACCTGCCGGACGGTTGGGGCGTCCCTGTGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcg       <  1:5368190/148‑1 (MQ=255)
                                                                                                                                              aGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgt      >  1:2852353/1‑150 (MQ=255)
                                                                                                                                              aGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgt      >  1:3044566/1‑150 (MQ=255)
                                                                                                                                               ggATCACCACTTGCCTGACGGTTGTGGCGGCCCGGGGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtc     <  1:609575/150‑1 (MQ=255)
                                                                                                                                                gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCTCACTGCtcgtcg    >  2:4864851/1‑150 (MQ=255)
                                                                                                                                                gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg    >  1:4307345/1‑150 (MQ=255)
                                                                                                                                                gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg    >  2:5464207/1‑150 (MQ=255)
                                                                                                                                                 aTCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa   >  2:3825092/1‑150 (MQ=255)
                                                                                                                                                  tCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGat  >  1:4335102/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGAT  >  NC_002947/4891014‑4891309

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: