Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,553,095	C→A	P4Q (CCG→CAG)	malT →	DNA‑binding transcriptional activator MalT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,553,095	0	C	A	100.0%	66.0 / NA	19	P4Q (CCG→CAG)	malT	DNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/5);  total (14/5)

GATCAGCGTCGTTTTAGGTGAGTTGTTAATAAAGATTTGGAATTGTGACACAGTGCAAATTCAGACACATAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCCGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGCTGGATGAAGGTGATAACCAGCAAGAGCGTTTCGCCAGCTATCTCATTGCCGCCGTG  >  NC_000913/3552903‑3553330                                                                                                                                                                                                 |                                                                                                                                                                                                                                            gATCAGCGTCGGTTTAGGTGAGTTGTTTATCAAGCTTTGGAATTGTGACACAGTGCAAATTCAGACACATAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAACCTAAGACGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCctgctg                                                                                                                                                                                   >  2:183781/1‑251 (MQ=255)    cAGCGTCGTTTTAGGTGAGTTGTTAATAAAGATTTGGAATTGTGACACAGTGCAAATTCAGACACATAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGcc                                                                                                                                                                                >  2:96554/1‑250 (MQ=255)     aGCGTCGTTTTAGGTGAGTTGTTAATAAAGATTTGGAATTGTGACACAGTGCAAATTCAGACACATAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTa                                                                                                                                                                               <  1:208574/251‑1 (MQ=255)           gtTTTAGGTGAGTTGTTAATAAAGATTTGGAATTGTGACACAGTGCAAATTCAGACAAATAAAAAAACGTCATCGCTTGCATTAGAAAGGATTCTGGCCGACCTTATAACCATTAATTACGACGCGCAAAAAACATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACttt                                                                                                                                                                         >  1:41703/1‑251 (MQ=255)                                               gACACAGTGCAAATTCAGACACATAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGACCAACTTCCGGCTGGCGCTGATCACg                                                                                                                                      >  1:116245/1‑250 (MQ=255)                                                         aaaTTCAGACACATAAACAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCggg                                                                                                                           >  2:223622/1‑251 (MQ=255)                                                           aTTCAGACACATAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCt                                                                                                                         >  1:222959/1‑251 (MQ=255)                                                                  cacaTAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGaa                                                                                                                   >  1:198625/1‑250 (MQ=255)                                                                            aaCGTCATCGCTTGCATTAGAACGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAACAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTc                                                                                                        >  1:159051/1‑251 (MQ=255)                                                                                          cATTAGAAAGGTTTCTGGCCGACCTTATAACCCTTAATTACGAAGCGTAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGCAAGACCACCCTCATTTCCCAGTGggc                                                                                          >  1:15414/1‑251 (MQ=255)                                                                                                     tttCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCaaa                                                                                >  2:218699/1‑250 (MQ=255)                                                                                                            ccGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGatat                                                                        >  1:182532/1‑251 (MQ=255)                                                                                                                                  cGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGCTGGATGa                                                   >  1:81971/1‑250 (MQ=255)                                                                                                                                                                  tCCACAGTGAAGTGATTCCCTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATAACGTGGTTCGTGAGCGACTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGATGGATGAAGGCGATAACCAGCAAGAGCGGTTCGCCAGCt                   >  2:158220/1‑250 (MQ=255)                                                                                                                                                                          gaagtgaTTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGATTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGCTGGATGAAGGTGATAACCAGCAAGAGCGTTTCGCCAGCTATCTCAtt           <  1:124928/250‑1 (MQ=255)                                                                                                                                                                          gaagtgaTTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGCTGGATGAAGGTGATAACCAGCAAGAGCGTTTCGCCAGCTATCTCAtt           <  1:124934/250‑1 (MQ=255)                                                                                                                                                                            agtgaTTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGCTGGATGAAGGTGATAACCAGCAAGAGCGTTTCGCCAGCTATCTCATTgcc        <  1:232804/251‑1 (MQ=255)                                                                                                                                                                               gaTTAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGCTGGATGAAGGTGATAACCAGCAAGAGCGTTTCGCCAGCTATCTCATTgccgcc     <  2:198625/251‑1 (MQ=255)                                                                                                                                                                                  tAACTATGCTGATTCAGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCCGGCAAAAACGATGTCGGCTGGTACTCGCTGGATGAAGGTGATAACCAGCAAGAGCGTTTCGCCAGCTATCTCATTGCCGCCGTg  >  1:76414/1‑251 (MQ=255)                                                                                                                                                                                                 |                                                                                                                                                                                                                                            GATCAGCGTCGTTTTAGGTGAGTTGTTAATAAAGATTTGGAATTGTGACACAGTGCAAATTCAGACACATAAAAAAACGTCATCGCTTGCATTAGAAAGGTTTCTGGCCGACCTTATAACCATTAATTACGAAGCGCAAAAAAAATAATATTTCCTCATTTTCCACAGTGAAGTGATTAACTATGCTGATTCCGTCAAAACTAAGTCGTCCGGTTCGACTCGACCATACCGTGGTTCGTGAGCGCCTGCTGGCTAAACTTTCCGGCGCGAACAACTTCCGGCTGGCGCTGATCACGAGTCCTGCGGGCTACGGAAAGACCACCCTCATTTCCCAGTGGGCGGCAGGCAAAAACGATATCGGCTGGTACTCGCTGGATGAAGGTGATAACCAGCAAGAGCGTTTCGCCAGCTATCTCATTGCCGCCGTG  >  NC_000913/3552903‑3553330

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