Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622,839	A→G	R1216R (CGA→CGG)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,622,839	0	A	G	100.0%	114.5 / NA	33	R1216R (CGA→CGG)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/11);  total (22/11)

CAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTACTAATGGT  >  NC_000913/3622747‑3622974                                                                                             |                                                                                                                                        cagcagCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTATAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa                                                                                          >  1:252010/1‑140 (MQ=9)           tCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCaa                                                                                >  1:55976/1‑140 (MQ=255)                  gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGata                                                                                    >  1:457249/1‑129 (MQ=255)                  gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGata                                                                                    <  2:457249/129‑1 (MQ=255)                    cGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTATAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCg                                                                       <  2:252010/140‑1 (MQ=255)                                       ggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGAccccc                                                    <  1:426353/140‑1 (MQ=255)                                        gagTCCGGCCTGTATTACAACCGCCACCGCTTTTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCt                                                   >  2:192764/1‑140 (MQ=255)                                        gagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCt                                                   <  2:329893/140‑1 (MQ=255)                                        gagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCt                                                   >  1:313865/1‑140 (MQ=255)                                           tCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTggg                                                >  2:393797/1‑140 (MQ=255)                                           tCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTggg                                                >  2:142215/1‑140 (MQ=255)                                           tCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTggg                                                >  1:86274/1‑140 (MQ=255)                                           tCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTggg                                                >  2:469889/1‑140 (MQ=255)                                                          aaCCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGgga                                 <  1:469889/140‑1 (MQ=255)                                                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTAt                                           <  1:46601/126‑1 (MQ=255)                                                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTAt                                           >  2:46601/1‑126 (MQ=255)                                                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGAt                             >  1:55590/1‑140 (MQ=255)                                                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGAt                             >  1:109535/1‑140 (MQ=255)                                                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGAt                             >  1:466965/1‑140 (MQ=255)                                                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGAt                             >  1:418369/1‑140 (MQ=255)                                                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGAt                             >  1:21010/1‑140 (MQ=255)                                                               ccACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATg                            >  1:528478/1‑140 (MQ=255)                                                                          atGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCgggg                 <  2:313865/140‑1 (MQ=255)                                                                             aCCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCa              >  1:59272/1‑140 (MQ=255)                                                                             aCCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCa              <  1:515942/140‑1 (MQ=255)                                                                               ccGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATg            <  2:109535/140‑1 (MQ=255)                                                                                cGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGt           >  1:190890/1‑140 (MQ=255)                                                                                  ttGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTAc         >  2:33908/2‑140 (MQ=255)                                                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGc               <  2:375263/134‑1 (MQ=255)                                                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGc               >  1:375263/1‑134 (MQ=255)                                                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTAc         >  1:533523/1‑140 (MQ=255)                                                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTAc         >  2:384303/1‑140 (MQ=255)                                                                                         ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTACTAATGGt  <  1:335503/140‑1 (MQ=255)                                                                                             |                                                                                                                                        CAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATATCAATTGAATCCGATATCAGACATCGACCCCCTGGGTTTATCTATGTGGGAGGATGCAAAATCGGGGGCATGTACTAATGGT  >  NC_000913/3622747‑3622974

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