Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	252,349	G→A	E17K (GAA→AAA)	yafO →	ribosome‑dependent mRNA interferase toxin YafO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	252,349	0	G	A	100.0%	115.7 / NA	32	E17K (GAA→AAA)	yafO	ribosome‑dependent mRNA interferase toxin YafO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/12);  total (20/12)

AATCGCCCCGCAGGATATCTCTTAAGTGCCAGCGCATTCGAAGCGTTAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAG  >  NC_000913/252104‑252589                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 aatCGCCCCGCAGGATATCTCTTAAGTGCCAGCGCATTCGAAGCGTTAATGGACATGCTTGCTGAACAAGAGGAGACAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACtt                                                                                                                                                                                                                                             >  2:202674/1‑251 (MQ=255)      ccccGCAGGATATCTCTTAAGTGCCAGCGCATTCGAAGCGTTAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGc                                                                                                                                                                                                                                        >  1:47500/1‑251 (MQ=255)        ccGCAGGATATCTCTTAAGTGCCAGCGCATTCGAAGCGTTAATGGACATGATTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGt                                                                                                                                                                                                                                      >  1:108853/1‑251 (MQ=255)              gATATCTCGTAAGTGCCAGCGCATTCGAAGCGTTAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGGGATGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAAcgg                                                                                                                                                                                                                                <  1:163076/251‑1 (MQ=255)                            gCCAGCGCATTCGAAGCGTTAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAATGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGTGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTGAAGGAATAAGGATGCGGGTATTCAAAACAAAACTGATTCGCCTGCAACTTTCAGCAGAGAAACTTGATGCGTTAACGGCGGAGTTTATTTc                                                                                                                                                                                                                   <  1:234998/250‑1 (MQ=255)                             ccAGCGCATTCGAAGCGTTAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGAtttttttttct                                                                                                                                                                                                                 >  2:223778/1‑248 (MQ=255)                                           gCGTTAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCTTGACGgt                                                                                                                                                                                                   <  2:257805/251‑1 (MQ=255)                                               tAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGGGttt                                                                                                                                                                                               >  1:207144/1‑251 (MQ=255)                                                aaTGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGtttt                                                                                                                                                                                              <  2:376800/251‑1 (MQ=255)                                                                                    cATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTacgacg                                                                                                                                                           >  1:62581/1‑250 (MQ=255)                                                                                            cgcgcTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGCATTTATTTCCTATAAGCGTGACGGGGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCttt                                                                                                                                                   >  1:414409/1‑250 (MQ=255)                                                                                                               tgcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGGGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcgg                                                                                                                               >  1:152337‑M1/1‑232 (MQ=255)                                                                                                                gcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGGGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTTggaaggtgcgaataagcggg                                                                                                                              >  1:387334‑M1/1‑231 (MQ=255)                                                                                                                        agagGAAATTACACGCCTCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTTTTTTCTATAAGCGTGAAGGTGTTTTGCCAGATATATTTTGTCGCGATGCTCTCTACGACGACTCCTTTAg                                                                                                                                                 >  2:231368‑M1/1‑223 (MQ=255)                                                                                                                        agagGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAg                                                                                                                                                 <  1:231368‑M1/224‑2 (MQ=255)                                                                                                                            gAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGGTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcgggtaaattcttctc                                                                                                                  <  2:414409‑M1/251‑33 (MQ=255)                                                                                                                                    acGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgact                                                                                                          <  2:62581‑M1/251‑41 (MQ=255)                                                                                                                                        cgcgCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggcgcgaataagcggggaaattcttctcggctgactcagt                                                                                                      <  1:159732‑M1/251‑45 (MQ=255)                                                                                                                                                   atatCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTac                                                                                                                                                               <  2:89742/183‑1 (MQ=255)                                                                                                                                                   atatCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTac                                                                                                                                                               >  1:89742/1‑183 (MQ=255)                                                                                                                                                   atatCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggggcgaataagcggggaaattcttctcggctgactcagtcatttcatttc                                                                                           >  1:57783‑M1/1‑196 (MQ=255)                                                                                                                                                         tAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcacttcatttctt                                                                                         <  2:331991‑M1/247‑58 (MQ=255)                                                                                                                                                             gATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttg                                                                                 <  1:202674‑M1/251‑66 (MQ=255)                                                                                                                                                                      gatgatgatTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAACCTTGATGCGTTAACGGCGGATTTTATTTTCTATAAGCGTGCCGGCGGTTTGCCAGCTATCTTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcattttttcatgtttgagccgattt                                                                        >  1:270142‑M1/1‑177 (MQ=255)                                                                                                                                                                                       cTTTAAGGACTAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAAATTGATGCGTTAACGGCGGTTTTTATTTCCTATAAGCGTGACGGGGTTGTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggggcgaataagcggggaaattcttctcggctgactcagtcatttcattttttcatgtttgagccgattttttctcccgtaaatgcc                                                       >  1:72358‑M1/1‑160 (MQ=255)                                                                                                                                                                                          tAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGAGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCCCTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcctttcatttcttcatgtttgagccgattttttctcccgtaaatgccttg                                                    >  1:272653‑M1/1‑157 (MQ=255)                                                                                                                                                                                                              aTTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgt                                >  1:347633‑M1/1‑137 (MQ=255)                                                                                                                                                                                                                 caaaacaaaaCTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgttt                              <  1:292043‑M1/250‑117 (MQ=255)                                                                                                                                                                                                                    aacaaaacTTATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTTCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactccgtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcccatttagaccgtttcttt                          >  1:361080‑M1/1‑131 (MQ=255)                                                                                                                                                                                                                            ttATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGGGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgttttttcgccatttt                  >  1:298312‑M1/1‑123 (MQ=255)                                                                                                                                                                                                                             tATTCGCCTGCAACTTACAGCAGAGAAACTTGATGCGTTAACGGCGGCTTTTATTTCCTATAAGCGTGAAGGTGTTTTGGCAGATATATTTGGTCGCGATGCACTCTACGACGAATCCTTTTggaaggtgcgaataagaggggcaatttttctcggctgactcagtcatttaatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccggttcttcgccattta                  >  1:266600‑M1/1‑122 (MQ=255)                                                                                                                                                                                                                                            tACAGCAGAGAAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCCGATATATTTTGTTGCGATGCACTCTACGACGACTCCTTTAggaaggtgcgaataagcggggaaattcttctcggccgactcagtcattttatttttttatgtttgagtcgattttttttcccgtaaatgccttgaatcagcctatttagaccgtttatttgccatttaatgcgttttccccagt  >  2:6568‑M1/1‑107 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 AATCGCCCCGCAGGATATCTCTTAAGTGCCAGCGCATTCGAAGCGTTAATGGACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAG  >  NC_000913/252104‑252589

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