Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,096 A→C 13.9% intergenic (+57/+95) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,0960AC13.9% 404.8 / 9.8 166intergenic (+57/+95)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (86/49);  new base C (7/18);  total (96/70)
Fisher's exact test for biased strand distribution p-value = 1.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

GCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCC  >  NC_002947/667948‑668244
                                                                                                                                                    |                                                                                                                                                     
gcgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAg                                                                                                                                                      <  1:844389/150‑1 (MQ=255)
  gTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGTTAGGGCCGCACAGCg                                                                                                                                                    >  1:4490063/1‑150 (MQ=255)
   tGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCAAAGCgg                                                                                                                                                   >  2:2065116/1‑150 (MQ=255)
     ccgccgccgCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                                                <  2:3913400/150‑1 (MQ=255)
          gccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaaa                                                                                                                                                       >  1:414315/1‑139 (MQ=255)
          gccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaaa                                                                                                                                                       <  2:414315/139‑1 (MQ=255)
           ccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGTGGCGCTAGGGCCGCCAATCGGccccccca                                                                                                                                          >  1:1766181/1‑146 (MQ=255)
             gcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGTGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAt                                                                                                                                        >  1:2680262/1‑150 (MQ=255)
             gcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTTCGGCGCGATAGGGCCCCAAAGCGGCCCCCACaaa                                                                                                                                        >  2:15912/1‑149 (MQ=255)
              cAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGCGCGGTAGGGCCGCCAAGCGGGCCCGCCaaat                                                                                                                                       >  1:286812/1‑148 (MQ=255)
              cAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGGGCGCTAGGGCCGCCAAGCGGCCCCGAccaaa                                                                                                                                       >  1:286728/1‑146 (MQ=255)
                 gTCGAATCTTGGTACCGGCCGAACTTGCTGGGGGTGGGCTTTGGGAATGGGTAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTAGGGGTAAAGTGTGGGAGCGGCCTTGCGTGGCTTTTGGGCCGCAAAGGGGCCCCGACAATTCaa                                                                                                                                    <  1:3484442/150‑1 (MQ=255)
                 gTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCaa                                                                                                                                    <  1:3005459/150‑1 (MQ=255)
                  tCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAg                                                                                                                                   >  2:4468035/1‑150 (MQ=255)
                   cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAgc                                                                                                                                  >  1:3377332/1‑150 (MQ=255)
                   cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATAAAgc                                                                                                                                  >  2:2223705/1‑150 (MQ=255)
                         gATGTACGGGCCTAACCTGCTGGCGGTGGGCGTTGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGGGGCCTTGCGTCGCTATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                                            <  1:3853847/150‑1 (MQ=255)
                          aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                                           <  2:2669877/150‑1 (MQ=255)
                          aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                                           >  2:2639768/1‑150 (MQ=255)
                               cgggCCGCACCTGCTGGCGGTGGGCGATGGGCATGGGGAGTTGACGGTGAACGCCCGCTTGAGCCCCGCTCTAGGGGTAATGTGTGGGTGCGGCGTTGCGTCGCGATTGGGCCGAAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGaa                                                                                                                      <  1:3815249/148‑1 (MQ=255)
                               ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGCGCGATAGGGCCGCAAAGCGGCCCCCACCAGTCAAGCGGCGAAGCTGaa                                                                                                                      >  1:3747798/1‑150 (MQ=39)
                               ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGTCGCGGTTGGGCCGCCACGCGGCCCCGACAATTGAAGCGGCGACGCGGaa                                                                                                                      >  1:2546690/1‑150 (MQ=255)
                                  gCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGGGCGGTAGGGCCGCAAAGCGGCCCCCCCAAATCAACCGGCGAAGCAGAAAGc                                                                                                                   >  1:2329152/1‑150 (MQ=255)
                                    cGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCt                                                                                                                 <  1:1752494/150‑1 (MQ=255)
                                     gacccTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                                <  2:3759867/147‑1 (MQ=255)
                                     gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                                >  2:286036/1‑150 (MQ=255)
                                     gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTGGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                                <  2:1991216/150‑1 (MQ=255)
                                      acccTGCTGGCGGTGGGCGTTGGGCATGGGGAGATGCGGGTGAACGCCCGCTAGCGACCCGCTCTCCGGGGAATGGGTGGGCGGGGCCTGGCGTCGCGTTCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                               <  1:2011663/148‑1 (MQ=255)
                                      aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                               <  2:2500319/150‑1 (MQ=255)
                                       aCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGACCATTAAAGCGGCGAAGCTGAAAGCCgggg                                                                                                              >  2:210075/1‑150 (MQ=255)
                                         cTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTggggg                                                                                                            >  2:2761655/1‑149 (MQ=255)
                                            cAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctg                                                                                                          <  2:1802241/150‑1 (MQ=255)
                                             aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                        >  1:754883/1‑150 (MQ=255)
                                             aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                        <  2:4269284/150‑1 (MQ=255)
                                             aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                        <  2:3200969/150‑1 (MQ=255)
                                              ggCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGctgctg                                                                                                       <  2:4280532/150‑1 (MQ=255)
                                                 cGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCCCAAAGCCGCCCCCACAAATCAAGCGGCGCAGCGGAAAGCCGGGGGCTGCTGTGc                                                                                                    >  2:2288292/1‑150 (MQ=255)
                                                  gTGGGCTCTGGTCATGGGGAGAGGTGGGTGACCGCCCGCTGGAGAACGGCTGGCCGCGGAAATTGTGGTTGCGGCCTTGGTACGCGATAGGGCCGCAAGGCGGGCCCGACATTTCAAGCTGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                   <  2:2099809/150‑1 (MQ=255)
                                                  gTGGGCGATGGTCAAGGGGATATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCTGAAAGTTTGGGAGCGTCATTGATTCGCGATAGGGCCTCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGACTGGGGCTGCTGTGCa                                                                                                   <  2:1912218/150‑1 (MQ=255)
                                                    gggCGTTGGGAATGGGGAGATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGGGTAATGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGAAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                 <  1:523813/150‑1 (MQ=37)
                                                    gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                 <  2:3953850/150‑1 (MQ=255)
                                                    gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                 <  2:4043823/150‑1 (MQ=255)
                                                     ggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                                <  2:754883/150‑1 (MQ=255)
                                                          tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                                           >  2:1412839‑M2/1‑148 (MQ=255)
                                                          tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTGGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                                           >  2:3435005‑M2/1‑148 (MQ=255)
                                                              cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                                       >  2:2184759‑M2/1‑144 (MQ=255)
                                                               aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCACAATTCAATCGGGGCACCGGAAAGCCAGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                                      <  2:2858110‑M2/150‑8 (MQ=255)
                                                               aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                                      <  2:2999241‑M2/150‑8 (MQ=255)
                                                               aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                                      <  2:3343301‑M2/150‑8 (MQ=255)
                                                               aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCCCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgac                                                                                      <  2:4468727‑M2/150‑8 (MQ=255)
                                                                  gggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgca                                                                                   <  2:2567316‑M2/150‑11 (MQ=255)
                                                                  gggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCCAAAGGGCCCCAAACCCGCCCCCACCAATCAAGCCGCCAACCGGAAAGCCAGGGGCGGCGGTGCCCCCCAGGCGCCACGcc                                                                                   >  2:61050/1‑149 (MQ=255)
                                                                   ggAGTTGCCGGTTAACGCCCGCTAGAGCCCCGCTCTACGGGGAAAGTGGGGGCGCGGCCTTGCGTCGCTTTCGGGCCCAAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                                                  >  1:1495633‑M2/1‑139 (MQ=255)
                                                                     agaTGACGGTGAAAGCCCGGTAGTGACGCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGTTAGGGCGGCATAGCGGCCCCGACAATTCAAGCGGCGATGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgctagg                                                                                >  2:366768‑M2/1‑137 (MQ=255)
                                                                     agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                                >  2:366592‑M2/1‑137 (MQ=255)
                                                                              tGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCCGGCCCATAGGGCCGCAAACCCGCCCCCCAAATTAAAGCGGCGGAGCGGAAAGCCAGGGGCTGGGGGGCCGCCCAGTCGCGACCCCCGGCCGCACCca                                                                       >  2:122044/1‑150 (MQ=255)
                                                                              tGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCCAAAATTCAAGCCGCGGAGCGGAAAGCCGGGGGCTGCGGTGCCGCCCAGgcgcgacgccaggcccctcccc                                                                       <  2:121742‑M2/150‑23 (MQ=255)
                                                                                              aCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAACCGGCGAAGCTGAAAGCCAGGGGCTGCTGTGCAGCCCAGTCGCGACGCAACGCCGCTCCCCCAAAGACCGCGCCGaa                                                       >  2:3433234/1‑149 (MQ=25)
                                                                                                 cGCCCTCCGCGTAATGTGTGGGAGCGGCCTTGGGGCGCGATAGGGCCGCAACGCGGCCCCGACAATTCACGCGGCGACGCTGAAAGCCTGGGGCGGCTGTGCAGCCCTGTCGCGCCGCACGGCCGCTCCCACAAAGACCGCGCAGATCAg                                                    <  2:2862236/150‑1 (MQ=255)
                                                                                                  gCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGGTAGGGCCGCAAAGCGGCCCCGGCAAATCAATCGGCGAAGCGGTAAGCCGGGGGCTGCGGTTCCGCCCGGTCGCCACGCAAGGCCGCACCCACAAAAAACGCCCacaagagt                                                   >  2:897886/1‑143 (MQ=19)
                                                                                                      tACGCGGAAAGTGTGGGAGCGGGCTTGCGGCGCGGTAGGGGCCCCAAGCCGCCCCGCCAAATCAAGCGGCGCATCGGCAAGCCTGGGGCGGGTGTGCAGCCCAGTCGCGACGCAAGGCCGCTCCCCCCAAGACCGCGCCGAACCGTGGAt                                               >  2:4189728/1‑150 (MQ=255)
                                                                                                      tACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGGGGAGCCCGGTCGCGACGCAAGGCCGCTCCCACAAAGACCCCGCAGATGAGTggct                                               >  2:127219/1‑148 (MQ=25)
                                                                                                         gcgGAAAGTGTGGGAGCGGCCTTTCGTCGCGGTATGGCCCCCAAGCGGCCCCCACAAATCAAGCGGCGGAGCGGAAAGCCGGGGGGGGCGGTGCCGCCCGGGCGCCGCGCAAGGCCCCTCCCCCAAAGACCCCCCCGAACAGTgggtaac                                            >  2:461388/1‑145 (MQ=255)
                                                                                                            tactgtggGGGCGGGGCCTTGCGTGGCTATCGGGCCGCAAAGCGGCCCCGAAAATTCAAGCGGCTAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCa                                         <  2:1083183/146‑1 (MQ=25)
                                                                                                             tctttgtgggagCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                                       >  1:1271088‑M2/5‑97 (MQ=255)
                                                                                                                          ggCCTGGCGTCGCGATTGGGCCGAAACGCGGCCCCGACAATTCACGGGGCGACGCTGAACGCCGGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCcaca                           <  2:3377350/150‑1 (MQ=38)
                                                                                                                          gccgtggCGTGGCTATCGGGCCGCAACGCGGCCCCTACAATTCACGCGGCGAAGCTGAATGCGTGGGGGTGCTGGGCCGCCCTGTCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCcaca                           <  2:735996/144‑1 (MQ=255)
                                                                                                                           gCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg                          >  2:2457464/1‑150 (MQ=255)
                                                                                                                           gCCTGGGGTCGCTATGGGCCCCCAACGCGCCCCCCACAATTCACGGGGCTACGCTTAAAGCGTGGGGCTGCTGTGCCGCCCGGTCGCGACTAAGGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg                          <  2:61077/150‑1 (MQ=255)
                                                                                                                            ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGc                         >  2:2637212/1‑150 (MQ=255)
                                                                                                                            ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGc                         >  2:2972252/1‑150 (MQ=255)
                                                                                                                              ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc                       >  2:4245825/1‑150 (MQ=255)
                                                                                                                              ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc                       >  2:4099979/1‑150 (MQ=255)
                                                                                                                                 cGTCGCGATTGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAg                    <  2:3747798/150‑1 (MQ=255)
                                                                                                                                  gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc                   >  1:3431844/1‑150 (MQ=255)
                                                                                                                                  gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCAACACCCACAGCCCCAGc                   >  1:4320172/1‑150 (MQ=255)
                                                                                                                                     gcgATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAcc                >  2:657934/1‑150 (MQ=255)
                                                                                                                                     gcgATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAcc                >  2:148092/1‑150 (MQ=255)
                                                                                                                                      cgATCGGGCCGCAACGGGGCCCCGCCAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCGGTGCAGCCGTGGCGCGACGAACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAcca               <  2:3527255/150‑1 (MQ=255)
                                                                                                                                      cgATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgccc               <  1:3332763‑M2/150‑79 (MQ=255)
                                                                                                                                      cgATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAcca               >  2:4356128/1‑150 (MQ=255)
                                                                                                                                       gATCGGGCCCCAACGCGGCCCCGACAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccag              <  2:2852092/150‑1 (MQ=255)
                                                                                                                                       gATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgccca              <  1:1846408‑M2/150‑80 (MQ=255)
                                                                                                                                       gATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccag              >  1:3227516/1‑150 (MQ=255)
                                                                                                                                          tgggCCGCACCGGGGCCCCTACAATCCACGCGCCGAAGCTTAAAGCGTGGGGCTGGTGTGCTGCCGTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca           <  2:67485/149‑1 (MQ=255)
                                                                                                                                            ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccatcgcc         <  1:1315823‑M2/150‑85 (MQ=255)
                                                                                                                                            ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga         >  1:154138/1‑150 (MQ=255)
                                                                                                                                            ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga         >  2:3548173/1‑150 (MQ=255)
                                                                                                                                            ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga         >  2:1538391/1‑150 (MQ=255)
                                                                                                                                             gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgc                                                             >  2:3718447‑M2/1‑65 (MQ=255)
                                                                                                                                             gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCa                     >  1:3338688/1‑137 (MQ=255)
                                                                                                                                             gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCa                     <  2:3338688/137‑1 (MQ=255)
                                                                                                                                             gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgat        >  2:1831152/1‑150 (MQ=255)
                                                                                                                                              ccGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatg       >  2:4043855/1‑150 (MQ=255)
                                                                                                                                              ccGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatg       >  1:1035267/1‑150 (MQ=255)
                                                                                                                                               cGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc                   <  1:1635202/137‑1 (MQ=255)
                                                                                                                                               cGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc                   >  2:1635202/1‑137 (MQ=255)
                                                                                                                                                gCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatgat     >  2:4133641/1‑150 (MQ=255)
                                                                                                                                                gCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatgat     >  2:1364810/1‑150 (MQ=255)
                                                                                                                                                   aaGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTcc  >  1:907504/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCC  >  NC_002947/667948‑668244

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