Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 233,545 C→T G548G (GGC→GGT PP_0180 → cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947233,5450CT94.4% 48.3 / NA 18G548G (GGC→GGTPP_0180cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/9);  minor base A (1/0);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGCTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGT  >  NC_002947/233420‑233673
                                                                                                                             |                                                                                                                                
gAAGTGCTCCCGTTCTACGAGACCCTGGGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCCGGCGGGGGGCGCGACGGCACTGGTGCTGTTGGGGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGc                                                                                                                        <  2:243794/136‑1 (MQ=255)
        ccTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTa                                                                                                                <  2:233270/136‑1 (MQ=255)
                    gaCCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGttcttc                                                                                                    >  1:105283/1‑136 (MQ=255)
                            ggCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCaa                                                                                            >  1:320176/1‑136 (MQ=255)
                                gcaggcagGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAcgcc                                                                                        <  1:209698/136‑1 (MQ=255)
                                    gcaggcCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCgcgc                                                                                    >  2:354120/1‑136 (MQ=255)
                                    gcaggcCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGAGCTGTTGGTGGGCCTCGCCTGGGTGATTCTGCGGGGATCGGCCAAGCAGCCGCTATCGCTCTACTTCAGCAACAAcgcccccc                                                                                    >  2:209698/1‑132 (MQ=255)
                                         ccccGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGctgctg                                                                               <  2:303682/136‑1 (MQ=255)
                                            cgccgGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTgc                                                                            <  1:20156/136‑1 (MQ=255)
                                            cgccgGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTgc                                                                            <  2:323358/136‑1 (MQ=255)
                                                                tGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGc                                                        >  2:332120/1‑136 (MQ=255)
                                                                tGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGc                                                        >  1:40224/1‑136 (MQ=255)
                                                                          cgcgACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGc                                              >  2:227409/1‑136 (MQ=255)
                                                                          cgcgACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGc                                              >  2:31692/1‑136 (MQ=255)
                                                                                     tGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGtt                                   <  1:31692/136‑1 (MQ=255)
                                                                                      ggTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTg                                  <  2:320176/136‑1 (MQ=255)
                                                                                             ttGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAgg                           <  1:16079/136‑1 (MQ=255)
                                                                                                                      tGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGt  >  1:314482/1‑136 (MQ=255)
                                                                                                                             |                                                                                                                                
GAAGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGCTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGT  >  NC_002947/233420‑233673

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: