Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,839,224 C→G 16.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,839,2240CG16.7% 66.3 / 9.3 30intergenic (+46/‑134)ccoP‑I/PP_4254cbb3‑type cytochrome c oxidase subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/18);  new base G (2/3);  total (9/21)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTCCG  >  NC_002947/4839077‑4839372
                                                                                                                                                    |                                                                                                                                                    
gcgcAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCcgcg                                                                                                                                                     >  3:115006/1‑149 (MQ=255)
   cAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                                  <  4:115006/149‑1 (MQ=255)
      gCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                               <  4:127413/149‑1 (MQ=255)
        cGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                             <  7:33175/149‑1 (MQ=255)
                 ggggAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGc                                                                                                                                    >  3:303457/1‑149 (MQ=255)
                   ggAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGccc                                                                                                                                  <  8:182863/149‑1 (MQ=255)
                    gAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGccc                                                                                                                                  >  5:235551/1‑148 (MQ=255)
                      aTGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCccac                                                                                                                               <  6:235551/149‑1 (MQ=255)
                        gACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCccac                                                                                                                               >  1:91712/1‑148 (MQ=255)
                             ggTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAAcc                                                                                                                        <  2:91712/149‑1 (MQ=255)
                              gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                              <  3:82759/126‑1 (MQ=255)
                              gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                              >  4:209020/1‑126 (MQ=255)
                              gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                              <  3:209020/126‑1 (MQ=255)
                              gTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                              >  4:82759/1‑126 (MQ=255)
                                    ctgctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAAc                                                                                                                 <  7:141329/149‑1 (MQ=255)
                                       ctgGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCtt                                                                                                              <  7:80537/149‑1 (MQ=255)
                                          gcggcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCAt                                                                                                                      >  2:330275/1‑138 (MQ=255)
                                          gcggcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCAt                                                                                                                      <  1:330275/138‑1 (MQ=255)
                                          gcggcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTaca                                                                                                           <  4:303457/149‑1 (MQ=255)
                                          gcggcTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTa                                                                                                             <  8:178788/147‑1 (MQ=255)
                                                   gtgtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATaa                                                                                                                          >  4:260186/1‑125 (MQ=255)
                                                   gtgtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCAAAACCAACAAACCTTACACGCcccccc                                                                                                  >  6:320725/1‑149 (MQ=255)
                                                   gtgtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCCCGCGGCGCTCGATATCCGCCCCACCATaa                                                                                                                          <  3:260186/125‑1 (MQ=255)
                                                     gtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCaaa                                                                                                                  >  7:316006/1‑131 (MQ=255)
                                                     gtATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCaaa                                                                                                                  <  8:316006/131‑1 (MQ=255)
                                                      tataGCCTGTCGCAGGATGGTTCTGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCccac                                                                                               <  5:72015/149‑1 (MQ=255)
                                                      tataGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCcca                                                                                                <  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)
                                                                                                                                                    cgcgGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTCCg  <  4:76864/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCAGGCCGAGATCCAGGGGAATGACAAGGTGCACCTGCTGGCGGCTTACGTGTATAGCCTGTCGCAGGATGGTTCGGCTGAGGCTGTGACGGCCAAGTGAGTTCGCCGAAAGGCATTTGCCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATATCCGCCCCACCATAACCATCAAACCTTACACCCCCGCCCCACCATCCCCCTCAAACCATACCCCCACCGCCGCCCCCTTCCCTTGCACCCCCACCGAACACAACTAAGCTTGTTGCCTCAGTGGGCTTCGCTCCG  >  NC_002947/4839077‑4839372

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: