Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	733,260	A→G	K1226K (AAA→AAG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	733,260	0	A	G	100.0%	232.0 / NA	62	K1226K (AAA→AAG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/36);  total (26/36)

TGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAGATATA  >  NC_000913/733170‑733392                                                                                           |                                                                                                                                     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                                                           >  2:222345/1‑140 (MQ=255)                           aCCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTaaaa                                                           <  1:650230/140‑1 (MQ=255)                           aCCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTaaaa                                                           <  1:656141/140‑1 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAaca                                                                 <  2:754999/131‑1 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAaca                                                                 >  1:754999/1‑131 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                                        >  1:330592/1‑140 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                                        >  1:929511/1‑140 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                                        >  2:736795/1‑140 (MQ=255)                              gCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGc                                                        >  2:441532/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                                 <  1:460132/140‑1 (MQ=255)                                                     cAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATggg                                 <  2:138388/140‑1 (MQ=255)                                                        gggCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCgg                              <  2:545558/140‑1 (MQ=255)                                                          gCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGaa                            <  1:124663/140‑1 (MQ=255)                                                          gCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGaa                            <  1:1211157/140‑1 (MQ=255)                                                            ggtatatCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACt                          >  1:17872/3‑140 (MQ=255)                                                            gATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACt                          <  1:87670/140‑1 (MQ=255)                                                                     cTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaaga                 >  2:1011409/1‑140 (MQ=255)                                                                       cAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaagaag               <  1:441532/140‑1 (MQ=255)                                                                       cAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaagaag               <  2:541949/140‑1 (MQ=255)                                                                       cAGGATCCGATTGGGCCGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaagaag               <  2:1135174/140‑1 (MQ=255)                                                                        aGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGaa                            >  1:319685/1‑126 (MQ=255)                                                                        aGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGaa                            <  2:319685/126‑1 (MQ=255)                                                                        aGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGc              >  1:488549/1‑140 (MQ=255)                                                                        aGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGc              >  1:654753/1‑140 (MQ=255)                                                                         ggATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCg             <  1:724370/140‑1 (MQ=255)                                                                          gATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCgg            <  2:214950/140‑1 (MQ=255)                                                                           aTCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGt           <  2:843203/140‑1 (MQ=255)                                                                           aTCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGt           <  1:644035/140‑1 (MQ=255)                                                                            tCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTc          <  1:1028240/140‑1 (MQ=255)                                                                            tCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTc          <  2:414646/140‑1 (MQ=255)                                                                            tCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTC‑AATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTcc         <  2:536509/140‑1 (MQ=255)                                                                               gATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAg       <  2:488549/140‑1 (MQ=255)                                                                               gATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAg       <  1:543957/140‑1 (MQ=255)                                                                                aTTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAGa      >  2:156814/1‑140 (MQ=255)                                                                                  tGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAGata    >  1:89025/1‑140 (MQ=255)                                                                                    ggCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaagaag               >  1:156666/1‑127 (MQ=255)                                                                                    ggCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaagaag               <  2:156666/127‑1 (MQ=255)                                                                                    ggCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAGatata  <  2:457887/140‑1 (MQ=255)                                                                                    ggCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAGatata  <  1:736795/140‑1 (MQ=255)                                                                                      cTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaagaag               <  2:556851/125‑1 (MQ=255)                                                                                      cTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAaagaag               >  1:556851/1‑125 (MQ=255)                                                                                           |                                                                                                                                     TGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCACTGCATTCAATGGGCGGAACTGGTGAAAGAAGCGGTCCAGATATA  >  NC_000913/733170‑733392

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: ‑