Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,645,530	T→C	E35G (GAG→GGG)	relE ←	Qin prophage; mRNA interferase toxin RelE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,645,530	0	T	C	100.0%	108.8 / NA	31	E35G (GAG→GGG)	relE	Qin prophage; mRNA interferase toxin RelE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/16);  total (15/16)

ACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACTCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGTTCATCCAGCGTCACA  >  NC_000913/1645401‑1645654                                                                                                                                  |                                                                                                                             aCAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTc                                                                                                                    >  2:170598/1‑140 (MQ=255)   aGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCaa                                                                                                                           >  2:520122/1‑131 (MQ=255)   aGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCaa                                                                                                                           <  1:520122/131‑1 (MQ=255)         cACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGct                                                                                                            <  2:315299/140‑1 (MQ=255)                   aCAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCaa                                                                                                  <  2:458256/140‑1 (MQ=255)                   aCAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCaa                                                                                                  <  2:184860/140‑1 (MQ=255)                    cAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAAc                                                                                                 <  1:170598/140‑1 (MQ=255)                    cAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAAc                                                                                                 <  2:462992/140‑1 (MQ=255)                                   taACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTc                                                                                  >  2:2944/1‑140 (MQ=255)                                         ggTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGccc                                                                            >  1:530942/1‑140 (MQ=255)                                         ggTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGccc                                                                            >  2:234870/1‑140 (MQ=255)                                          gTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCa                                                                           >  1:332029/1‑140 (MQ=255)                                            ataCAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGc                                                                         <  2:60328/140‑1 (MQ=255)                                                 aGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCg                                                                    >  1:122313/1‑140 (MQ=255)                                                    cGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCa                                                                 <  1:2944/140‑1 (MQ=255)                                                     gATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCAt                                                                <  1:409561/140‑1 (MQ=255)                                                       tAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTc                                                              >  1:342166/1‑140 (MQ=255)                                                                   cGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCg                                                  <  1:234870/140‑1 (MQ=255)                                                                    ggAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGc                                                 >  2:432182/1‑140 (MQ=255)                                                                     gAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCt                                                >  1:209290/1‑140 (MQ=255)                                                                          tAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCa                                           <  2:332029/140‑1 (MQ=255)                                                                           aaTCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCaa                                          >  2:439654/1‑140 (MQ=255)                                                                              cTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAAt                                       <  1:418049/140‑1 (MQ=255)                                                                                      aaTCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGaaaa                               <  2:122313/140‑1 (MQ=255)                                                                                               tACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATc                      <  2:530942/140‑1 (MQ=255)                                                                                                     ggAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGtt                >  1:472261/1‑140 (MQ=255)                                                                                                      gAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGTTc               >  2:16865/1‑140 (MQ=255)                                                                                                      gAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGTTc               >  2:356243/1‑140 (MQ=255)                                                                                                      gAGCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGTTc               >  1:190799/1‑140 (MQ=255)                                                                                                        gCTTGTTTGCTTCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGTTCAt             <  2:472261/140‑1 (MQ=255)                                                                                                                   tCAATCCGGGGTGACCCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGTTCATCCAGCGTcaca  <  2:477007/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                             ACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAATCTTGTAACAATCAGGCATACCACGGAGCTTGTTTGCTTCAATCCGGGGTGACTCAAGTACTTCAACCAGCTTCTTTTTCAACTGTTCACGTACCGTCGAGCCCAGCTTTCGCCATTCCTTTAGTGCCCGCTCGTCAAAATCCAGAAAATACGCCATCAGAGTTCATCCAGCGTCACA  >  NC_000913/1645401‑1645654

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