Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 697,103 T→G intergenic (+27/‑724) mmsA‑I → / → PP_16SD methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947697,1030TG83.3% 11.1 / ‑5.1 13intergenic (+27/‑724)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base G (10/0);  minor base C (2/0);  total (12/1)
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

GCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCA‑‑GGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCG  >  NC_002947/696973‑697219
                                                                                                                                    |                                                                                                                         
gCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGCacgg                                                                                                                      >  1:1987656/1‑133 (MQ=255)
        ggTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGGA‑‑GGGCGcccccg                                                                                                          >  1:106793/1‑134 (MQ=255)
                     ataCCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGGCA‑‑GGCCTCCCCGCGTTGTTTTTgggg                                                                                            >  2:1235832/1‑136 (MQ=255)
                       aCCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGGCG‑‑GGCCGCCCCCCGGTGTTTTTggggg                                                                                           >  1:146006/1‑136 (MQ=255)
                              gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGGCAGCGGCCCCCCCGGTGTGTTTTTggggtggtgg                                                                                      >  1:1187819/1‑130 (MQ=255)
                              gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCCGGCG‑‑GGCGCCCCCCCGGTGTTTTTTGGgggtggttt                                                                                    >  1:1202372/1‑135 (MQ=255)
                              gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGGA‑‑CGCCGCCCCCCGTTTTTTTTTTGGggttgtttt                                                                                    >  1:2039618/1‑135 (MQ=255)
                              gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGCCGGCCG‑‑GGCGCCCCCCCATGGTTTTTTGGggtgggttt                                                                                    >  1:1604936/1‑135 (MQ=255)
                              gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGCCGGCCC‑‑GGCGCCCCCCGCTGTTTTTTGGTggttggttt                                                                                    >  1:2442037/1‑135 (MQ=255)
                                         aCCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCCGGCG‑‑GCCCCCCCCCGGTGGGTTTTGGGGGTGGGTTTTGGTGGTgggg                                                                        >  1:2739203/1‑136 (MQ=255)
                                          cccAGCGTTGGCCGCAGCGGGCCAGCCATGCAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGCC‑‑GGCCGAGCCCGGTTGTTGTTTGGGGGGGGGTTGTTGGGGTGgggg                                                                        >  2:1507288/1‑136 (MQ=255)
                                          cccAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGCA‑‑GGGCGCCCCCCGTTGTTTTTTGGGGGTGGTTGTTTGGGTTggggg                                                                       >  2:1585208/1‑135 (MQ=255)
                                                                                                               ataTGAAGCGGGGAGGCCGTCA‑‑GGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCg  <  2:2130889/134‑1 (MQ=255)
                                                                                                                                    |                                                                                                                         
GCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCA‑‑GGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCG  >  NC_002947/696973‑697219

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