Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,475	G→A	L631L (CTG→CTA)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	731,475	0	G	A	100.0%	169.4 / NA	61	L631L (CTG→CTA)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/36);  total (25/36)

ACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCG  >  NC_000913/731339‑731587                                                                                                                                         |                                                                                                                 acacGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcg                                                                                                               <  2:321877/140‑1 (MQ=255)   acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  1:1039142/140‑1 (MQ=255)   acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  2:175847/140‑1 (MQ=255)   acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  2:408458/140‑1 (MQ=255)   acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  2:523048/140‑1 (MQ=255)   acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  1:850628/140‑1 (MQ=255)   acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  1:343228/140‑1 (MQ=255)   acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  1:1232376/140‑1 (MQ=255)     gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  1:427539/138‑1 (MQ=255)     gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             >  2:427539/1‑138 (MQ=255)        gggCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCATAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTa                                                                                                        <  2:196074/140‑1 (MQ=11)        gggCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTa                                                                                                        <  2:1043915/140‑1 (MQ=255)          gCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             <  1:477698/133‑1 (MQ=255)          gCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAAcgcg                                                                                                             >  2:477698/1‑133 (MQ=255)          gCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATg                                                                                                      <  1:341312/140‑1 (MQ=255)          gCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATg                                                                                                      <  1:526688/140‑1 (MQ=255)                 aaCGCGGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACg                                                                                               <  2:585809/140‑1 (MQ=255)                   cgcgGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACACGCAGTATGGAATACGAt                                                                                             >  2:152869/1‑140 (MQ=255)                     cgGTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtgc                                                                                           >  2:329345/1‑140 (MQ=255)                      ggTATGAATACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtgct                                                                                          >  2:425524/1‑140 (MQ=255)                            aaTACAACGCCGCTGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGa                                                                                    <  1:580225/140‑1 (MQ=255)                              tACAACGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACg                                                                                  <  2:1265882/140‑1 (MQ=255)                                  aCGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACgggac                                                                              >  2:1064347/1‑138 (MQ=255)                                  aCGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc                                                                              >  2:718468/1‑140 (MQ=12)                                  aCGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc                                                                              >  1:811354/1‑140 (MQ=12)                                  aCGCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc                                                                              >  2:1020832/1‑140 (MQ=12)                                   cgccgcCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCa                                                                             <  1:641886/140‑1 (MQ=12)                                   cgccgcCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCa                                                                             >  2:1011265/1‑140 (MQ=9)                                                      gccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGa                                                                    >  1:809784/2‑130 (MQ=9)                                                      gccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGa                                                                    <  2:809784/129‑1 (MQ=9)                                                        cGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg                                                        <  1:124380/140‑1 (MQ=17)                                                          tCATTGCCCCGGACGGCAGCAGAAACGGGACACGGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                      <  1:485544/140‑1 (MQ=12)                                                           cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCa                                                     <  1:718468/140‑1 (MQ=17)                                                              tGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGcc                                                  <  1:754196/140‑1 (MQ=17)                                                                    ggACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                      <  2:44644/130‑1 (MQ=17)                                                                    ggACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                      >  1:775602/1‑130 (MQ=17)                                                                    ggACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                      <  2:775602/130‑1 (MQ=17)                                                                    ggACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                      >  1:44644/1‑130 (MQ=17)                                                                     gACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccac                                           >  2:925218/1‑140 (MQ=17)                                                                     gACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccac                                           >  2:1266381/1‑140 (MQ=17)                                                                     gACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccac                                           >  1:69986/1‑140 (MQ=17)                                                                     gACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccac                                           >  1:1099899/1‑140 (MQ=17)                                                                      aCGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACaccacc                                          <  2:103897/140‑1 (MQ=17)                                                                          cagcagAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACGACCTTcc                                      >  2:861255/1‑140 (MQ=12)                                                                          cagcagAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCtt                                        >  2:1184304/1‑138 (MQ=17)                                                                          cagcagAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCtt                                        <  1:1184304/138‑1 (MQ=17)                                                                              agAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTa                                  <  2:976948/140‑1 (MQ=17)                                                                                aaaCGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACg                                <  1:776069/140‑1 (MQ=17)                                                                                aaaCGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGCATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACg                                <  2:1099899/140‑1 (MQ=12)                                                                                          caGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGa                      <  1:955970/140‑1 (MQ=21)                                                                                               cGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGc                 <  1:1020832/140‑1 (MQ=25)                                                                                               cGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGc                 <  1:72895/140‑1 (MQ=25)                                                                                                  tGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGa                      >  1:86359/1‑132 (MQ=255)                                                                                                  tGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGa                      <  2:86359/132‑1 (MQ=21)                                                                                                    cGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGc                 <  1:959443/135‑1 (MQ=25)                                                                                                    cGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGc                 >  2:959443/1‑135 (MQ=255)                                                                                                      tGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACa          >  2:282349/1‑140 (MQ=25)                                                                                                      tGGGGAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACa          >  2:477135/1‑140 (MQ=255)                                                                                                       ggggAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAg         >  1:1107140/1‑140 (MQ=255)                                                                                                       ggggAAAAGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAg         >  2:17568/1‑140 (MQ=255)                                                                                                              aGCCATCTGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCg  <  2:1053581/140‑1 (MQ=18)                                                                                                                                         |                                                                                                                 ACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCG  >  NC_000913/731339‑731587

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