Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,351,540 G→A 16.9% intergenic (+6/‑36) acsA‑II → / → PP_4703 acetyl‑CoA synthetase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,351,5400GA16.9% 49.5 / 3.2 24intergenic (+6/‑36)acsA‑II/PP_4703acetyl‑CoA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base A (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 = 1.00e+00

CGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACC  >  NC_002947/5351421‑5351669
                                                                                                                       |                                                                                                                                 
cggcgTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGgccg                                                                                                                   >  2:459058/1‑136 (MQ=255)
    gTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCt                                                                                                               <  1:156294/136‑1 (MQ=255)
             gcAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGCGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTAcagcag                                                                                                      >  2:342326/1‑136 (MQ=255)
              aaaaaTCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGTGCGTGGTAGAGGGCCTGATCGACTAACGCCTCTACCAGTAACCGGCGGCCGCAATTCTGGGCGCCCTGTACAGCAgg                                                                                                     <  1:342326/135‑1 (MQ=255)
               aaaaTCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCACCCAGTAACCGGCAGCCGCAAATCTGGCCGCCCTGTACAGCAGGt                                                                                                    >  1:114357/1‑136 (MQ=255)
               aaaaTCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGt                                                                                                    >  1:286364/1‑136 (MQ=255)
               aaaaTCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCAGCCGCAAATCTGGCCGCCCTGTACAGCAGGt                                                                                                    >  1:116963/1‑136 (MQ=255)
                aaaTCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTAcagca                                                                                                       >  1:368334/1‑132 (MQ=255)
                aaaTCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTAcagca                                                                                                       <  2:368334/132‑1 (MQ=255)
                                  tGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTtcgt                                                                                 >  2:54662/1‑136 (MQ=255)
                                       aaCCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCTGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCg                                                                            >  1:287794/1‑136 (MQ=255)
                                                                ccGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTc                                                   >  2:543152/1‑136 (MQ=255)
                                                                    cccGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCg                                               <  1:543152/136‑1 (MQ=255)
                                                                                  aGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGc                                 >  2:371990/1‑136 (MQ=255)
                                                                                            cGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGAcc                       >  1:458107/1‑136 (MQ=255)
                                                                                            cGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGAcc                       >  1:450897/1‑136 (MQ=255)
                                                                                              aCAAACGCCTCAACCAGTACCCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCTTATCTAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                     <  2:286364/136‑1 (MQ=255)
                                                                                              aCAAACGCCTCAACCAGTACCCGGCACCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGTAAGCCCTTCGTCGCCGTATCGAAACCCCGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                     <  2:114357/136‑1 (MQ=255)
                                                                                              aCAAACGCCTCAACCAGTAACCGGCTGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCTACCTg                     <  2:287794/136‑1 (MQ=255)
                                                                                              aCAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                     <  2:128930/136‑1 (MQ=255)
                                                                                              aCAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                     <  2:458107/136‑1 (MQ=255)
                                                                                              aCAAACGCCTCAACCAGTAACCGGCAGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGACGCCGAATCTAAACCCCGGTCATGAGCCCCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                     <  2:116963/136‑1 (MQ=255)
                                                                                                         aaCCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCaa          <  2:450897/136‑1 (MQ=255)
                                                                                                                 aCCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCACAAGGGCGAcc  >  1:30609/1‑136 (MQ=255)
                                                                                                                       |                                                                                                                                 
CGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATACCTCGACCCTGGCCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACC  >  NC_002947/5351421‑5351669

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: