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%	177.3 / NA	55	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 (26/29);  total (26/29)

TGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACG  >  NC_000913/524025‑524290                                                                                                                                      |                                                                                                                                     tGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATg                                                                                                                                 <  1:952510/139‑1 (MQ=255)            cccGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgac                                                                                                                      <  1:259731/139‑3 (MQ=255)                     ctTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCAt                                                                                                                    >  1:87703/1‑132 (MQ=255)                     ctTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCAt                                                                                                                    <  2:87703/132‑1 (MQ=255)                             tgttcttctGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGTGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGc                                                                                                     >  2:808201/3‑139 (MQ=255)                                 cttctGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGc                                                                                                 <  1:481524/139‑1 (MQ=255)                                 cttctGACAGTTCCCGCCCTCTCGCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGc                                                                                                 <  1:1215969/139‑1 (MQ=255)                                   tctGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTg                                                                                               <  1:808201/139‑1 (MQ=255)                                       aCAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGcc                                                                                           >  1:344264/1‑139 (MQ=255)                                           ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACAGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                                           <  2:883078/103‑1 (MQ=255)                                           ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACAGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                                           >  1:883078/1‑103 (MQ=255)                                            tCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                       >  2:999687/1‑106 (MQ=255)                                            tCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                       <  1:999687/106‑1 (MQ=255)                                              ccGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    <  2:344264/139‑1 (MQ=255)                                                                agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATa                                                                                                                   <  2:782302/90‑1 (MQ=255)                                                                agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATa                                                                                                                   >  1:782302/1‑90 (MQ=255)                                                                agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGaa                                                                  >  1:746482/1‑139 (MQ=255)                                                                                   cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCgg                                               >  2:130194/1‑139 (MQ=255)                                                                                        tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTatat                                          >  1:1153029/1‑139 (MQ=255)                                                                                         aCCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    <  2:988204/96‑1 (MQ=255)                                                                                         aCCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    >  1:988204/1‑96 (MQ=255)                                                                                           cGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    <  2:241521/94‑1 (MQ=255)                                                                                           cGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    >  1:241521/1‑94 (MQ=255)                                                                                             aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                       >  2:29862/1‑57 (MQ=255)                                                                                             aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                       <  1:29862/57‑1 (MQ=255)                                                                                             aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCGTGCCCGCTGCGCCACTtgtg                                                                                    <  1:913951/92‑1 (MQ=255)                                                                                             aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCGGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCg                                     <  2:1153029/139‑1 (MQ=255)                                                                                             aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    >  2:913951/1‑92 (MQ=255)                                                                                             aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCg                                     <  2:746482/139‑1 (MQ=255)                                                                                                  ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCgg                                                                                 >  1:848768/1‑90 (MQ=255)                                                                                                  ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCgg                                                                                 <  2:848768/90‑1 (MQ=255)                                                                                                  ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCgcag                                  <  1:717599/137‑1 (MQ=255)                                                                                                  ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCgcag                                  >  2:717599/1‑137 (MQ=255)                                                                                                    ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGt                                                                          >  1:163823/1‑95 (MQ=255)                                                                                                    ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGt                                                                          <  2:163823/95‑1 (MQ=255)                                                                                                    ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGa                                                                   <  2:1190083/102‑1 (MQ=255)                                                                                                    ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGa                                                                   >  1:1190083/1‑102 (MQ=255)                                                                                                    ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAAt                              >  2:452707/1‑139 (MQ=255)                                                                                                       cGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    >  1:535797/1‑82 (MQ=255)                                                                                                       cGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    <  2:535797/82‑1 (MQ=255)                                                                                                       cGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTac                                                                              <  1:389011/88‑1 (MQ=255)                                                                                                       cGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTac                                                                              >  2:389011/1‑88 (MQ=255)                                                                                                                  tCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATa                                                                                                                   <  1:116257/40‑1 (MQ=38)                                                                                                                  tCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATa                                                                                                                   >  2:116257/1‑40 (MQ=38)                                                                                                                  tCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACg                           >  2:1338099/1‑128 (MQ=255)                                                                                                                  tCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACg                           <  1:1338099/128‑1 (MQ=255)                                                                                                                   ccGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                       >  2:41484/1‑35 (MQ=255)                                                                                                                   ccGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                       <  1:41484/35‑1 (MQ=37)                                                                                                                       cTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTg                                                                                                      <  2:587132/48‑1 (MQ=255)                                                                                                                       cTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTg                                                                                                      >  1:587132/1‑48 (MQ=255)                                                                                                                        tttCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    <  2:825147/65‑1 (MQ=255)                                                                                                                        tttCGGCGGTGTGGTTGATGAACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                                    >  1:825147/1‑65 (MQ=255)                                                                                                                         ttCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGt                                                                                     >  2:547020/1‑63 (MQ=255)                                                                                                                         ttCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGt                                                                                     <  1:547020/63‑1 (MQ=255)                                                                                                                                gtgtGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGAcg  <  1:970043/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     TGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACG  >  NC_000913/524025‑524290

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