Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,328,603 G→T 28.5% K66N (AAG→AAT PP_5631 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,328,6030GT28.5% 26.2 / 4.2 14K66N (AAG→AATPP_5631hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base T (2/2);  total (7/7)
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

CGGTGCAGGTACTTCTCTCGACCTACGCCAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCTTGCACAGCCAACGACATACCAACTTTTAGTGAGCCTCAACGTGAAAACATCCCTG  >  NC_002947/4328470‑4328716
                                                                                                                                     |                                                                                                                 
cGGTGCAGGTACTTCTCTCGACCTACGCCAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAATGc                                                                                                                 >  2:43953/1‑136 (MQ=255)
cGGTGCAGGTACTTCTCTCGACCTACGCCAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAATGc                                                                                                                 >  2:43955/1‑136 (MQ=255)
    gCAGGTACTTCTCTCGACCTACGCCAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTc                                                                                                             <  1:561773/136‑1 (MQ=255)
                ctcGACCTACGCCAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCa                                                                                                 >  1:626595/1‑136 (MQ=255)
                 ccGACCTACGCCAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGTGACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAg                                                                                                <  2:620500/135‑1 (MQ=255)
                           ccAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTATGCGCCCTTTAGGCATAAGGAATTCGACAGCAATTCAGCAATACTCCa                                                                                      >  2:478534/1‑136 (MQ=255)
                            cAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCaa                                                                                     >  2:263162/1‑136 (MQ=255)
                                                         cGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCtt                                                        <  1:263162/136‑1 (MQ=255)
                                                         cGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCtt                                                        <  2:596526/136‑1 (MQ=255)
                                                                         ggAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAATGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCTTGCACAGCCAACGACAt                                        <  1:43953/136‑1 (MQ=255)
                                                                         ggAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAATGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCTTGCACAGCCAACGACAt                                        <  1:43955/136‑1 (MQ=255)
                                                                         ggAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCTTGCACAGCCAACGACAt                                        >  2:148663/1‑136 (MQ=255)
                                                                                                        atCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCTTGCACAGCCAACGACATACCAACTTTTAGTGAGCCTCAACGTGAAAAc         <  1:273533/136‑1 (MQ=255)
                                                                                                               aCTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCTTGCACAGCCAACGACATACCAACTTTTAGTGAGCCTCAACGTGAAAACATCCCTg  >  1:661570/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                 
CGGTGCAGGTACTTCTCTCGACCTACGCCAAGTGGTTGAATTCTGCCAACGATTGGACGGAACTGGCGAAGCTGGAAAAGAACGTAATGGGTACAGGATTAAATATCGTTCACTTTGCGCCCTTTAGGCATAAGGCATTCGACAGCAATTCAGCCATACTCCAAAATGCATCGGTTTTAGGGGGAAAACCCTTGCACAGCCAACGACATACCAACTTTTAGTGAGCCTCAACGTGAAAACATCCCTG  >  NC_002947/4328470‑4328716

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