Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 958,524 C→A 28.5% intergenic (+168/‑25) ykgJ → / → PP_0820 ferredoxin/acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947958,5240CA28.5% 28.0 / 3.5 14intergenic (+168/‑25)ykgJ/PP_0820ferredoxin/acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base A (2/2);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTGCAAGGCTCCTGGCCTGTGGTGAACCTGCCGGCGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCATGCGCCATGCCGATGCCGCGCAGTGGCTGGCGATCATG  >  NC_002947/958425‑958653
                                                                                                   |                                                                                                                                 
aGCGTGCAAGGCTCCTGGCCTGTGGTGAACCTGCCGGCGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGaccccaccc                                                                                               <  2:384319/136‑1 (MQ=255)
   gTGCAAGGCTCCTGGCCTGTGGTGAACCTGCCGGCGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGc                                                                                            >  2:14090/1‑136 (MQ=255)
     gCAAGGCTCCTGGCCTGTGGTGAACCTGCCGGCGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCCCCCCGCAg                                                                                          >  1:242646/1‑136 (MQ=255)
                         tgAACCTGCCGGCGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCACCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGcc                                                                      >  2:257569/1‑136 (MQ=255)
                         tgAACCTGCCGGCGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCACCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGcc                                                                      >  2:338441/1‑136 (MQ=255)
                                     cGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCACCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCAc                                                          <  1:257569/136‑1 (MQ=255)
                                     cGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCACCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCAc                                                          <  1:338441/136‑1 (MQ=255)
                                                       cccGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCa                                        >  1:187004/1‑136 (MQ=255)
                                                       cccGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCa                                        >  1:425899/1‑136 (MQ=255)
                                                                    aaCCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCATGCGCCatgccga                           >  1:352583/1‑136 (MQ=255)
                                                                    aaCCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCATGCGCCatgccga                           >  1:453814/1‑136 (MQ=255)
                                                                                 aCGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCATGCGCCATGCCGATGCCGCGCAGtgg              <  1:476031/136‑1 (MQ=255)
                                                                                         cTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCATGCGCCATGCCGATGCCGCGCAGTGGCTGGCGAt      >  2:372723/1‑136 (MQ=255)
                                                                                             cATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCATGCGCCATGCCGATGCCGCGCAGTGGCTGGCGATCATg  <  1:14090/136‑1 (MQ=255)
                                                                                                   |                                                                                                                                 
AGCGTGCAAGGCTCCTGGCCTGTGGTGAACCTGCCGGCGACGGGCTGCAAAGCTGCCCGAGATGCCTGAACCCTGCACTTGACGTTAAACTGTCATCTCCCCAGCAATGGACAGGAGATCCGGCGTGACCCCACCCCGCAGTTTCCCCAGTGACCTCTGCCTCGACAGCCCACGCCTGCAACTGCGCCCCATGCGCCATGCCGATGCCGCGCAGTGGCTGGCGATCATG  >  NC_002947/958425‑958653

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: