Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,917,234 G→T 16.6% A158A (GCC→GCA ccmG ← holocytochrome c synthetase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,917,2340GT16.6% 64.9 / 2.2 24A158A (GCC→GCAccmGholocytochrome c synthetase subunit
Reads supporting (aligned to +/- strand):  ref base G (11/9);  new base T (2/2);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

GATAGGTGTCGATCGCCGCCTTGGCCACACCGGCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGCCCTGCTCGTCACGAATATCCAGCTGAT  >  NC_002947/4917099‑4917361
                                                                                                                                       |                                                                                                                               
gATAGGTGTCGATCGCCGCCTTGGCCACACCGGCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTg                                                                                                                                 >  2:493752/1‑136 (MQ=255)
    ggTGTCGATCGCCGCCTTGGCCACACCGGCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCAt                                                                                                                             >  2:31678/1‑136 (MQ=255)
                  ccTTGGCCACACCGGCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACg                                                                                                               >  2:492076/1‑136 (MQ=255)
                  ccTTGGCCACACCGGCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGACCTGATACAACGTCGCCAGCTGTTCAAGCCAGACCGTGGCATCGACGATGCATACg                                                                                                               >  2:97751/1‑136 (MQ=255)
                               ggCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTTGCATCGACGATGCCTACGTGCTTGTAGCGga                                                                                                  >  2:497175/1‑136 (MQ=255)
                               ggCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTTGCATCGACGATGCCTACGTGCTTGTAGCGga                                                                                                  >  2:12194/1‑136 (MQ=255)
                                                         cTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATcga                                                                                                                          <  2:514776/86‑1 (MQ=255)
                                                         cTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATcga                                                                                                                          >  1:514776/1‑86 (MQ=255)
                                                             cagccagcGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGttt                                                                    <  1:294439/136‑1 (MQ=255)
                                                              agccagcGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCACACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTc                                                                   >  1:397883/1‑136 (MQ=255)
                                                               gccagcGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTcc                                                                  <  2:337413/136‑1 (MQ=255)
                                                                ccagcGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCg                                                                 <  1:493752/136‑1 (MQ=255)
                                                                 cagcGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCgg                                                                <  1:351453/136‑1 (MQ=255)
                                                                          aTGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATa                                                       >  1:107240/1‑136 (MQ=255)
                                                                                 ggCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCCTACACCCCa                                                >  1:87068/1‑136 (MQ=255)
                                                                                   ccTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAg                                              <  2:107240/136‑1 (MQ=255)
                                                                                    cTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAgg                                             >  2:329399/1‑136 (MQ=255)
                                                                                               ggCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTccaggcccag                                  >  1:253489/1‑136 (MQ=255)
                                                                                                     gATACAACGGCGCCAGCTGTTCACGCCAGACCGTTGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGcc                            <  1:497175/136‑1 (MQ=255)
                                                                                                     gATACAACGGCGCCAGCTGTTCACGCCAGACCGTTGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGcc                            <  1:12194/136‑1 (MQ=255)
                                                                                                         cAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCATACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGCCCTGc                        >  1:95739/1‑136 (MQ=255)
                                                                                                             ggCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCGTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGCCCTGCTCGt                    <  2:55845/136‑1 (MQ=255)
                                                                                                              gcgcCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGCCCTGCTCGta                   >  2:202656/1‑135 (MQ=255)
                                                                                                                       gTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGCCCTGCTCGTCACGAATATc          <  1:437659/136‑1 (MQ=255)
                                                                                                                               cAGACCGTGGCATCGACGATGCATACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGCCCTGCTCGTCACGAATATCCAGCTGAt  <  2:95739/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                               
GATAGGTGTCGATCGCCGCCTTGGCCACACCGGCCAGGCTCATGCCCAGCACTGCAGCTGCCAGCCAGCGCCTCATGGCTTGGCCTCGTCGATCAGGCCCTGATACAACGGCGCCAGCTGTTCACGCCAGACCGTGGCATCGACGATGCCTACGTGCTTGTAGCGGATGATGCCCTTGGCATCGATCAGGAAGGTTTCCGGCGCCCCATACACCCCAAGGTCCAGGCCCAGGCTGCCCTGCTCGTCACGAATATCCAGCTGAT  >  NC_002947/4917099‑4917361

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