Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,488,579 (C)7→8 coding (1485/1503 nt) yjgR ← DUF853 domain‑containing protein YjgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,488,5721.C95.0% 61.6 / ‑3.0 20coding (1492/1503 nt)yjgRDUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (10/9);  total (10/10)
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.79e-01

AGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAA  >  NC_000913/4488446‑4488699
                                                                                                                                |                                                                                                                               
aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTcccccc                                                                                                                           >  2:566982/1‑133 (MQ=255)
aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTcccccc                                                                                                                           >  2:355887/1‑133 (MQ=255)
  cAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTAACTGACCTTCTTCccccccccc                                                                                                                         >  1:231689/1‑132 (MQ=255)
                 tCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCCCCCCCCCCAGCAAACTCCCCAAc                                                                                                           <  2:411364/133‑1 (MQ=255)
                        tGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCcgc                                                                                                    >  1:361841/1‑133 (MQ=255)
                                                       tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACGTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAgcggc                                                                     >  1:90726/1‑133 (MQ=255)
                                                       tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAgcggc                                                                     >  1:101931/1‑133 (MQ=255)
                                                         gggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGgcgc                                                                   >  2:338538/1‑133 (MQ=255)
                                                           gtgggCATGATAAGACGCGCCAGCATCGTATCTGGCAAAGAGAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCtt                                                                 <  2:231689/133‑1 (MQ=255)
                                                              ggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTg                                                              >  1:134163/1‑133 (MQ=255)
                                                                 aTGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGcc                                                           <  2:134163/133‑1 (MQ=255)
                                                                  tGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCa                                                          <  1:355887/133‑1 (MQ=255)
                                                                        gACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGttt                                                    >  2:377016/1‑133 (MQ=255)
                                                                          cgcgCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGc                                                  <  2:526809/133‑1 (MQ=255)
                                                                             gcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGtt                                                     <  1:112038/127‑1 (MQ=255)
                                                                             gcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGtt                                                     >  2:112038/1‑127 (MQ=255)
                                                                                  cGTCGCATCAGGCAAAGACAAATCGCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGcaccacacc                                          >  2:149728/1‑133 (MQ=255)
                                                                                                  gACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCccgccg                          <  2:422142/133‑1 (MQ=255)
                                                                                                            gCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCgg                >  1:435500/1‑133 (MQ=255)
                                                                                                             ccTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGt               >  1:499116/1‑133 (MQ=255)
                                                                                                                 aCTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTg           <  1:338538/133‑1 (MQ=255)
                                                                                                                 aCTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTg           <  2:90726/133‑1 (MQ=255)
                                                                                                                        ttcttcTCCCCCCCCAGCAAGCTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACa    <  1:149728/133‑1 (MQ=255)
                                                                                                                          cttctcCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACaaa  <  1:6557/133‑1 (MQ=255)
                                                                                                                                |                                                                                                                               
AGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAA  >  NC_000913/4488446‑4488699

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