Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,455,246	C→T	C66C (TGC→TGT)	unknown →	IncQ plasmid conjugative transfer protein TraQ (RP4 TrbM homolog)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,455,246	0	C	T	100.0%	41.8 / NA	14	C66C (TGC→TGT)	unknown	IncQ plasmid conjugative transfer protein TraQ (RP4 TrbM homolog)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/4);  total (10/4)

ACCGGTGATACCAAAACCGCCTGTGAAGCTCTACTCTGCTTATCATCACCAAAAGACGCGGCGAAAACGCCGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAACTTTTTAAAACTATGCCCGGATGCTAGTGATGATTCAGTGAACCAGGCTTTATCTACAAAATGCGCCCGCAAGCCGTTCAGTCGAAAATGTATATTTAAATAAAAACATAAACACATATAAAAATAAAAGCATAAATTGATATTTATGCTTTTATTG‑TTTTTGAGGTGATATAAAGCAACGGCTATATTACAAATTCACGCCCTTTTGCCA  >  1/1455082‑1455441                                                                                                                                                                     |                                                                                                                                                                                                     cccGGTGATAAGAAAACCGCTTGTGAAGCTCTACTCTGCTTATCATCACCAAAAGACGCGGCGAAAACGCCGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTtaaata                                                                                                                <  1:589770/250‑1 (MQ=255)                     ttgtgAAGCTCTACTCTGCTTATCATCACCAAAAGACGCGGCGAAAACGCCGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATaaatatacttaag                                                                                            <  1:62019/250‑11 (MQ=255)                                      gCTTATCATCACCAAAAGACGCGGCGAAAACGCCGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATaaatat                                                                                                   >  1:659119/1‑224 (MQ=255)                                          atcatcACAAAAAGACGCGGCAAAAACGCCGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATAAATATACTTAAGAATAAAAGCATAAAtatttat                                                                       <  1:856024/251‑7 (MQ=255)                                                                    cGCCGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATaaat                                                                                                     >  1:537831/1‑194 (MQ=255)                                                                      ccGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAGAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACAt                                                                                                         >  1:489645/1‑189 (MQ=255)                                                                             gTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATaaatatac                                                                                                 >  1:668120/1‑185 (MQ=255)                                                                                aTCCGCCGTTAAATAAGGTCTATTCAAATCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATAAATATACTTAAGAATAAAAGCATAAATATTTATTTATACTTTTATTGTTTTTTGAGGTGATATAACGCacc                                 <  1:357263/251‑3 (MQ=255)                                                                                  ccgccgTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATAAATATACTTAAGAATAAAAGCATAAAtatttatttatac                                                                 >  1:573286/1‑205 (MQ=255)                                                                                     ccgTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATaaatatacttaagaataaaag                                                                                    >  1:106465/1‑177 (MQ=255)                                                                                           aaTAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATAAATATACTTAAGAATAAAAGCATAAATATTTATTTATACTTTTATTGTTTTTTGAGGTGATATAACGCACCGGCTATATTa                       >  1:603365/1‑250 (MQ=255)                                                                                            aTAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATAAATATACTTAAGAATAAAAGCATAAATATTTATTTATACTTTTATTGTTTTTTGAGGTGATATaa                                       >  1:344237/1‑233 (MQ=255)                                                                                                     aTTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATAAATATACTTAAGAATAAAAGCATAAATATTTATTTATACTTTTATTGTTTTTTGAGGTGATATAACGCACCGGCTATATTACAAATTcccg             >  1:724935/1‑247 (MQ=255)                                                                                                               ggCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAATTTTTTAAAATTGTGTCCAAATTCTAACGAGGACTCTATTAATCAAGCATTATCGAGCAAATGCGCCCGCAAGCCGTTTAGTCGAAGATGTATGTTTAAATAAAAACATAAATATACTTAAGAATAAAAGCATAAATATTTATTTATACTTTTATTGTTTTTTGAGGTGATATAACGCACCGGCTATATTACAAATTCCCGCCCTTTTGcta  >  1:17085/1‑249 (MQ=255)                                                                                                                                                                     |                                                                                                                                                                                                     ACCGGTGATACCAAAACCGCCTGTGAAGCTCTACTCTGCTTATCATCACCAAAAGACGCGGCGAAAACGCCGGAATGTCATCCGCCGTTAAATAAGTTCTATTCAATTCGGGCTAAAAAACCACACGAAACATTAAAGAAACGTCGTAACTTTTTAAAACTATGCCCGGATGCTAGTGATGATTCAGTGAACCAGGCTTTATCTACAAAATGCGCCCGCAAGCCGTTCAGTCGAAAATGTATATTTAAATAAAAACATAAACACATATAAAAATAAAAGCATAAATTGATATTTATGCTTTTATTG‑TTTTTGAGGTGATATAAAGCAACGGCTATATTACAAATTCACGCCCTTTTGCCA  >  1/1455082‑1455441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑