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%	140.4 / NA	44	*111* (TAG→TAA)	ygiZ	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/24);  total (20/24)

ATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATA  >  NC_000913/3171744‑3172005                                                                                                                                        |                                                                                                                               aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTaa                                                                                                                             <  1:134348/139‑1 (MQ=255)       cAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                       <  2:191065/139‑1 (MQ=255)       cAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                       <  2:626600/139‑1 (MQ=255)                                            ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                                                                  >  1:1252534/1‑139 (MQ=255)                                            ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                                                                  >  2:574868/1‑139 (MQ=255)                                             tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGTTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAAGATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 <  2:759556/139‑1 (MQ=255)                                             tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 <  2:186591/139‑1 (MQ=255)                                             tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 <  2:1252534/139‑1 (MQ=255)                                              gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                                                                                >  1:1156257/1‑139 (MQ=255)                                              gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                                                                                >  2:1508849/1‑139 (MQ=255)                                                           aaGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                                   <  1:1508849/139‑1 (MQ=255)                                                           aaGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                                   <  2:1156257/139‑1 (MQ=255)                                                            aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                  <  1:574868/139‑1 (MQ=255)                                                             ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                 <  2:429843/139‑1 (MQ=255)                                                                       aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                           >  1:1179072/1‑103 (MQ=255)                                                                       aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                           <  2:1179072/103‑1 (MQ=255)                                                                             aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTaa                                                                                                                             >  2:795568/1‑63 (MQ=255)                                                                             aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTaa                                                                                                                             <  1:795568/63‑1 (MQ=255)                                                                                              ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 >  2:366681/1‑90 (MQ=255)                                                                                              ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 <  1:366681/90‑1 (MQ=255)                                                                                                           agCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtt                   >  2:575130/1‑139 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                           <  1:1221229/65‑1 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                           >  2:1221229/1‑65 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 >  1:57090/1‑75 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 <  2:57090/75‑1 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 >  1:1247621/1‑75 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 <  2:1247621/75‑1 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                 <  1:570331/91‑1 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                 >  2:570331/1‑91 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                 >  1:240917/1‑139 (MQ=255)                                                                                                             ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                 >  1:1289633/1‑139 (MQ=255)                                                                                                                ttAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt              <  1:264191/139‑1 (MQ=255)                                                                                                                ttAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt              >  2:264191/1‑139 (MQ=255)                                                                                                                 tAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                 <  2:78353/87‑1 (MQ=255)                                                                                                                 tAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                 >  1:78353/1‑87 (MQ=255)                                                                                                                      cGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata        <  2:240917/139‑1 (MQ=255)                                                                                                                       gACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataa       <  2:1289633/139‑1 (MQ=255)                                                                                                                        aCTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                                                                      >  1:1339402/1‑59 (MQ=255)                                                                                                                        aCTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                                                                      <  2:1339402/59‑1 (MQ=255)                                                                                                                          tctcATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 >  2:514769/1‑62 (MQ=255)                                                                                                                          tctcATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                 <  1:514769/62‑1 (MQ=255)                                                                                                                            tcATAGATTTTATTAATAGCAAAATATAACCCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaata  >  1:1026004/1‑139 (MQ=255)                                                                                                                              aTAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                                   >  1:1044362/1‑72 (MQ=255)                                                                                                                              aTAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                                   <  2:1044362/72‑1 (MQ=255)                                                                                                                                        |                                                                                                                               ATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATA  >  NC_000913/3171744‑3172005

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