Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,854 A→G intergenic (‑68/‑309) yaaI ← / → dnaK DUF2541 domain‑containing protein YaaI/chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,8540AG92.3% 30.1 / NA 13intergenic (‑68/‑309)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (8/4);  minor base T (0/1);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAA  >  NC_000913/11616‑12070
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                        
gagTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTa                                                                                                                                                                                                              <  2:59455/251‑1 (MQ=255)
 agTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTGCAGCGCCAGATCGTTGGTTTCGTTACGCGGCAGTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAggttcggaa                                                                                                                                                                                                             <  2:75349/251‑8 (MQ=255)
                     gCGACAAACAGGAGGTTTCAGAGCCCGGTCGTTGGTTTCGTTACGCGGCATTGCAGTGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAAACGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGCTTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTTTTTTTATTACAAt                                                                                                                                                                                         >  2:114825/1‑251 (MQ=255)
                                                                                                   gTTTGCCTGCGCCGTGCAGCACAGCATCAGGCGAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCat                                                                                                           <  2:151691/251‑1 (MQ=255)
                                                                                                    tttGCCTGCGCCGTGCAGCACAGCATCAGGCTACTCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCat                                                                                                           >  2:84552/1‑250 (MQ=255)
                                                                                                                  gcagcaCAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTg                                                                                            >  1:149111/1‑251 (MQ=255)
                                                                                                                                          ccggCTGGCGGAAATCGTAAAAACGGATTTCATACGGATTCTCTTAGTGGGAAGCGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAACTTTTTg                                                                     >  2:228742/3‑250 (MQ=255)
                                                                                                                                              cTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCAtct                                                                >  2:117869/1‑251 (MQ=255)
                                                                                                                                              cTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCAtct                                                                >  2:81326/1‑251 (MQ=255)
                                                                                                                                                ggCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTcc                                                              <  1:90857/251‑1 (MQ=255)
                                                                                                                                                         cGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATTAATACCCACTAGTTTACTGCTGATAAAGAGAAGTTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTTATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTGTTGCATCTCCCCCTtgatg                                                     <  1:228742/251‑1 (MQ=255)
                                                                                                                                                                                                        ggATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCt      >  1:47701/1‑251 (MQ=255)
                                                                                                                                                                                                             aaTACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCaa  >  1:253755/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                        
GAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAA  >  NC_000913/11616‑12070

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