Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,481 G→A intergenic (+126/‑229) rnlB → / → yfjP antitoxin RnlB/putative GTP‑binding protein YfjP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,4810GA90.9% 21.5 / ‑3.8 11intergenic (+126/‑229)rnlB/yfjPantitoxin RnlB/putative GTP‑binding protein YfjP
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (4/6);  total (4/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 = 5.20e-01

CGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTAATCATCCTGAGTACCTAGAAGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGT  >  NC_000913/2767238‑2767726
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
cGTTCCGACACATAGCCTGAGCAAACTGACTTTTACTTTACGTAATCATCCTGAGTACCTAGAAGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGGCAAGGGATTTTTGAATCTCTTCTCGCTAAACGAACCAGTATAGAATTCGCGCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTATGCa                                                                                                                                                                                                                                                   <  1:286799/248‑1 (MQ=255)
                                                 ccTGAGTACCTAGAAGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCAATAACCCCCCCTCACAACGTCCTGCATAGTGAACACGTTTGCCCCTGTGCTCCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAAttt                                                                                                                                                                                                <  1:32239/250‑1 (MQ=255)
                                                      gTACCTAGAAGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACACGCACTCATAACCTCATGCATAGGGAAAATGTTTTCCCATATTCACCCACGACTAAATTTAACCACCTTTACACTTATCAATAATGCCACCAATTTTCTgg                                                                                                                                                                                           >  1:270198/1‑250 (MQ=255)
                                                                                                                                 ctAAACGAACCAGTATAGAATTCGCCCGTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCGCACAACGTCCTGCATAGGGAACACGTTTTCCCCTATGCACCCACGTCAAAATTTCCCCCCCGTTAAACGTCTCAATAATGTCACCAATGGTCGGGAAATCCTATGATGCCTTTCCCg                                                                                                                                                                     <  1:409611/197‑1 (MQ=255)
                                                                                                                                                                                                                 cacaACGTCCTGCATAGGGAACACGTTTTCCCCTATGCACCCACGACTATATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAAATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                               <  1:256865/251‑1 (MQ=255)
                                                                                                                                                                                                                 cacaACGTCCTGCATAGGGAACACGTTTTCCCCTATGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCGCCAATTTTCTGGAAAGCCTATGATGCCTTTCCCGGTGAATGGTGGTGTGAGGACAGAAAGATGCTGCAATCTATGTGCATGGAAAAAATCTTCAGATAGATACCATGCTCCTGAACAATCCGATGAACGGCAGAACCTcg                                                           <  1:415273/223‑1 (MQ=255)
                                                                                                                                                                                                                                aGGGAACANGNNNTNNCCTATGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAANACTGTAACCAATTTTNTGNCAATNCTCNGATGc                                                                                                                                                                             >  2:351696/1‑94 (MQ=255)
                                                                                                                                                                                                                                aGGGAACACGGTTTCCCCTATGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                                                                                                                                                             <  1:351696/94‑1 (MQ=255)
                                                                                                                                                                                                                                        cGTTTTCCCCTATGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAAAAATGTCACCAATTTTCGGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCAGTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                          <  1:490497/201‑1 (MQ=255)
                                                                                                                                                                                                                                              ccccTATGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGt  >  1:177275/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ccccTATGCACCCACGACTAAATTTCCCCCCCTTTAAACTTATCAATAATGTCCCCACTTATCTGCAAATCCTATGCTCCCTTTCCCGGCAACTGGCGCTGTAAGTTCAGAAAGATTCTGCAATCTATGTGCATTGACAAAATCTTCAGATAGATATCATTCTCATGAACAACCCGCTGAACGGCAGAACCTCGCGGGTCTGCCGATTTTg                                          >  1:230069/1‑211 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
CGTTCCGACACATATCCTGAGCAAACTGACTTTTACTTTACGTAATCATCCTGAGTACCTAGAAGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGT  >  NC_000913/2767238‑2767726

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