Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,443,961	G→A	A46V (GCC→GTC)	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,961	0	G	A	100.0%	52.9 / NA	19	A46V (GCC→GTC)	secY	Sec translocon subunit SecY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/7);  total (12/7)

CCCTTCTTTCTTAATTTCTGCCAACGTTGGGTGAACCACCGTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGGCATCAATACCAGGGATCGGAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTTAGTAACACGCAGGCCACGAACAGTTACCGGAGTCGTTACTTCG  >  NC_000913/3443771‑3444198                                                                                                                                                                                               |                                                                                                                                                                                                                                               cccTTCTTTCTTAATTTCTGCCAACGTTGGGTGAACCACCGTCAGCAGCTGGATAATGATAGACGCCGAAATATACGGCATGATCCCCTGAGCAAAGATAGAAGCACGGCTGAGAACACCACCAGAGAACATGTTAAAAATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAACTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCAACGATa                                                                                                                                                                                    >  1:315638/1‑250 (MQ=255)                                  aaCCACCGTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCGAAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCGCGCCTAAGc                                                                                                                                                  <  1:136369/251‑1 (MQ=255)                                    ncaccGTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATCg                                                                                                                                                                                                                              >  1:109957/2‑174 (MQ=255)                                         gTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATCAATACCAGGGATcn                                                                                                                                                                                                                              <  1:117520/169‑2 (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)                                                                                                     naaGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGACATc                                                                                                                                                                                                                                            >  1:93601/2‑95 (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)                                                                                                                                                                                               |                                                                                                                                                                                                                                               CCCTTCTTTCTTAATTTCTGCCAACGTTGGGTGAACCACCGTCAGCAGCTGGATAATGATCGACGCCGAAATATACGGCATGATCCCCAGAGCAAAGATAGAAGCACGGCTGAGAGCACCACCAGAGAACATGTTAAACATCTCAATGATGGTGCCTCGCTGTTGCTCAAGCAGTTTGGCAAGTACAGCGGCATCAATACCAGGGATCGGAATAAAAGA‑GCCAATACGGAACACAATCAGCGCACCGATAACAAACAGCAGTCTGCGTTTCAGCTCGCCTAAGCCACCTTTGGCACTTTGAAAATCTAATCCCGGTTGTTTAGCCATCTGCTACTTATTCCTCGATTTTACCGCCAGCAGCTTCGATAGCAGCACGAGCGCCTTTAGTAACACGCAGGCCACGAACAGTTACCGGAGTCGTTACTTCG  >  NC_000913/3443771‑3444198

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: ‑