Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,171,879	C→T	*111* (TAG→TAA)	ygiZ ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,171,879	0	C	T	100.0%	84.2 / NA	27	*111* (TAG→TAA)	ygiZ	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/21);  total (6/21)

TCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATA  >  NC_000913/3171748‑3171999                                                                                                                                    |                                                                                                                         tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                   <  1:108536/139‑1 (MQ=255)  ccAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                  <  1:804421/139‑1 (MQ=255)  ccAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                  <  1:451901/139‑1 (MQ=255)   cAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                 <  1:453230/139‑1 (MQ=255)   cAGAGAATCGCGAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                 <  1:697060/139‑1 (MQ=255)          tCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGA‑TTTATTAATAGCAAAATATaaa                                                                                                         <  1:361793/139‑1 (MQ=255)                   tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTcccc                                                                                                 <  1:777550/139‑1 (MQ=255)                         gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                           <  2:459424/139‑1 (MQ=255)                                         tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                  <  2:323884/100‑1 (MQ=255)                                         tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                  >  1:323884/1‑100 (MQ=255)                                         tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                           <  2:724437/139‑1 (MQ=255)                                               tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                                                                     >  2:760159/1‑139 (MQ=255)                                                       aaGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                             <  2:511098/139‑1 (MQ=255)                                                        aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                            <  1:760159/139‑1 (MQ=255)                                                        aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                            <  1:437873/139‑1 (MQ=255)                                                         ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                           <  2:425068/139‑1 (MQ=255)                                                              aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                           <  1:68759/134‑1 (MQ=255)                                                              aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                           >  2:68759/1‑134 (MQ=255)                                                                                             tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                          >  2:862561/1‑104 (MQ=255)                                                                                             tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                          <  1:862561/104‑1 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                     >  2:159082/1‑65 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                     <  1:159082/65‑1 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt           >  1:534518/1‑139 (MQ=255)                                                                                                                 aCGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat   <  1:691634/139‑1 (MQ=255)                                                                                                                 aCGACTCTCATAGACTTTATTAATAGCAAAATATAAACCGCCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat   <  2:14808/139‑1 (MQ=255)                                                                                                                  cGACTCTCATAGATTTTATTAATAGCAAAATAAAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata  <  2:281289/139‑1 (MQ=255)                                                                                                                  cGACTCTCACAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata  <  2:633648/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                         TCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATA  >  NC_000913/3171748‑3171999

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