Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,527,451 G→A S29S (AGC→AGT sgcE ← putative epimerase SgcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,527,4510GA91.7% 20.5 / NA 12S29S (AGC→AGTsgcEputative epimerase SgcE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (4/7);  minor base C (1/0);  total (5/7)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CGAAGGATAGTCCAGCGTTTCTGCATGAACGAAGATCCACGCGGGTCTGATTTCAGCCAGCGCGTTAAACCAGCGCTGTGGCCGGGCAAGCATGAAATGAAAGGACAAGGGATGAGGCGTCTGTCTGGCGACGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGAAGCTGCTGTCCTCAATATCCAGATGCAGGCTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGC  >  NC_000913/4527250‑4527634
                                                                                                                                                                                                         |                                                                                                                                                                                       
cGAAGGATAGTCCAGCGTTTCTGCATGAACGAAGATCCACGCGGGTCTGATTTCAGCCAGCGCGTTAAACCAGCGCTGTGGCCGGGCAAGCATGAAATGAAAGGACAAGGGATGAGGCGTCTGTCTGGCGACGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGAAGCTGCTGTCCTCAATATCCAGATGCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGc                                                                                                                                                        >  1:106197/1‑235 (MQ=255)
                                                      agccagcGCTTTAAACCAGCGCTGTGGCCGGGCAAGCATGAAATGAAAGGACAAGGGATGAGGCGTCTGTCTGGCGACGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGAAGCTGCTGTCCTCAATATCCAGATGCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAAcc                                                                                  <  1:15612/251‑1 (MQ=255)
                                                                               ngCCGGGCAAGCATGAAATGAAAGGACAAGGGATGAGGCGTCTGTCTGGCGACGGCCCGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGCAGCTGCTGTCCTCCATATCCAGATGCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGCGGATTTGCGCTGGCGCGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGAGCAAACAGTGCATCg                                                         >  1:30222/2‑251 (MQ=255)
                                                                                gCCGGGCAAGCATGAAATGAAAGGACAAGGGATGAGGCGTCTGTCTGGCGACGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGCAGATGCTGTCCTCAATATCCAGATGCAGCCCCCCACAATCCAGATTGTCCa                                                                                                                                                                   >  1:147044/1‑144 (MQ=255)
                                                                                                              gATGAGGCGTCTGTCTGGCGACGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGAAGCTGCTGTCCTCAATATCCAGATGCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTTATACCAGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGc                          <  1:113652/251‑1 (MQ=255)
                                                                                                                              ggCGACGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGAAGCTGCTGTCCTCAATATCCAGATGCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGCGCAAACAg                                                                >  1:101288/1‑197 (MQ=255)
                                                                                                                                   cGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTGAATGAAGCTGCTGTCCTCAATATCCAAATGCAGACTCCCAAAATCCAGATTGGCCAATTCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCAt                                                                                                   <  1:146553/157‑1 (MQ=255)
                                                                                                                                      ccTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGAAGCTGCTGTCCTCAATATCCAGATGCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGCGCAAACAGTGCATCGTCTGGGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAAGCAATAAACTGGc  <  1:89118/251‑1 (MQ=255)
                                                                                                                                                    aTCCCAAACGTAATGTTGTGAATGAAGCTGCTGGCCTCAATATCCAGATGCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTg                                                                                                                               <  1:230337/112‑1 (MQ=255)
                                                                                                                                                                        nangangCTGCNGTCCNCAATATCCAGATGCANACTCCCAATATCCAGATTTNNCNATGCGNTCAGc                                                                                                                                                        <  2:106197/66‑1 (MQ=255)
                                                                                                                                                                                                     gCAGACTCCCAAAATCCAGATTGTCCAATGCCNTCAGCTCACNNCCNTAGTTNANNNGNNNTGNGCTGGCGAGAGAAGGGTGCNGGNNCATAAGNGCCTCCGTTAACCTGATACCAGCGCa                                                                     <  2:61249/121‑1 (MQ=255)
                                                                                                                                                                                                     gCAGACTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGCGCa                                                                     >  1:61249/1‑121 (MQ=255)
                                                                                                                                                                                                         |                                                                                                                                                                                       
CGAAGGATAGTCCAGCGTTTCTGCATGAACGAAGATCCACGCGGGTCTGATTTCAGCCAGCGCGTTAAACCAGCGCTGTGGCCGGGCAAGCATGAAATGAAAGGACAAGGGATGAGGCGTCTGTCTGGCGACGGCCTGCACCGTTTTCATCCCAAACGTAATGTTGTTAATGAAGCTGCTGTCCTCAATATCCAGATGCAGGCTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGC  >  NC_000913/4527250‑4527634

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