Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,839,221 C→G 16.0% intergenic (+43/‑137) 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,839,2210CG16.0% 70.6 / 7.2 31intergenic (+43/‑137)ccoP‑I/PP_4254cbb3‑type cytochrome c oxidase subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/18);  new base G (2/3);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCGGCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGC  >  NC_002947/4839074‑4839363
                                                                                                                                                    |                                                                                                                                              
ccGGCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGccgccc                                                                                                                                               <  5:372853/149‑1 (MQ=255)
   gcgcAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCcgcg                                                                                                                                            >  3:115006/1‑149 (MQ=255)
      cAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                         <  4:115006/149‑1 (MQ=255)
         gCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                      <  4:127413/149‑1 (MQ=255)
           cGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                    <  7:33175/149‑1 (MQ=255)
                    ggggAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGc                                                                                                                           >  3:303457/1‑149 (MQ=255)
                      ggAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGccc                                                                                                                         <  8:182863/149‑1 (MQ=255)
                       gAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGccc                                                                                                                         >  5:235551/1‑148 (MQ=255)
                         aTGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCccac                                                                                                                      <  6:235551/149‑1 (MQ=255)
                           gACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCccac                                                                                                                      >  1:91712/1‑148 (MQ=255)
                                ggTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAAcc                                                                                                               <  2:91712/149‑1 (MQ=255)
                                 gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                     <  3:209020/126‑1 (MQ=255)
                                 gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                     <  3:82759/126‑1 (MQ=255)
                                 gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                     >  4:209020/1‑126 (MQ=255)
                                 gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                     >  4:82759/1‑126 (MQ=255)
                                       ctgctgGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAAc                                                                                                        <  7:141329/149‑1 (MQ=255)
                                          ctgGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCtt                                                                                                     <  7:80537/149‑1 (MQ=255)
                                             gcggcTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCAt                                                                                                             >  2:330275/1‑138 (MQ=255)
                                             gcggcTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCAt                                                                                                             <  1:330275/138‑1 (MQ=255)
                                             gcggcTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTaca                                                                                                  <  4:303457/149‑1 (MQ=255)
                                             gcggcTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTa                                                                                                    <  8:178788/147‑1 (MQ=255)
                                                      gtgtATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATaa                                                                                                                 >  4:260186/1‑125 (MQ=255)
                                                      gtgtATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCAAAACCAACAAACCTTACACGCcccccc                                                                                         >  6:320725/1‑149 (MQ=255)
                                                      gtgtATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCCCGCGGCGCTCGATATCCGCCCCACCATaa                                                                                                                 <  3:260186/125‑1 (MQ=255)
                                                        gtATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCaaa                                                                                                         <  8:316006/131‑1 (MQ=255)
                                                        gtATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCaaa                                                                                                         >  7:316006/1‑131 (MQ=255)
                                                         tataGCCTGTCGCAGGAT‑GGTTCTGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCccac                                                                                      <  5:72015/149‑1 (MQ=255)
                                                         tataGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCcca                                                                                       <  3:209810/148‑1 (MQ=255)
                                                                                                 gCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCcccttcc                                               <  3:121419/149‑1 (MQ=255)
                                                                                                     aGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTg                                           <  5:283524/149‑1 (MQ=255)
                                                                                                                               gccTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCtt                 <  1:230310/149‑1 (MQ=255)
                                                                                                                                        aTCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATa                                                                    <  1:12821/89‑1 (MQ=255)
                                                                                                                                               ggcgcccgcGCGGCGGTCGTTTTCCGGCCCAACCTAACCCTCAAAACATTCACCCCCGCCCCACCATGCCCCTCAAACCCTTCCCCCACCGGCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGc  <  5:320725/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCGGCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGAT‑GGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGC  >  NC_002947/4839074‑4839363

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: ‑