Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,443,980	G→A	P40S (CCG→TCG)	secY ←	Sec translocon subunit SecY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,980	0	G	A	100.0%	54.0 / NA	20	P40S (CCG→TCG)	secY	Sec translocon subunit SecY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/8);  total (12/8)

CCCTTCTTTCTTAATTTCTGCCAACGTTGGGTGAACCACCGTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGGCATCAATACCAGGGATCGGAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTTAGTAACACGCAGGCCACGAACAGTTACCGGAGTCGTTACTTCGCCAGCCAGGATCACTTTCGCG  >  NC_000913/3443771‑3444219                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 cccTTCTTTCTTAATTTCTGCCAACGTTGGGTGAACCACCGTCAGCAGCTGGATAATGATAGACGCCGAAATATACGGCATGATCCCCTGAGCAAAGATAGAAGCACGGCTGAGAACACCACCAGAGAACATGTTAAAAATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAACTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCAACGATa                                                                                                                                                                                                         >  1:315638/1‑250 (MQ=255)                                  aaCCACCGTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCGCGCCTAAGc                                                                                                                                                                       <  1:136369/251‑1 (MQ=255)                                             gcagcTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGc                                                                                                                                                             >  1:71502/1‑250 (MQ=255)                                             gcagcTGGATAATGATCGACGCCGAAATATACGGCATGATACCCAGAGCAACGCTAGCAGCACGGCTGAGAGCACTACCAGAGACCATGTTAAACATCTCACTGATGGTGCCTCGCTGTTGCTCAAGCCGTTTGGCAAGTACAGCGACCTCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAac                                                                                                                                                                                                                           >  1:294094/1‑188 (MQ=255)                                                                                     ccccAGAGCAAAGATAGAAGCACGGCTGGGATGACCAGCAGAGAACATGTTAATCATCTCAATGATGGTGCCTCGCTGTTTCTCAAGCAGGTTGGCAAGTAGAGCTACAGCGATACCAAGGGTCGAAGGGAAATATGACG‑TTCGGAGCTCAATCAGCGCTACGATAACAAACAGCCATCTGCGTTTCAGCTCGCCTAAGCCAACTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACt                                                                                                                    <  1:174232/251‑1 (MQ=255)                                                                                      nccAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACCTCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTc                                                                                                                                                                                            >  1:93597/2‑178 (MQ=255)                                                                                          nagCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAac                                                                                                                                                                                                                           >  1:52493/2‑143 (MQ=255)                                                                                                   tAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGc                                                                                                                        >  1:17651/1‑233 (MQ=255)                                                                                                            ggCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCACTACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCcagc                                                                                             >  1:185015/1‑251 (MQ=255)                                                                                                                        caccaGAGAACATGTTAAACATCTCANTGATGGNNCCNCGCTGTTGCTCAAGCAGTTNGGNNNGNNNANCGACANCANNACCAGGGATNGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGNCTGCGTTTCNNCTCGCCTAAGCCACCTTTGGCNCTTNNANNANCTANTc                                                                                                                                            >  2:98925/1‑192 (MQ=255)                                                                                                                        caccaGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCTCCTTTGGCACTTTGAAAATCTAATccc                                                                                                                                          <  1:98925/194‑1 (MQ=255)                                                                                                                                              ntcAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTTAGTAac                                                            >  1:268434/2‑250 (MQ=255)                                                                                                                                                                                  ggCAAGTACAGNGACATCAATACCAGGGATCGAAATAAAANA‑GCCAATACGGNACACANNNannn                                                                                                                                                                                                                 >  2:59153/1‑62 (MQ=255)                                                                                                                                                                                  ggCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAgcg                                                                                                                                                                                                                 <  1:59153/65‑1 (MQ=255)                                                                                                                                                                                   gccgatctaGCGACATCAATACCAGGGATCGAAAGAAAAGA‑GCCAATACGGAACACAATCAGTGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCGTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTTAGTAACACGCAGGCCACGAACAGTTACCGGAGTCGTTACTTCg                       <  1:147001/242‑1 (MQ=255)                                                                                                                                                                                    cAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCATCGCACCGATAGCAAACAGCAGTCTTCGTTTc                                                                                                                                                                                    <  1:372733/92‑1 (MQ=255)                                                                                                                                                                                             cgacATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGTACGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGc                                                                                                                               <  1:191871/133‑1 (MQ=255)                                                                                                                                                                                             cgacANCAATANCNNCGCTCGAAATTAAAGA‑GCCNATACGGAACACAATCNGCGCNCCGATnacn                                                                                                                                                                                                      >  2:191871/4‑64 (MQ=255)                                                                                                                                                                                                cATCAATACCAGGGATCGAAATAAAAGA‑GCAAATACGGAACACAATCAGCGCACCGAAAACAAACAGCAGTCTGCGGTTCAGCCCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAAGCCCGGTTGTTTAGCCATCTGCTACTTAGTCCTCAATTGTTCCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTGAGTAACACGCAGGCCACGAACAGTTACCgg                                    <  1:189139/224‑1 (MQ=255)                                                                                                                                                                                                        ccAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTTAGTAACACGCAGGCCACGAACAGTTACCGGAGTCGTTACTTCGCCAGCCAGGATCACTTTcgcg  >  1:12716/1‑250 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 CCCTTCTTTCTTAATTTCTGCCAACGTTGGGTGAACCACCGTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGGCATCAATACCAGGGATCGGAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTTAGTAACACGCAGGCCACGAACAGTTACCGGAGTCGTTACTTCGCCAGCCAGGATCACTTTCGCG  >  NC_000913/3443771‑3444219

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