Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	732,962	C→A	A1127E (GCG→GAG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	732,962	0	C	A	100.0%	136.5 / NA	53	A1127E (GCG→GAG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (26/27);  total (26/27)

GATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGAT  >  NC_000913/732873‑733095                                                                                          |                                                                                                                                      gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTc                                                                                     >  2:713863/1‑140 (MQ=14) gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTc                                                                                     <  1:35671/140‑1 (MQ=17) gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTc                                                                                     >  2:89712/1‑140 (MQ=9)   tGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCAc                                                                                   <  1:1127230/140‑1 (MQ=18)      tCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTg                                                                                <  2:738427/140‑1 (MQ=18)       cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTc                                                                                     <  2:146028/134‑1 (MQ=17)       cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTc                                                                                     >  1:146028/1‑134 (MQ=14)         aCCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTc                                                                                     <  2:612053/132‑1 (MQ=17)         aCCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTc                                                                                     >  1:612053/1‑132 (MQ=14)         aCCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTAc                                                                             >  2:913675/1‑140 (MQ=17)         aCCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTAc                                                                             >  1:323131/1‑140 (MQ=17)                  aaaGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGAc                                                                    <  2:345579/140‑1 (MQ=17)                    aGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCa                                                                  >  2:1143862/1‑140 (MQ=17)                      tGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATc                                                                <  1:130676/140‑1 (MQ=18)                        aaaTCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATcgc                                                              >  2:1157169/1‑140 (MQ=17)                           tCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGc                                                           <  1:1033299/140‑1 (MQ=17)                            ccTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGcc                                                          <  1:1143862/140‑1 (MQ=17)                            ccTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGcc                                                          <  1:825585/140‑1 (MQ=17)                            ccTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGcc                                                          <  1:1219032/140‑1 (MQ=17)                             cTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCt                                                         <  2:455826/140‑1 (MQ=17)                                gCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGcc                                                      >  1:806131/1‑140 (MQ=17)                                    aCCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTg                                                  <  2:806131/140‑1 (MQ=17)                                       gggTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGcc                                               <  1:643815/140‑1 (MQ=17)                                         gtgagtgaGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCt                                             <  2:83417/140‑1 (MQ=17)                                          tgagtgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCtt                                            >  1:1156256/1‑140 (MQ=17)                                          tgagtgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCtt                                            >  1:1292090/1‑140 (MQ=17)                                           gagtgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTa                                           >  2:1243857/1‑140 (MQ=17)                                              tgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCa                                        >  1:416906/1‑140 (MQ=17)                                              tgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCa                                        >  1:265169/1‑140 (MQ=17)                                                  gAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCGTCGCGGCCTGCCGCTGGCCCTTATCAGCAc                                    >  2:833267/1‑140 (MQ=16)                                                  gAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCAc                                    >  2:640826/1‑140 (MQ=18)                                                   aaaGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACg                                   >  2:349102/1‑140 (MQ=18)                                                   aaaGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACACCGGCGCGAAAAATTAACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGAACg                                   >  2:365688/1‑140 (MQ=11)                                                        cgccgcTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAgg                              >  1:417770/1‑140 (MQ=21)                                                          ccgcTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGa                            >  1:1027861/1‑140 (MQ=21)                                                           cgcTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGac                           <  1:349102/140‑1 (MQ=25)                                                           cgcTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGac                           >  2:278719/1‑140 (MQ=21)                                                           cgcTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGac                           >  2:1136553/1‑140 (MQ=21)                                                                  tggcATCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCg                    <  2:116158/140‑1 (MQ=255)                                                                       tCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGc                     >  1:894245/1‑134 (MQ=16)                                                                       tCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGc                     <  2:894245/134‑1 (MQ=255)                                                                       tCGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTa               <  2:588229/140‑1 (MQ=255)                                                                        cGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTAc              <  1:1243857/140‑1 (MQ=255)                                                                        cGTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTAc              >  2:858853/1‑140 (MQ=255)                                                                         gTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGtgcg             <  2:1356422/140‑4 (MQ=255)                                                                         gTGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACg             <  1:89712/140‑1 (MQ=255)                                                                          tGCGGCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGc            >  1:1323996/1‑140 (MQ=255)                                                                              gCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGaa        <  1:1263305/140‑1 (MQ=255)                                                                              gCCTGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGaa        <  2:265169/140‑1 (MQ=255)                                                                                 tGACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATAc     <  2:1156256/140‑1 (MQ=255)                                                                                  gACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACg    <  1:858853/140‑1 (MQ=255)                                                                                  gACCGTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACg    <  2:1027861/140‑1 (MQ=255)                                                                                   accgTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGa   >  1:1322194/4‑140 (MQ=255)                                                                                    ccgTGGAGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGAt  >  1:61334/3‑140 (MQ=255)                                                                                          |                                                                                                                                      GATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGAT  >  NC_000913/732873‑733095

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