Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	524,158	C→T	L300L (CTG→TTG)	rhsD →	Rhs protein with DUF4329 family putative toxin domain, putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	524,158	0	C	T	100.0%	187.3 / NA	59	L300L (CTG→TTG)	rhsD	Rhs protein with DUF4329 family putative toxin domain, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/28);  total (31/28)

CCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCT  >  NC_000913/524037‑524277                                                                                                                          |                                                                                                                        ccGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCa                                                                                                        <  2:321315/139‑1 (MQ=255)    cACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATAc                                                                                                     >  1:74819/1‑139 (MQ=255)       ctctTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCg                                                                                                  <  1:933344/139‑1 (MQ=255)                  cttcttctGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                  <  1:890700/112‑1 (MQ=255)                  cttcttctGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                  >  2:890700/1‑112 (MQ=255)                           aCAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGcc                                                                              >  2:779035/1‑139 (MQ=255)                               ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                          >  1:424149/1‑139 (MQ=255)                               ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                          <  2:23419/139‑1 (MQ=255)                               ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                          >  2:74316/1‑139 (MQ=255)                               ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                          >  1:1063857/1‑139 (MQ=255)                                tCCCGCCCGCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTgg                                                                         <  2:1063857/139‑1 (MQ=255)                                  ccGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       <  1:74316/139‑1 (MQ=255)                                       ctctctCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTa                                                                  <  2:471599/139‑1 (MQ=255)                                         ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTaca                                                                >  2:751423/1‑139 (MQ=255)                                          tctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATg                                                                                                                    <  1:36314/86‑1 (MQ=255)                                          tctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATg                                                                                                                    >  2:36314/1‑86 (MQ=255)                                                cctcagcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                  <  2:111980/82‑1 (MQ=255)                                                cctcagcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                  >  1:111980/1‑82 (MQ=255)                                                    agcGTTGCCCGACACAATGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCTCCTTTCGGCGGTTTTGTTGATGCACGACCCGGCATACCCGGAGAGTCTGCCCGCTGCGCCACTGGTTCGGTACACGTAAACGGaa                                                     >  2:435435/1‑139 (MQ=255)                                                    agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       >  2:398470/1‑121 (MQ=255)                                                    agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       >  1:1067811/1‑121 (MQ=255)                                                    agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       <  2:1067811/121‑1 (MQ=255)                                                    agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       <  1:398470/121‑1 (MQ=255)                                                    agcGTTCCCCGACACACTGCCCGGTACCAAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCgg                                                                    <  2:56632/124‑1 (MQ=255)                                                    agcGTTCCCCGACACACTGCCCGGTACCAAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCgg                                                                    >  1:56632/1‑124 (MQ=255)                                                                    cTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     <  1:779035/139‑1 (MQ=255)                                                                    cTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     >  2:867552/1‑139 (MQ=255)                                                                     tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATGCGGAAGCCGGTGAACTGATGGc                                    >  2:976295/1‑139 (MQ=255)                                                                     tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     <  1:696372/138‑1 (MQ=255)                                                                     tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     >  2:696372/1‑138 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                              <  2:197132/63‑1 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                              >  1:197132/1‑63 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTg                                                                                  <  2:527033/91‑1 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTg                                                                                  >  1:527033/1‑91 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       >  1:762069/1‑102 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       <  2:762069/102‑1 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     >  2:244321/1‑136 (MQ=255)                                                                       cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     <  1:244321/136‑1 (MQ=255)                                                                            tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGa                                                                                                               <  1:586500/57‑1 (MQ=255)                                                                            tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGa                                                                                                               >  2:586500/1‑57 (MQ=255)                                                                            tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGc                                    <  2:354993/132‑1 (MQ=255)                                                                            tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGc                                    >  1:354993/1‑132 (MQ=255)                                                                               cGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTg                                                                                  <  2:272350/83‑1 (MQ=255)                                                                               cGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTg                                                                                  >  1:272350/1‑83 (MQ=255)                                                                               cGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGAc                          <  2:424149/139‑1 (MQ=255)                                                                                   tACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     >  2:403942/1‑124 (MQ=255)                                                                                   tACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     <  1:403942/124‑1 (MQ=255)                                                                                      ggCCCCGACAGGGGTATCCGCCTTTCGGCGGAGTGGTTGATGCa                                                                                                                  >  2:116031/1‑44 (MQ=255)                                                                                      ggCCCCGACAGGGGTATCCGCCTTTCGGCGGAGTGGTTGATGCa                                                                                                                  <  1:116031/44‑1 (MQ=255)                                                                                        ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAAt                 >  2:438138/1‑139 (MQ=255)                                                                                             aCAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTac                                                                 >  1:679964/1‑86 (MQ=255)                                                                                             aCAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTac                                                                 <  2:679964/86‑1 (MQ=255)                                                                                                  ggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     >  2:748107/1‑109 (MQ=255)                                                                                                  ggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                     <  1:748107/109‑1 (MQ=255)                                                                                                       ccGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCt  <  1:867552/139‑1 (MQ=255)                                                                                                       ccGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCt  >  1:768851/1‑139 (MQ=255)                                                                                                           cTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       >  2:496564/1‑66 (MQ=255)                                                                                                           cTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       <  1:496564/66‑1 (MQ=255)                                                                                                                       tgGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                       >  1:978178/1‑54 (MQ=255)                                                                                                                          |                                                                                                                        CCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCT  >  NC_000913/524037‑524277

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