Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,719,691 G→A 10.9% G13S (GGC→AGC)  PP_2379 → cytochrome oxidase biogenesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,719,6910GA10.9% 146.2 / 7.0 55G13S (GGC→AGC) PP_2379cytochrome oxidase biogenesis protein
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base A (3/3);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

ATCCCCGCCATGGTCGGGTAGAACCCCTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAACGTGCGCACCCTGTCGA  >  NC_002947/2719547‑2719838
                                                                                                                                                |                                                                                                                                                   
aTCCCCGCCATGGTCGGGTAGAACCCCTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCAt                                                                                                                                                 >  3:67045/1‑149 (MQ=255)
   cccGCCATGGTCGGGTAGAACCCCTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATgg                                                                                                                                               <  5:114364/148‑1 (MQ=255)
     cGCCATGGTCGGGTAGAACCCCTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCg                                                                                                                                            >  8:164786/1‑149 (MQ=255)
                 gTAGAACCCCTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTg                                                                                                                                <  7:164786/149‑1 (MQ=255)
                 gTAGAACCCCTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTg                                                                                                                                >  5:191187/1‑149 (MQ=255)
                        cccTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCctgct                                                                                                                         <  3:156451/149‑1 (MQ=255)
                        cccTCTGGGTGTCATCCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCctgct                                                                                                                         <  7:226897/149‑1 (MQ=255)
                             tGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTAACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAg                                                                                                                    >  2:96415/1‑149 (MQ=255)
                               ggTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGc                                                                                                                  <  1:279481/149‑1 (MQ=255)
                               ggTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGc                                                                                                                  <  2:284204/149‑1 (MQ=255)
                                  gtCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCAt                                                                                                                                                 >  6:126452/1‑115 (MQ=255)
                                  gtCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCAt                                                                                                                                                 <  5:126452/115‑1 (MQ=255)
                                     aTGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAgc                                                                                                            >  7:34661/1‑149 (MQ=255)
                                           agcggCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAgcccggcc                                                                                                      <  6:14329/145‑1 (MQ=255)
                                            aaTGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCAGCATGGGCGTTCTCGGGCTTGGCCTGCTGGc                                                                                                                      >  7:260527/1‑132 (MQ=255)
                                            aaTGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCAGCATGGGCGTTCTCGGGCTTGGCCTGCTGGc                                                                                                                      >  7:260531/1‑132 (MQ=255)
                                            aaTGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCAGCATGGGCGTTCTCGGGCTTGGCCTGCTGGc                                                                                                                      <  8:260535/132‑1 (MQ=255)
                                            aaTGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCAGCATGGGCGTTCTCGGGCTTGGCCTGCTGGc                                                                                                                      <  8:260531/132‑1 (MQ=255)
                                            aaTGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCAGCATGGGCGTTCTCGGGCTTGGCCTGCTGGc                                                                                                                      <  8:260527/132‑1 (MQ=255)
                                            aaTGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCAGCATGGGCGTTCTCGGGCTTGGCCTGCTGGc                                                                                                                      >  7:260535/1‑132 (MQ=255)
                                                 cccGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCgg                                                                                                                                     <  5:124124/112‑1 (MQ=255)
                                                 cccGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCgg                                                                                                                                     >  6:124124/1‑112 (MQ=255)
                                                 cccGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCggggg                                                                                                >  5:59509/1‑149 (MQ=255)
                                                 cccGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCggggg                                                                                                <  6:183924/149‑1 (MQ=255)
                                                  ccGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGc                                                                                               >  8:323861/1‑149 (MQ=255)
                                                    ggTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTc                                                                                             <  6:338038/149‑1 (MQ=255)
                                                      tttttCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGa                                                                                           <  4:67045/149‑1 (MQ=255)
                                                             aCCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGcccggcccgg                                                                                                   <  4:294206/134‑1 (MQ=255)
                                                             aCCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGcccggcccgg                                                                                                   >  3:294206/1‑134 (MQ=255)
                                                                 gCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACg                                                                                >  3:248332/1‑149 (MQ=255)
                                                                   ccGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGc                                                                              <  6:59509/149‑1 (MQ=255)
                                                                      tCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAg                                                                           <  8:134593/149‑1 (MQ=255)
                                                                        ccGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGaa                                                                         >  7:160308/1‑149 (MQ=255)
                                                                          gAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGa                                                                          <  2:152516/146‑1 (MQ=255)
                                                                          gAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGa                                                                          >  1:152516/1‑146 (MQ=255)
                                                                                  gCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAAc                                                               <  8:34661/149‑1 (MQ=255)
                                                                                       ttCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGc                                                                                                               >  1:306801/1‑96 (MQ=255)
                                                                                       ttCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGc                                                                                                               <  2:306801/96‑1 (MQ=255)
                                                                                       ttCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAAt                                                          >  5:75617/1‑149 (MQ=255)
                                                                                          ccTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATcc                                                        <  2:155403/148‑1 (MQ=255)
                                                                                                   gATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCac                                              >  7:264892/1‑148 (MQ=255)
                                                                                                            atgaatgaTCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGcc                                      >  3:275653/1‑148 (MQ=255)
                                                                                                                  gaTCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTc                                                                                       >  5:59494/1‑93 (MQ=255)
                                                                                                                  gaTCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTc                                                                                       <  6:59494/93‑1 (MQ=255)
                                                                                                                      ttttGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCa                                   <  4:77929/141‑1 (MQ=255)
                                                                                                                      ttttGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCa                                   >  3:77929/1‑141 (MQ=255)
                                                                                                                           aCCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAg                       <  4:248332/148‑1 (MQ=255)
                                                                                                                             ccGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGt                     >  3:331235/1‑148 (MQ=255)
                                                                                                                             ccGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGt                     >  2:248457/1‑148 (MQ=255)
                                                                                                                               ggCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAAc                  <  2:67682/149‑1 (MQ=255)
                                                                                                                                 cGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAACGt                <  3:32752/149‑1 (MQ=255)
                                                                                                                                   tGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAACGTgc              <  8:160308/149‑1 (MQ=255)
                                                                                                                                        ttGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAACGTGCGCAcc         >  4:190857/1‑149 (MQ=255)
                                                                                                                                             gCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAACGTGCGCACCCTGTc    <  1:96415/149‑1 (MQ=255)
                                                                                                                                               cGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAACGTGCGCACCCTGTCGa  <  7:323861/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATCCCCGCCATGGTCGGGTAGAACCCCTCTGGGTGTCATGCCAGAATGCCCCGGTTTTTCGACCGGCCCGTCCCGAGGGCCGGCACCTTCCCTGTAAAGGATGGTTTCATGAATGATCTTTTGACCCGGCGTGCGGTTGTCGCCGGCATGGGCGTTCTCGGGCTTGGCCTGCTGGCAGGCTGCAGCCCGGCCCGGGGGCTCGAGTTCAAGTACGGCAAGAACATGAGCAACGAAATCCTCGGGCGCAAGTTCAGCCTCAAAGACCCGCAAGGTAACGTGCGCACCCTGTCGA  >  NC_002947/2719547‑2719838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑