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%	119.6 / NA	38	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 (20/18);  total (20/18)

GCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCG  >  NC_000913/524039‑524292                                                                                                                        |                                                                                                                                       gCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATa                                                                                                                     <  2:658883/139‑1 (MQ=255)             tATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 <  2:987808/115‑1 (MQ=255)             tATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 >  1:145232/1‑115 (MQ=255)             tATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 <  2:145232/115‑1 (MQ=255)             tATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 >  1:987808/1‑115 (MQ=255)                             ttCCCGCCCTCTCTCAGCCTCATCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 <  2:794370/99‑1 (MQ=255)                             ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 >  1:794370/1‑99 (MQ=255)                              tCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCAGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTgg                                                                                        <  2:1015196/139‑1 (MQ=255)                                       ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAAGACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTaca                                                                               >  1:726566/1‑139 (MQ=255)                                                tcagcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTaca                                                                               <  1:256580/130‑1 (MQ=255)                                                tcagcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTaca                                                                               >  2:256580/1‑130 (MQ=255)                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                      <  2:997761/121‑1 (MQ=255)                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                      >  1:997761/1‑121 (MQ=255)                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGa                                                                     >  1:127422/1‑138 (MQ=255)                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGa                                                                     <  2:127422/138‑1 (MQ=255)                                                                   tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                      <  1:783324/104‑1 (MQ=255)                                                                   tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                      >  2:783324/1‑104 (MQ=255)                                                                   tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGc                                                   >  2:421598/1‑139 (MQ=255)                                                                   tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGc                                                   <  2:393851/139‑1 (MQ=255)                                                                     cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCgg                                                 >  1:430733/1‑139 (MQ=255)                                                                           aCCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 <  1:729983/53‑1 (MQ=255)                                                                           aCCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                                 >  2:729983/1‑53 (MQ=255)                                                                               aaTACGGCCCCGACCGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCGGCCCGCTGCGCCCCTGGTGCGGTACAGGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCg                                       >  1:675399/1‑139 (MQ=255)                                                                               aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCg                                       <  1:421598/139‑1 (MQ=255)                                                                                    ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTa                                                                                 >  1:1002599/1‑92 (MQ=255)                                                                                    ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTa                                                                                 <  2:1002599/92‑1 (MQ=255)                                                                                      ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                                    >  2:1055163/1‑119 (MQ=255)                                                                                      ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTgg                                                    <  1:1055163/119‑1 (MQ=255)                                                                                         cGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGc                                                                                                   <  2:236990/69‑1 (MQ=255)                                                                                         cGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGc                                                                                                   >  1:236990/1‑69 (MQ=255)                                                                                          gACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGc                            >  2:705159/1‑139 (MQ=255)                                                                                                 gTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTgcg                     <  1:46580/139‑1 (MQ=255)                                                                                                            tCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCAc                                                                                           >  1:900505/1‑58 (MQ=255)                                                                                                            tCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCAc                                                                                           <  2:900505/58‑1 (MQ=255)                                                                                                            tCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGt          >  1:865014/1‑139 (MQ=255)                                                                                                                   tgtgGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGAcgc   >  2:1060278/1‑139 (MQ=255)                                                                                                                   tgtgGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGAcgc   >  2:745242/1‑139 (MQ=255)                                                                                                                    gtgGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGAcgcg  <  2:726566/139‑1 (MQ=255)                                                                                                                        |                                                                                                                                       GCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCG  >  NC_000913/524039‑524292

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