Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,570 T→C 9.8% Y269Y (TAT→TAC rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947542,5700TC 9.8% 536.6 / 33.8 184Y269Y (TAT→TACrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (82/84);  new base C (7/11);  total (89/95)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 186 total aligned reads displayed.

CTGAAGCTGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTC  >  NC_002947/542424‑542719
                                                                                                                                                  |                                                                                                                                                     
cTGAAGCTGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGc                                                                                                                                                    >  2:3737737/1‑150 (MQ=255)
   aaGCTGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCgg                                                                                                                                                 <  1:316193/150‑1 (MQ=255)
   aaGCTGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCgg                                                                                                                                                 <  1:4980502/150‑1 (MQ=255)
     gCTGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGt                                                                                                                                               <  2:2583059/150‑1 (MQ=255)
      cTGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTg                                                                                                                                              >  2:3440322/1‑150 (MQ=255)
       tGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGa                                                                                                                                             <  2:814809/150‑1 (MQ=255)
           ggAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCaa                                                                                                                                         >  1:5199183/1‑150 (MQ=255)
           ggAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCaa                                                                                                                                         >  1:2666285/1‑150 (MQ=255)
               gCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGt                                                                                                                                     >  2:3476465/1‑150 (MQ=255)
               gCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGt                                                                                                                                     <  2:550761/150‑1 (MQ=255)
               gCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGt                                                                                                                                     >  1:3578436/1‑150 (MQ=255)
                cTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTa                                                                                                                                    >  2:2438690/1‑150 (MQ=255)
                 tttCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTaa                                                                                                                                   <  2:3819388/150‑1 (MQ=255)
                 tttCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTaa                                                                                                                                   >  1:4706621/1‑150 (MQ=255)
                     cAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTaacaac                                                                                                                               <  1:1954760/150‑1 (MQ=255)
                       ggaatccggCAACCGTCCTGAGTGGATGATTATGAAGATCGTTCCGGTTATTCCACCGGATCTGCGTCCTCTGGTTCCGCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCg                                                                                                                             <  1:3409270/143‑1 (MQ=255)
                        ggCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGt                                                                                                                            <  1:1162331/150‑1 (MQ=255)
                        ggCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGt                                                                                                                            <  2:3555919/150‑1 (MQ=255)
                         gCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTc                                                                                                                           >  1:3354394/1‑150 (MQ=255)
                          cTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCt                                                                                                                          >  1:2791176/1‑150 (MQ=255)
                          cTCGGGCAACCTGCCTGAGAGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGACCGCTGGTACCGCTGGATGGTGGCCGTTTCGCGACTACCGACCTCAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCt                                                                                                                          >  2:3124197/1‑150 (MQ=255)
                             gggCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCgg                                                                                                                                                 <  2:4160237/124‑1 (MQ=255)
                             gggCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGaa                                                                                                                       <  1:3584380/150‑1 (MQ=255)
                             cggCAACCGTCCTGAGTGGATGATTATGAAGATCGTTCCGGTTATTCCACCGGATCTGAGTCCTCTGGTTCCGCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGaa                                                                                                                       <  2:5299384/149‑1 (MQ=255)
                              ggCAACCGTCCTGAGTGGATGATTATGAAGATCGTTCCGGTTATTCCACCGGATCTGCGTCCTCTGGTTCCGCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAg                                                                                                                      >  2:2971305/1‑150 (MQ=255)
                                   cctgcctgAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCct                                                                                                                 <  2:3277643/150‑1 (MQ=255)
                                   cctgcctgAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCATAACAACCGTCTGAAGCGCct                                                                                                                 <  2:4405117/150‑1 (MQ=255)
                                      gcctgAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCctgct                                                                                                              <  2:2659101/150‑1 (MQ=255)
                                      gcctgAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCctgct                                                                                                              <  1:4599299/150‑1 (MQ=255)
                                         tgAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGa                                                                                                           >  1:1042906/1‑150 (MQ=255)
                                          gAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGac                                                                                                          <  1:2204199/150‑2 (MQ=255)
                                          gAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGAt                                                                                                          <  1:2339089/150‑1 (MQ=255)
                                               gATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTc                                                                                                     <  1:3440322/150‑1 (MQ=255)
                                                 tGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCgg                                                                                                   >  2:675642/1‑150 (MQ=255)
                                                    tCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGgcgc                                                                                                >  2:1845168/1‑150 (MQ=255)
                                                      cTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCg                                                                                              >  1:5389665/1‑150 (MQ=255)
                                                        gACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGa                                                                                            >  1:253844/1‑150 (MQ=255)
                                                         aCCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACACCCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGAc                                                                                           <  2:4010344/150‑1 (MQ=255)
                                                            gTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCt                                                                                                                          <  2:5217645/116‑1 (MQ=255)
                                                            gTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCt                                                                                                                          >  1:5217645/1‑116 (MQ=255)
                                                             tCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCgg                                                                                                                                                 <  1:2068898/92‑1 (MQ=255)
                                                               cTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACatcatc                                                                                     <  1:4714777/150‑1 (MQ=255)
                                                               cTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACatcatc                                                                                     <  2:364128/150‑1 (MQ=255)
                                                               cTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACatcatc                                                                                     <  2:2471730/150‑1 (MQ=255)
                                                                tGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCg                                                                                    >  1:441981/1‑150 (MQ=255)
                                                                   cAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTgc                                                                                 <  1:788684/150‑1 (MQ=255)
                                                                   cAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTgc                                                                                 <  2:1288368/150‑1 (MQ=255)
                                                                              ccgGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGt                                                                                                                            <  1:5068443/96‑1 (MQ=255)
                                                                              ccgGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGt                                                                                                                            >  2:5068443/1‑96 (MQ=255)
                                                                               cgGATCTGCGTCCTCTGGTTCCGCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAgc                                                                     <  2:3362016/150‑1 (MQ=255)
                                                                               cgGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAgc                                                                     <  1:5020508/150‑1 (MQ=255)
                                                                                    cTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGTAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCAtg                                                                >  2:4493748/1‑150 (MQ=255)
                                                                                      gCGTCCTCTGGTTCCGCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTAtgtt                                                              >  2:2840544/1‑148 (MQ=255)
                                                                                         tCCTCTGGTTCCGCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTATGTTGCa                                                           <  1:767917/150‑1 (MQ=255)
                                                                                         tCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCa                                                           >  2:2945287/1‑150 (MQ=255)
                                                                                          ccGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAg                                                          >  1:480776/1‑150 (MQ=255)
                                                                                            gCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGa                                                        >  1:1301905/1‑150 (MQ=255)
                                                                                             cTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGaa                                                       <  1:1578263/150‑1 (MQ=255)
                                                                                             cTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGaa                                                       >  1:1970458/1‑150 (MQ=255)
                                                                                             cTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGaa                                                       <  1:4075539/150‑1 (MQ=255)
                                                                                             cTGGTTCCGCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTATGTTGCAGGaa                                                       <  2:536004/150‑1 (MQ=255)
                                                                                              tGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAg                                                      >  1:2988441/1‑150 (MQ=255)
                                                                                               ggTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGc                                                     >  1:1730363/1‑150 (MQ=255)
                                                                                                gTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCg                                                    >  1:656099/1‑150 (MQ=255)
                                                                                                gTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCg                                                    >  1:1119178/1‑150 (MQ=255)
                                                                                                     gCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGa                                               >  1:2436150/1‑150 (MQ=255)
                                                                                                     gCTCGATGGCGGTCGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTATGTTGCAGGAAGCTGTTGa                                               >  1:3976556/1‑150 (MQ=255)
                                                                                                         gATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGct                                           <  2:248358/150‑1 (MQ=255)
                                                                                                         gATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGct                                           >  2:3342710/1‑150 (MQ=255)
                                                                                                           tggtggCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGTGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGctct                                         <  2:1612234/150‑1 (MQ=255)
                                                                                                              cggtcGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTATGTTGCAGGAAGCTGTTGACGCGctgtt                                      <  2:470/146‑3 (MQ=255)
                                                                                                               ggtcGTTTTGCTACATCCGATCTCAACGACCTTTACCGTCGTGTGATCAACCGTAACAACCGTCTGAAGCGTCTTATTGAACTGCGCGCTCCTGGCATCATTATCCGTAACGAAAAGCGTATGTTGCAGGAAGCTGTTGACGCGctgttc                                     <  2:150863/147‑4 (MQ=255)
                                                                                                               ggCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTctgctg                                     <  1:5252378/150‑1 (MQ=255)
                                                                                                               ggCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTctgctg                                     <  1:443535/150‑1 (MQ=255)
                                                                                                                gCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGaa                                                       <  1:1567857/131‑1 (MQ=255)
                                                                                                                gCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGaa                                                       <  1:1567718/131‑1 (MQ=255)
                                                                                                                gCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGaa                                                       >  2:1567857/1‑131 (MQ=255)
                                                                                                                gCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTgg                                    <  2:4153918/150‑1 (MQ=255)
                                                                                                                 ccGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGa                                   <  2:1119178/150‑1 (MQ=255)
                                                                                                                   gTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACa                                 >  1:4007540/1‑150 (MQ=255)
                                                                                                                   gTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACa                                 >  1:4318217/1‑150 (MQ=255)
                                                                                                                    tttCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACaa                                <  2:3578436/150‑1 (MQ=255)
                                                                                                                      tCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACg                              >  1:1544898/1‑150 (MQ=255)
                                                                                                                       cgcgACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACgg                             <  1:3885865/150‑1 (MQ=255)
                                                                                                                       cgcgACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGCGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACgg                             >  1:4467210/1‑150 (MQ=255)
                                                                                                                         cgACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGcc                           >  2:5010837/1‑150 (MQ=255)
                                                                                                                          gACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTctgc                                       <  2:2813511/137‑1 (MQ=255)
                                                                                                                          gACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTctgc                                       >  1:2813511/1‑137 (MQ=255)
                                                                                                                             ttCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCcgtcg                       >  1:1374650/1‑150 (MQ=255)
                                                                                                                               ccGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTg                     >  1:1004409/1‑150 (MQ=255)
                                                                                                                               ccGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTg                     <  2:1642609/150‑1 (MQ=255)
                                                                                                                               ccGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTg                     <  2:2302756/150‑1 (MQ=255)
                                                                                                                                    cTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTcgc                >  2:4721567/1‑150 (MQ=255)
                                                                                                                                    cTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTcgc                >  1:854451/1‑150 (MQ=255)
                                                                                                                                       aaCGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGcc             >  2:1398962/1‑150 (MQ=255)
                                                                                                                                       aaCGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGcc             <  1:3491207/150‑1 (MQ=255)
                                                                                                                                           aCCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCa         <  2:2203802/150‑1 (MQ=255)
                                                                                                                                             cTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACt       >  2:3197162/1‑150 (MQ=255)
                                                                                                                                              tGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTg      >  2:1989221/1‑150 (MQ=255)
                                                                                                                                               gTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTgg     >  2:1602191/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTGAAGCTGATGGAAGCTTTCCAGGGCTCGGGCAACCTGCCTGAGTGGATGGTCCTGACCGTCCTGCCAGTGCTGCCGCCGGACCTGCGTCCGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCCGACCTGAACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTC  >  NC_002947/542424‑542719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: