Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,841,191 C→G 8.7% intergenic (+46/‑134) ccoP‑I → / → PP_4254 cbb3‑type cytochrome c oxidase subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,841,1910CG 8.7% 446.2 / 26.3 138intergenic (+46/‑134)ccoP‑I/PP_4254cbb3‑type cytochrome c oxidase subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (50/76);  new base G (2/10);  total (52/86)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 150 total aligned reads displayed.

GGCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTCCGA  >  NC_002947/4841043‑4841340
                                                                                                                                                    |                                                                                                                                                     
ggCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCcgcg                                                                                                                                                      <  2:5764887/150‑1 (MQ=255)
ggCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCcgcg                                                                                                                                                      <  1:750861/150‑1 (MQ=255)
 gcgcAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                                     >  1:956938/1‑150 (MQ=255)
    cAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                                  <  2:7063546/150‑1 (MQ=255)
    cAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                                  <  2:3369182/150‑1 (MQ=255)
    cAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                                  <  2:4159681/150‑1 (MQ=255)
       gCCGAGATCCAGGGGAATGACACGGTGCACCGGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGCCGGCCAAGTGAGGTCGCCGAAAGGCATGTGCCGCCGGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                               <  2:1407744/150‑1 (MQ=39)
       gCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                               >  2:8542292/1‑150 (MQ=255)
       gCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                               <  2:3971221/150‑1 (MQ=255)
        ccGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                              <  1:7126288/150‑1 (MQ=255)
        ccGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                              <  1:5871882/150‑1 (MQ=255)
         cGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                             >  1:181171/1‑150 (MQ=255)
         cGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGAGCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                             >  1:6738547/1‑150 (MQ=255)
          gagaTCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                            <  1:8395438/150‑1 (MQ=255)
          gagaTCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                            <  1:5552634/150‑1 (MQ=255)
          gagaTCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                            <  1:4786556/150‑1 (MQ=255)
          gagaTCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGGGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                            <  1:5552452/150‑1 (MQ=255)
              tCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATc                                                                                                                                        <  2:2554644/150‑1 (MQ=255)
              tCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATc                                                                                                                                        >  1:8176736/1‑150 (MQ=255)
               ccAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctca                                                                                                                                       <  1:6695107/150‑2 (MQ=255)
               ccAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATcc                                                                                                                                       <  1:4740797/150‑1 (MQ=255)
                   gggAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGccc                                                                                                                                   >  1:1394189/1‑150 (MQ=255)
                     gAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCcca                                                                                                                                 >  1:2487041/1‑150 (MQ=255)
                      aaTGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCccac                                                                                                                                <  1:9120775/150‑1 (MQ=255)
                        tGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCAGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCccacca                                                                                                                              >  2:6712774/1‑150 (MQ=255)
                         gACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAAGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCAt                                                                                                                             >  1:2103860/1‑150 (MQ=255)
                          aCAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATa                                                                                                                            >  1:7874340/1‑150 (MQ=255)
                          aCAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATa                                                                                                                            <  1:3426322/150‑1 (MQ=255)
                           cAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATaa                                                                                                                           >  1:8988570/1‑150 (MQ=255)
                            aaGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAAc                                                                                                                          <  2:7764957/150‑1 (MQ=255)
                              ggTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCa                                                                                                                        >  1:6234283/1‑150 (MQ=255)
                               gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCAt                                                                                                                       >  1:2891569/1‑150 (MQ=255)
                               gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCAt                                                                                                                       >  1:8903931/1‑150 (MQ=255)
                               gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCAt                                                                                                                       >  1:2891520/1‑150 (MQ=255)
                                    ccTGCTGGCGGCTTATGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCCCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAAc                                                                                                                  >  2:1826221/1‑150 (MQ=255)
                                    ccTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAAc                                                                                                                  >  2:6357851/1‑150 (MQ=255)
                                    ccTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAAc                                                                                                                  <  1:393847/150‑1 (MQ=255)
                                     ctgctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCCAAcc                                                                                                                 >  1:4002297/1‑150 (MQ=255)
                                     ctgctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCACAcc                                                                                                                 >  1:2352744/1‑150 (MQ=255)
                                      tgctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCt                                                                                                                >  1:4509392/1‑150 (MQ=255)
                                       gctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCtt                                                                                                               <  1:197466/150‑1 (MQ=255)
                                       gctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCtt                                                                                                               <  1:197287/150‑1 (MQ=255)
                                        ctgGCGGCTTACGTGTATAGCCTGTCGCAGTATGGTTCGGCTGAGGCTGTGACGGCCAAGTGTGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTa                                                                                                              <  1:2613973/150‑1 (MQ=255)
                                        ctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTa                                                                                                              <  2:6870296/150‑1 (MQ=255)
                                        ctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTa                                                                                                              <  2:7007371/150‑1 (MQ=255)
                                        ctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTa                                                                                                              <  2:2692287/150‑1 (MQ=255)
                                         tgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTac                                                                                                             >  2:7449600/1‑150 (MQ=255)
                                           gcggcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCCCCATAACCCTCAAACCTTacac                                                                                                           >  1:5212112/1‑150 (MQ=255)
                                             ggcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCACACCTTACAccc                                                                                                         >  2:6297910/1‑150 (MQ=255)
                                              gcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTAcccccc                                                                                                        >  1:4131592/1‑150 (MQ=255)
                                              gcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACAcccc                                                                                                        <  2:2005889/150‑1 (MQ=255)
                                              gcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACAcccc                                                                                                        >  1:4131578/1‑150 (MQ=255)
                                                ttACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCACACCGTACAcccccc                                                                                                      >  1:4848703/1‑149 (MQ=255)
                                                  aCGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCCTCAAACCTTACCCCCCCGcc                                                                                                    >  2:6476237/1‑150 (MQ=255)
                                                  aCGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACAcccccccc                                                                                                    >  2:608256/1‑150 (MQ=255)
                                                    gtgtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGcccc                                                                                                  <  2:8632357/150‑1 (MQ=255)
                                                     tgtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGccccc                                                                                                 >  1:4744843/1‑149 (MQ=255)
                                                      gtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGcccccc                                                                                                >  2:112327/1‑148 (MQ=255)
                                                       tataGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCACACCTTACACcccccccccccc                                                                                               >  2:5647357/1‑150 (MQ=255)
                                                       tataGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACcccccccccccc                                                                                               >  2:8362542/1‑150 (MQ=255)
                                                           gCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATcc                                                                                           <  2:605584/150‑1 (MQ=255)
                                                             cTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATcccc                                                                                         <  2:1462854/150‑1 (MQ=255)
                                                             cTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATcccc                                                                                         <  1:1573105/150‑1 (MQ=255)
                                                             cTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATcccc                                                                                         <  2:8403135/150‑1 (MQ=255)
                                                             cTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATcccc                                                                                         <  1:336803/150‑1 (MQ=255)
                                                             cTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATcccc                                                                                         <  1:6983103/150‑1 (MQ=255)
                                                              tGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATCACCATCAAACCTTACACcccccccccccccccgccc                                                                                        >  2:5221580/1‑143 (MQ=39)
                                                              tGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATccccc                                                                                        <  2:8174690/150‑1 (MQ=255)
                                                                tCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTc                                                                                      <  1:7384735/150‑1 (MQ=255)
                                                                  gCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGGTATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATCACCAGCAACCCTTCCCCCcccgccccccccgcccccgccc                                                                                    >  2:6405156/1‑146 (MQ=255)
                                                                   cAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCCCCCCCCCCGcccccgccca                                                                                   >  2:5461490/1‑145 (MQ=39)
                                                                    aGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCCAcccccgccccc                                                                                  >  2:3279564/1‑144 (MQ=255)
                                                                        tGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATa                                                                              <  2:362458/150‑1 (MQ=255)
                                                                            tCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATAcccc                                                                          <  2:5885237/150‑1 (MQ=255)
                                                                             cGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATAccccc                                                                         <  1:8542292/150‑1 (MQ=255)
                                                                                 tGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCCCCAGCCCCCCCCCACCCAAccccccccc                                                                     >  2:4939730/1‑149 (MQ=38)
                                                                                   aGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAAc                                                                                                                  <  1:6451995/103‑1 (MQ=255)
                                                                                     gCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGATATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                               >  1:5714067/1‑72 (MQ=25)
                                                                                     gCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGATATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                               <  2:5714067/72‑1 (MQ=25)
                                                                                     gCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGATATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                               <  2:1126244/72‑1 (MQ=25)
                                                                                        gtgACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGcccc                                                              <  2:404393/150‑1 (MQ=255)
                                                                                             ggCCACGTGAGTTCGCCGAACGGCATTTGCCGCCTGTGAGCTCGCGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCCTCAACCCTGACCCCCCCGCCCCACCATCCCCCTCAAACCCTACCCCCACCGCCGCCcccttcc                                                         <  2:8789831/150‑1 (MQ=38)
                                                                                             ggCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCCCCAGCCCCCCGCAGCCACTCCCCCCCCCCCGccccccggc                                                         >  2:5332657/1‑146 (MQ=37)
                                                                                              gCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCcccttccc                                                        <  1:2178755/150‑1 (MQ=255)
                                                                                                   gTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCa                                                   <  2:5192538/150‑1 (MQ=255)
                                                                                                    tGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCAc                                                  <  2:1551060/150‑1 (MQ=255)
                                                                                                        ttCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCAccccc                                              <  1:4545355/150‑1 (MQ=255)
                                                                                                          cGCCGAAAGGCCTTTGCCGCCTGTGAGCTCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCCTACCCCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCCTGCCCCCCCAc                                            <  2:5093929/150‑1 (MQ=39)
                                                                                                           gCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCAcc                                           <  1:794250/150‑1 (MQ=255)
                                                                                                                         gccgccTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGc                             <  1:2141274/150‑1 (MQ=255)
                                                                                                                                     aTCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAg                 <  2:6161800/150‑1 (MQ=255)
                                                                                                                                         aGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTggg             <  1:1943302/150‑1 (MQ=255)
                                                                                                                                           cgccgcCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCt           <  1:5247565/150‑1 (MQ=255)
                                                                                                                                            gccgccCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCGTCCCTTGCACCCCCACCGAACACCACTAAGCTTGTTGCCTCAGTGGGCtt          <  2:3595064/150‑1 (MQ=255)
                                                                                                                                               gcccgcGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGc       <  1:1111538/150‑1 (MQ=255)
                                                                                                                                               gcccgcGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGc       <  1:1111369/150‑1 (MQ=255)
                                                                                                                                                 ccgcGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTc     <  1:8669754/150‑1 (MQ=255)
                                                                                                                                                 ccgcGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTc     <  2:4191453/150‑1 (MQ=255)
                                                                                                                                                 ccgcGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTGCCACCCCCGCCCCCCCATCCCCCTCAAACCAGACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTc     <  2:1378220/150‑1 (MQ=255)
                                                                                                                                                    cgcgGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTCCGa  <  2:7874340/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTCCGA  >  NC_002947/4841043‑4841340

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