Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,586,030 2 bp→TC intergenic (+113/+101) PP_4061 → / ← PP_4063 hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,586,0300CT87.5% 15.1 / ‑3.5 9intergenic (+113/+102)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (1/6);  total (1/8)
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
*NC_0029474,586,0310TC88.9% 17.1 / ‑3.1 9intergenic (+114/+101)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (1/7);  minor base G (0/1);  total (1/8)
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.73e-01

CCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCC  >  NC_002947/4585888‑4586116
                                                                                                                                               ||                                                                                        
ccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgg                                                                                     <  4:140430/148‑1 (MQ=255)
                              ggggAGACCACGGCACCGGCGGGAACAGCTGCTCTGCGCTGCTCCATTGAAAGGGTGAGGATGACCGGGGGCTGGGCACGGGGACGCGCGTGACTGGCTCGGGCCGGCCGCGTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCg                                                       <  1:131751/148‑1 (MQ=255)
                                    aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                                  <  8:62460/147‑1 (MQ=39)
                                      cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGaa                                               <  5:141793/148‑1 (MQ=39)
                                                         gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  4:188551/148‑1 (MQ=39)
                                                                 gcgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTTAAACCCCTCCGAAAAGGGGCCCGGGCGGGCGTGAATGCGGCGCGGCCCCCGGGTGGCGCggg                    >  3:52534/1‑148 (MQ=255)
                                                                   gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGc                       <  8:179048/143‑1 (MQ=38)
                                                                   gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg                     <  5:147953/145‑1 (MQ=38)
                                                                                   gtggcggATGACCGGGTGCTTGCCACGGTGACTCTCTTGCATTGATCGGGCCGGCCTCTGCGCGGGTAACCCCGCTCCACCCCGGGACCGGGCCGGCCTGGAATGCGGCGCGGCCCCGGGTAGGAGCGGGTTTACCCGCGAAGAggcc  <  8:105729/143‑1 (MQ=19)
                                                                                                                                               ||                                                                                        
CCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCC  >  NC_002947/4585888‑4586116

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: