Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,619,272	T→C	R27R (CGT→CGC)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,619,272	0	T	C	100.0%	84.6 / NA	31	R27R (CGT→CGC)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/15);  total (16/15)

TTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCA  >  NC_000913/3619135‑3619356                                                                                                                                          |                                                                                     ttcttcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCAt                                                                                    >  1:280283/1‑140 (MQ=255)    ttcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTgg                                                                                 >  1:11400/1‑140 (MQ=255)    ttcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGAGACATGACGCAGTATGGCGGCAGCAATGTTCAGGGATCAACCGGGGTAAGCATTgg                                                                                 >  1:328248/1‑140 (MQ=25)                 aTGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGc                                                                    <  2:557055/140‑1 (MQ=255)                 aTGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGc                                                                    <  2:280283/140‑1 (MQ=255)                 aTGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGc                                                                    <  1:687175/140‑1 (MQ=255)                     aGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTgg                                                                >  1:213488/1‑140 (MQ=255)                           tGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTgg                                                                >  1:876932/1‑134 (MQ=255)                           tGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTgg                                                                <  2:876932/134‑1 (MQ=255)                              cATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCgg                                                       >  2:884458/1‑140 (MQ=255)                               aTTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTTTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGGGGCCTGTTCCgt                                                      >  1:982060/1‑140 (MQ=255)                               aTTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTTAGGGTTCAGCCGGGGGACGCATTGGTGCCCCCACCGGCGGGGGCCGTTTggg                                                      >  2:250638/1‑139 (MQ=18)                               aTTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGtt                                                          <  2:563049/136‑1 (MQ=255)                               aTTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGtt                                                          >  1:563049/1‑136 (MQ=255)                                ttAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgtg                                                     <  1:250638/140‑1 (MQ=255)                                 tAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgtgt                                                    >  1:975408/1‑140 (MQ=255)                                  aaaaTAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGtgtg                                                   <  1:884458/140‑1 (MQ=255)                                       aataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcg                                              >  1:879593/1‑140 (MQ=255)                                         taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggc                                            >  2:748278/1‑140 (MQ=255)                                              aaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAgg                                       >  1:461614/1‑139 (MQ=255)                                                      tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGTCCg                               >  1:726056/1‑140 (MQ=255)                                                           tGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCAt                          <  1:1001395/140‑1 (MQ=21)                                                             aGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGAAGCATTGTTCATGGGTCAGTCAGGGTACGCATTGGGGCCCCCACCCGCGTCGCCCGTTTGGCGGGCCCCCGCGGCGAGCCCGCCCGCCATcc                        >  2:1057671/1‑140 (MQ=2)                                                               cGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGGCCGGCCATCCgg                      >  1:174869/1‑140 (MQ=11)                                                                  aaaaCCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCa                   <  2:975408/140‑1 (MQ=21)                                                                       cggcggCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATccc              >  1:850077/1‑140 (MQ=25)                                                                       cggcggCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATccc              >  2:329427/1‑140 (MQ=25)                                                                       cggcggCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGACGCATTGGTGCCCCCACCGGCGTGGCCCGTTCGGGGTGCCCCGGCGGGGGGACGTCCCGCCCTCCGGTCAATccc              >  2:359898/1‑140 (MQ=11)                                                                         gcggcgCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCCTTGGTGCCCCCCCCGGCGTGGCCCGTTCGGTGTGCCCCGGGGGAGGGACGTCCCGCCATCCGGTCAATccccc            >  2:719766/1‑139 (MQ=11)                                                                         gcggcgCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCt            <  2:850077/140‑1 (MQ=21)                                                                            gcgcGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCt         <  1:782656/140‑1 (MQ=17)                                                                            gcgcGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCt         <  2:879593/140‑1 (MQ=17)                                                                            gcgcGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCt         <  1:329427/140‑1 (MQ=17)                                                                            gcgcGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCAGGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCt         <  1:1057671/140‑1 (MQ=17)                                                                                   aGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGTGGTACGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCa  <  2:146165/140‑1 (MQ=17)                                                                                                                                          |                                                                                     TTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCA  >  NC_000913/3619135‑3619356

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