Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	733,230	A→G	R1216R (CGA→CGG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	733,230	0	A	G	97.2%	125.2 / ‑4.1	36	R1216R (CGA→CGG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/18);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACT  >  NC_000913/733170‑733368                                                             |                                                                                                                                           tgatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAg                                                             <  1:1199406/140‑1 (MQ=255) tgatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAg                                                             >  1:139982/1‑140 (MQ=255) tgatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAg                                                             <  1:285992/140‑1 (MQ=255)   atgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGAt                                                           <  2:709171/140‑1 (MQ=255)       ggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCAt                                                       <  2:907805/140‑1 (MQ=255)       ggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCAt                                                       <  2:842600/140‑1 (MQ=255)       ggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCAt                                                       <  1:707338/140‑1 (MQ=255)        gagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCAtt                                                      >  1:414646/1‑140 (MQ=255)         agTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTa                                                     <  1:865347/140‑1 (MQ=255)         agTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTa                                                     >  2:1018419/1‑140 (MQ=255)                     aTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAaca                                         >  1:295677/1‑140 (MQ=255)                       tACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACt                                       <  2:295677/140‑1 (MQ=255)                           aCCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTaaaa                                   <  1:650230/140‑1 (MQ=255)                           aCCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTaaaa                                   <  1:656141/140‑1 (MQ=255)                           aCCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTaaaa                                   >  2:222345/1‑140 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAaca                                         <  2:754999/131‑1 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAaca                                         >  1:754999/1‑131 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                >  2:736795/1‑140 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                >  2:441532/1‑140 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                >  1:929511/1‑140 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                >  1:330592/1‑140 (MQ=255)                                cACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCAt                              >  1:541949/1‑140 (MQ=255)                                       attatGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCa                       <  1:794933/140‑1 (MQ=255)                                             aCCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCAt                 >  2:460132/1‑140 (MQ=255)                                              cccGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCAtt                >  2:644035/1‑140 (MQ=255)                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAAt            >  2:124663/1‑140 (MQ=255)                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTC‑AATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATg           >  1:536509/1‑140 (MQ=255)                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCATTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAAt            >  1:966347/1‑140 (MQ=255)                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAA‑TTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCAc                      >  2:473080/1‑129 (MQ=255)                                                  cTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAA‑TTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCAc                      <  1:473080/129‑1 (MQ=255)                                                     cAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATggg         <  2:138388/140‑1 (MQ=255)                                                     cAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATggg         <  1:460132/140‑1 (MQ=255)                                                        gggCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCgg      <  2:545558/140‑1 (MQ=255)                                                          gCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGaa    <  1:1211157/140‑1 (MQ=255)                                                          gCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGaa    <  1:124663/140‑1 (MQ=255)                                                            gATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACt  <  1:87670/140‑1 (MQ=255)                                                             |                                                                                                                                           TGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACT  >  NC_000913/733170‑733368

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