Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,010:1 +C 5.2% intergenic (‑43/+37) pth ← / ← rplY peptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947838,0101.C 5.2% 761.4 / 33.0 172intergenic (‑43/+37)pth/rplYpeptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc
Reads supporting (aligned to +/- strand):  ref base . (73/90);  new base C (6/3);  total (79/93)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 176 total aligned reads displayed.

TGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGAGCTTTCA  >  NC_002947/837861‑838159
                                                                                                                                                      |                                                                                                                                                     
tGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGcgcgcg                                                                                                                                                        >  2:142786/1‑150 (MQ=255)
 gggCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCgg                                                                                                                                                      <  2:3933041/150‑1 (MQ=255)
 gggCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCgg                                                                                                                                                      <  1:831039/150‑1 (MQ=255)
  ggCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGc                                                                                                                                                     >  1:3670609/1‑150 (MQ=255)
     gcTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACt                                                                                                                                                  <  2:2763416/150‑1 (MQ=255)
     gcTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCCGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACt                                                                                                                                                  <  2:3167507/150‑1 (MQ=255)
        ggCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTa                                                                                                                                               <  1:5081170/150‑1 (MQ=255)
         gCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTaa                                                                                                                                              <  2:4212334/150‑1 (MQ=255)
         gCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTCGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTaa                                                                                                                                              <  2:1199149/150‑1 (MQ=255)
          cAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCCGCGCGCACTGTAAc                                                                                                                                             <  2:288005/150‑1 (MQ=255)
           aaGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAAcc                                                                                                                                            <  2:958603/150‑1 (MQ=255)
           aaGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAAcc                                                                                                                                            >  1:4515465/1‑150 (MQ=255)
            aGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAAccc                                                                                                                                           <  1:3650736/150‑1 (MQ=255)
            aGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAAccc                                                                                                                                           >  2:2978342/1‑150 (MQ=255)
            aGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTAGTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAAccc                                                                                                                                           <  2:5118748/150‑1 (MQ=255)
               ccgaataCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGAt                                                                                                                                        >  1:4599270/7‑150 (MQ=255)
               cGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGAt                                                                                                                                        <  2:1794723/150‑1 (MQ=255)
                gTTCAACGAAAAGAGCCCCTGCGTTCTGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCACCTGGAGGGCGGTCACGTCAGGGGCTCGTCCTTGGAGTTGGTGGGGTACAGTGCGCGCGGCACTGTAACCCGATc                                                                                                                                       <  2:4533183/150‑1 (MQ=255)
                  tCAACGAAACGAGCCCCTGCGGTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGGTACCCAGGCCAACGATCACCTGGCTGGCGGTCACGTCAGGGGCTCGTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCaa                                                                                                                                     <  2:4637988/150‑1 (MQ=255)
                    aaCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACg                                                                                                                                   >  1:3147130/1‑150 (MQ=255)
                    aaCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACg                                                                                                                                   >  2:618600/1‑150 (MQ=255)
                     aCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCCATCAACGc                                                                                                                                  >  1:475702/1‑150 (MQ=255)
                      cGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGCTGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGcc                                                                                                                                 <  2:1940026/150‑1 (MQ=255)
                         aaaGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGc                                                                                                                              >  1:1555902/1‑150 (MQ=255)
                             agCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGc                                                                                                                          >  2:2819995/1‑150 (MQ=255)
                               ccccTGCGTTATGCCGGGTCTGTTTGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGa                                                                                                                        >  2:3596415/1‑150 (MQ=255)
                               ccccTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGa                                                                                                                        <  2:357829/150‑1 (MQ=255)
                               ccccTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGa                                                                                                                        >  1:2506337/1‑150 (MQ=255)
                                cccTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCCACGCCGGCGTGCGAg                                                                                                                       >  2:780773/1‑150 (MQ=255)
                                  cTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTg                                                                                                                     >  1:4886035/1‑150 (MQ=255)
                                  cTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTg                                                                                                                     >  1:4885941/1‑150 (MQ=255)
                                  cTGCGTTATGCCGGGTCTGTTCGTAGTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTg                                                                                                                     <  2:3417120/150‑1 (MQ=255)
                                     cGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGAt                                                                                                                  <  1:4363625/150‑1 (MQ=255)
                                            ccGGGTCTGTTCGTATACGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCgg                                                                                                           >  2:823337/1‑150 (MQ=255)
                                              gggTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCa                                                                                                         >  2:34110/1‑150 (MQ=255)
                                                     ttCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCctt                                                                                                  <  2:170459/150‑1 (MQ=255)
                                                      tCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttc                                                                                                 >  1:319438/1‑150 (MQ=255)
                                                      tCGGATTCGGGGCCGGGGGTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGTGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttc                                                                                                 <  1:2090113/150‑1 (MQ=255)
                                                       cGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttct                                                                                                >  2:2318566/1‑150 (MQ=255)
                                                        gTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctg                                                                                               >  1:3944157/1‑150 (MQ=255)
                                                        gTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctg                                                                                               <  2:3227782/150‑1 (MQ=255)
                                                         tATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgc                                                                                              <  2:2165873/150‑1 (MQ=255)
                                                            tCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttc                                                                                           >  1:3788850/1‑150 (MQ=255)
                                                              ggggCCGGGGTTACCCAGGCCAACGATCACCTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                         <  2:1873476/150‑1 (MQ=255)
                                                              ggggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                         >  1:1663229/1‑150 (MQ=255)
                                                              ggggCCGGGGTTACCCAGGCAAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                         <  1:3972922/150‑1 (MQ=255)
                                                                ggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                                       <  2:448517/150‑1 (MQ=255)
                                                                ggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                                       <  1:3039714/150‑1 (MQ=255)
                                                                  ccGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGc                                                                                     >  2:354872/1‑150 (MQ=255)
                                                                       gTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAAcgc                                                                                <  2:73745/150‑1 (MQ=255)
                                                                       gTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAAcgc                                                                                <  2:4807343/150‑1 (MQ=255)
                                                                        ttACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAAcgcg                                                                               <  1:2819995/150‑1 (MQ=255)
                                                                        ttACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAAcgcg                                                                               <  1:2819878/150‑1 (MQ=255)
                                                                          aCCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTg                                                                             <  1:917317/150‑1 (MQ=255)
                                                                           cccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTgg                                                                            >  2:522005/1‑150 (MQ=255)
                                                                            ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa                                                                           >  2:2877758/1‑150 (MQ=255)
                                                                                  cAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGa                                                                     >  1:1243100/1‑150 (MQ=255)
                                                                                   aaCGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGaa                                                                    >  1:3579380/1‑150 (MQ=255)
                                                                                   aaCGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGaa                                                                    >  1:4126134/1‑150 (MQ=255)
                                                                                   aaCGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGaa                                                                    <  1:2728960/150‑1 (MQ=255)
                                                                                       aTCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGtt                                                                <  1:609356/150‑1 (MQ=255)
                                                                                       aTCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGtt                                                                <  2:530512/150‑1 (MQ=255)
                                                                                          aaCTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGTGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGc                                                             >  2:182931/1‑150 (MQ=255)
                                                                                              ggATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACa                                                         <  2:103270/150‑1 (MQ=255)
                                                                                              ggATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACa                                                         >  1:3913782/1‑150 (MQ=255)
                                                                                                   gCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                                    <  1:618600/150‑1 (MQ=255)
                                                                                                    cGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTg                                                   >  1:5202656/1‑150 (MQ=255)
                                                                                                    cGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCATCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCtttt                                                   >  1:5202861/1‑149 (MQ=255)
                                                                                                      gTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGtc                                                 <  2:5036108/150‑1 (MQ=255)
                                                                                                          cGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGtcatca                                             <  1:5271848/150‑1 (MQ=255)
                                                                                                           gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCAc                                            >  2:4976275/1‑150 (MQ=255)
                                                                                                               ggggCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTg                                        <  1:1203966/150‑1 (MQ=255)
                                                                                                               ggggCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTg                                        <  2:4173173/150‑1 (MQ=255)
                                                                                                                  gCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGc                                     <  2:744382/150‑1 (MQ=255)
                                                                                                                  gCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGc                                     <  2:3366127/150‑1 (MQ=255)
                                                                                                                      ttCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTgg                                       >  2:5216330/1‑144 (MQ=255)
                                                                                                                      ttCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTgg                                       <  1:5216330/144‑1 (MQ=255)
                                                                                                                      ttCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCaga                                 >  2:2342352/1‑150 (MQ=255)
                                                                                                                                ttGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAAct                       <  1:4724271/150‑1 (MQ=255)
                                                                                                                                   gtggGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTa                    <  1:3473656/150‑1 (MQ=255)
                                                                                                                                   gtggGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTa                    <  1:4949059/150‑1 (MQ=255)
                                                                                                                                   gtggGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTa                    <  2:1576738/150‑1 (MQ=255)
                                                                                                                                   gtggGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTa                    >  1:2330947/1‑150 (MQ=255)
                                                                                                                                    tggGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTAc                   >  2:3735685/1‑150 (MQ=255)
                                                                                                                                    tggGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTAc                   >  2:1587233/1‑150 (MQ=255)
                                                                                                                                     gggTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACg                  <  1:2877758/150‑1 (MQ=255)
                                                                                                                                     gggTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                                    <  2:3108155/116‑1 (MQ=255)
                                                                                                                                      ggTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTAc                   >  2:4993388/1‑148 (MQ=255)
                                                                                                                                       gTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGcc                >  1:4008768/1‑150 (MQ=255)
                                                                                                                                        ttACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCt               >  1:2778804/1‑150 (MQ=255)
                                                                                                                                           cAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTCCGCCTTTc            <  2:3147130/150‑1 (MQ=255)
                                                                                                                                           cAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTc            <  2:3670609/150‑1 (MQ=255)
                                                                                                                                            aGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCg           <  1:1807825/150‑1 (MQ=255)
                                                                                                                                            aGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCg           <  2:4885941/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCagag                                >  1:999153/1‑127 (MQ=255)
                                                                                                                                              tGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGa         <  2:1339102/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGa         <  2:4435539/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGa         <  1:1887133/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGa         <  2:4355247/150‑1 (MQ=255)
                                                                                                                                              tGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGa         <  1:1066986/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGAGCTTTCA  >  NC_002947/837861‑838159

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: