Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,439 C→G S116S (TCG→TCC)  yaaI ← DUF2541 domain‑containing protein YaaI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,4390CG92.3% 73.1 / ‑1.9 27S116S (TCG→TCC) yaaIDUF2541 domain‑containing protein YaaI
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base G (12/12);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCG  >  NC_000913/11190‑11683
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
ggATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCg                                                                                                                                                                                                                                                     <  2:167415/251‑1 (MQ=255)
      ccTTCCATCGATTTAAACATTTCGTTGCGTATCAGTCCGCTGGAGATGCGAGTGTTGCCGTTTTCATTATGGGTGAGCAATCGGGCGAAATTGGCCAACTGTTCCTCTCTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGTCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTGCACCGTATTACCTGAAAAgg                                                                                                                                                                                                                                               <  2:382089/251‑1 (MQ=255)
                            cGTTGCGTATCAGTACGCTGTAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCATATCTTAATCATTCAGGTAATTAACATGCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCGAAGAGGTGATTTTTGAGACGCAGCGtt                                                                                                                                                                                                                          <  1:429831/250‑1 (MQ=255)
                             gTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTa                                                                                                                                                                                                                        <  1:41473/251‑1 (MQ=255)
                                gcgtaTCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAAGCGGGCGACATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAAGTTTCAGCGTG‑GCCATATCCGATGATTTCACCCGATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTAtt                                                                                                                                                                                                                      <  1:427005/246‑1 (MQ=255)
                                  gtaTCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTc                                                                                                                                                                                                                   <  1:455425/249‑1 (MQ=255)
                                      cAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTa                                                                                                                                                                                                               <  1:55003/251‑1 (MQ=255)
                                            gCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTg                                                                                                                                                                                                         <  1:353144/251‑1 (MQ=255)
                                                              gCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCcagtca                                                                                                                                                                                       <  1:479964/251‑1 (MQ=255)
                                                                           tGGGTGAGCAATCGGGCGAAATTTGCCAACTGATCCTCACTACAATGCTGCAGAAAATCCAGATCTGAGTCCTTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCctt                                                                                                                                                                          >  2:186508/1‑251 (MQ=255)
                                                                                                       aCTGTTCCTCACTACAATGCTGAAGAAACTCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGAGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAgact                                                                                                                                              >  1:349491/1‑248 (MQ=255)
                                                                                                                  cTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTCTATTATGGCGTTAGTCGTCGCCGTTAATTTTCCGCGTG‑GCCATATCCGATGAGTTCACCGTCTGACCGGAAAAGGTGATTTTTGAGACGAAGCGTTTATTGTCGTTATCCCTGTTACTGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTCTTTCTGAAGGATTATTCAGACTCTGACTGgcgc                                                                                                                                   >  2:482718/1‑251 (MQ=255)
                                                                                                                             aagaaAAGCCAGATCTGAATCATGCAGGTAATTAACATTCATTTTTTGTGGTTGCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGGG‑GCCATATCCGATGAGTTCACCGTATGAACGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCt                                                                                                                        <  2:161444/251‑1 (MQ=255)
                                                                                                                              agaaAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCtt                                                                                                                       >  1:174731/1‑251 (MQ=255)
                                                                                                                                  aaTCCAGATCTGAATCATTCAGGTACTTAAATTTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCTTCGCCGATAAGTTTCAGCGTG‑GCCATAGCCGATGAGTTCACCGTATGACCGTAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTTATCCAGGCAGTGGTTTGCCCTTCTTTTAGTTCGGAAGGAATATGCAGACTCTGACTGGCGCTACGGGCGGCTTTGaa                                                                                                                   <  1:398725/251‑1 (MQ=255)
                                                                                                                                  aaTCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGTCGGCTTTGa                                                                                                                    >  1:169598/1‑250 (MQ=255)
                                                                                                                                              aaTCATTCAGGGAATTAACATTCATTTTTTGTGGTTTCTATATTCTTGCGTTAGTCGTCGCCGATAATTTTCAGC‑TGTGCCATATCCGATTAGTTCACCGTATGACCGGAAAATGTGATTTTTGAGACGCAGCGTTTATTTTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGTGCGGCTTTGAAATAAACCGATGc                                                                                                       <  1:454853/251‑1 (MQ=255)
                                                                                                                                                              aaCATTCATTTTTTGTGGTTTATATATTCTGGCGTTCGTCGTAGCCGATACTTTTCAGCGTG‑GCCATATCCGATGAGTTAACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCATTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAGCCGATGCACTGCTTAACTgcaa                                                                                        >  1:389602/1‑247 (MQ=255)
                                                                                                                                                                     attTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCAt                                                                                >  1:402309/1‑251 (MQ=255)
                                                                                                                                                                              ggtttCTATATTCTGGCGTTAGTCGTCTCCGATAATTTTCAGCGTG‑GCCATATCCGATTAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGTTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCa                                                                        <  1:467771/247‑1 (MQ=255)
                                                                                                                                                                                     atatTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTCTAAATCGCCATGGTCGGCAGAGAGTTg                                                                >  2:332934/1‑251 (MQ=255)
                                                                                                                                                                                                         gCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGAc                                            <  2:133913/251‑1 (MQ=255)
                                                                                                                                                                                                           cGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACa                                           >  2:165601/1‑250 (MQ=255)
                                                                                                                                                                                                                aTTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAgg                                     >  1:213035/1‑251 (MQ=255)
                                                                                                                                                                                                                           g‑gCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAgc                           <  1:377716/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aCCGTATGACCGGACAAGGCGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATACAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGCCTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGATTAACTGTAAATCGCCATGGACGGCAGAGAGTTGTATTCGTTTCACAATGCGACAAACAGGACGTCTAAGCGCCAGATCGTTGGTTTCGtt       >  2:471351/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    aTGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTAcgcg  >  1:274365/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
GGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCG  >  NC_000913/11190‑11683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑