Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,830	G→A	R1216R (CGG→CGA)	rhsA →	rhs element protein RhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,830	0	G	A	100.0%	77.8 / NA	27	R1216R (CGG→CGA)	rhsA	rhs element protein RhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/8);  total (19/8)

CGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCC  >  NC_000913/3765757‑3765958                                                                          |                                                                                                                                 cGGGGCAGCAGTATGATCAGGAGTCCGGCCCGTCTAACATCCGCCACCCCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCa                                                                <  2:412596/140‑1 (MQ=11)   gggCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGt                                                              <  2:525902/140‑1 (MQ=16)       agcagTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACg                                                          <  2:807137/140‑1 (MQ=25)             atgatgAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACa                                                    >  2:20752/1‑140 (MQ=255)              tgatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAGTTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAg                                                   <  2:496849/140‑1 (MQ=255)                   aggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCTCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATACGTGGAATCCAGTTACGAATACAGATCct                                              >  1:379123/1‑140 (MQ=14)                   aggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGCAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCct                                              >  2:589142/1‑140 (MQ=255)                     gagTCCGGCCTGTATGACAACCTCCGCCGTTATTATGGCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCctct                                            <  1:404812/140‑1 (MQ=14)                        tCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTggg                                         >  2:517398/1‑140 (MQ=255)                               tGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTCTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGaa                                  >  2:439028/1‑140 (MQ=255)                                  aTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGtt                               >  2:479686/1‑140 (MQ=255)                                  aTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGtt                               <  2:199905/140‑1 (MQ=255)                                  aTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGtt                               >  1:437609/1‑140 (MQ=255)                                      cAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATTCAGTTACAAATACAGACCCTCCGGAGTTAGAAGTTTTTc                           >  2:323019/1‑140 (MQ=17)                                           gCCACCGCTATTATGACCTGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACATATCTTCGGGGGTTGGAAGTTTTTCCTAGa                      >  1:274680/1‑140 (MQ=17)                                           gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGa                      >  2:114940/1‑140 (MQ=255)                                           gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGa                      >  2:685972/1‑140 (MQ=255)                                           gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGa                      >  2:723208/1‑140 (MQ=255)                                             cACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTATGGGGTTAGAAGTTTTTCCTAGAcc                    >  1:242816/1‑140 (MQ=255)                                              aCCGCTATTATGACCTGTTGCAAGGGCGATATATCTCTCTGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCa                   <  2:274680/140‑1 (MQ=21)                                              aCCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCa                   >  1:604241/1‑140 (MQ=255)                                                   tattatGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTccc              <  1:517398/140‑1 (MQ=255)                                                            ccGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAAt     >  1:20100/1‑140 (MQ=255)                                                             cGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAAtt    >  1:458587/1‑140 (MQ=255)                                                               cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGTGTTCGAATTTTTTCCTGGACCATTCCCCTTGCCAATTcc  >  2:230798/1‑140 (MQ=21)                                                               cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTCGAATTTTTTCCTTGACCGTTCCCCTTGCCAATTcc  >  2:345972/1‑140 (MQ=21)                                                               cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTcc  >  2:701460/1‑140 (MQ=255)                                                                          |                                                                                                                                 CGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCC  >  NC_000913/3765757‑3765958

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