Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,489,389 G→A 24.9% A225A (GCC→GCT yjgR ← DUF853 family protein with NTPase fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,489,3890GA24.9% 24.9 / 3.1 16A225A (GCC→GCTyjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base A (2/2);  total (9/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

CAGCAGATGTGCCTCGTCGAAGAAAAACACCAGTTTTGGTTTCTCCAGATCGCCCGCTTCCGGCAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTGCGCCTTGCTGTTCCAGCGACAACAGCCCGCGCTGGATGGCACCAACCGATGCGCT  >  NC_000913/4489251‑4489535
                                                                                                                                          |                                                                                                                                                  
cagcagATGTGCCTCGTCGAAGAAAAACACCAGTTTTGGTTTCTCCAGATCGCCCGCTTCCGGTAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGa                                                                                                                                          >  4:12690/1‑149 (MQ=255)
             tcgtcgAAGAAAAACACCAGTTTTGGTTTCTCCAGATCGCCCGCTTCCGGCAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGcc                                                                                                                             <  3:73774/149‑1 (MQ=255)
                       aaaaCACCAGTTTTGGTTTCTCCAGATCGCCCGCTTCCGGCAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTg                                                                                                                   <  4:13209/149‑1 (MQ=255)
                                 ttttGGTTTCTCCAGATCGCCCGCTTCCTGCAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGTCATCGGTgc                                                                                                         <  4:87097/149‑1 (MQ=255)
                                              aGATCGCCCGCTTCCGGCAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCttt                                                                                            <  4:47971/149‑1 (MQ=255)
                                                     ccGCTTCCGGCAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATcc                                                                                     >  3:23080/1‑149 (MQ=255)
                                                                aaTTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCAGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTc                                                                          >  2:9295/1‑149 (MQ=255)
                                                                aaTTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCAGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCAACGGCTc                                                                          >  4:40309/1‑149 (MQ=255)
                                                                                    aaGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTgcgc                                                      >  2:54914/1‑149 (MQ=255)
                                                                                     aGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCAGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTGCGcc                                                     <  1:9295/149‑1 (MQ=255)
                                                                                     aGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCAGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCAACGGCTCACCAAAGAAGTGTGCTGCGcc                                                     <  3:40309/149‑1 (MQ=255)
                                                                                                                 gTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTGCGCCTTGCTGTTCCAGCGACAACAGCCCgcgc                         >  3:1930/1‑149 (MQ=255)
                                                                                                                          tCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTGCGCCTTGCTGTTCCAGCGACAACAGCCCGCGCTGGATGGCa                >  4:42837/1‑149 (MQ=255)
                                                                                                                              aTAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTGCGCCTTGCTGTTCCAGCGACAACAGCCCGCGCTGGATGGCACCaa            >  2:54896/1‑149 (MQ=255)
                                                                                                                                  aGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTGCGCCTTGCTGTTCCAGCGACAACAGCCCGCGCTGGATGGCACCAACCGa        <  3:12690/149‑1 (MQ=255)
                                                                                                                                        tcGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCTTCGTTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGAGCTGCGCGTTGCTGTTCCAGCGACAACAGCCCGCGCTGGATGGCACCAACCGATCCGCt  >  2:5667/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CAGCAGATGTGCCTCGTCGAAGAAAAACACCAGTTTTGGTTTCTCCAGATCGCCCGCTTCCGGCAATTGTTCATACAACTCTGAAAGCATCCACAGCAGGCTGGCGGCGTACAGTTTCGGCATCTGATAAAGTTTCTCGGCGCTGAGGATATTGATAACGCCTTTACCGTTGGCATCGGTGCGCATCCAGTCTTTGATATCCAGCATCGGCTCACCAAAGAAGTGTGCTGCGCCTTGCTGTTCCAGCGACAACAGCCCGCGCTGGATGGCACCAACCGATGCGCT  >  NC_000913/4489251‑4489535

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