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

TGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAG  >  NC_000913/3171746‑3172010                                                                                                                                      |                                                                                                                                    tGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACCAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATCGATTTTATTAGTa                                                                                                                                >  2:890824/1‑139 (MQ=255) tGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATa                                                                                                                                >  2:1019590/1‑139 (MQ=255)   tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                              <  2:515844/139‑1 (MQ=255)             cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATaaa                                                                                                                    >  1:228632/1‑139 (MQ=255)                     tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTcccc                                                                                                            >  1:1139537/1‑139 (MQ=255)                           gcatgAATGTTGAATTTGGGAATGCGGAGCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAATCGACTCTCATAGATTTTATTAATAGCa                                                                                                                             >  1:853960/1‑116 (MQ=255)                           gcatgAATGTTGAATTTGGGAATGCGGAGCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                             <  2:853960/116‑1 (MQ=255)                              tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                                                                                   >  1:24595/1‑139 (MQ=255)                                 aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                <  2:866723/139‑1 (MQ=255)                                 aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                <  2:228632/139‑1 (MQ=255)                                           tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                      <  2:1139537/139‑1 (MQ=255)                                           tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                      <  1:1019590/139‑1 (MQ=255)                                           tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                      <  2:1344173/139‑1 (MQ=255)                                                          aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                       <  2:24595/139‑1 (MQ=255)                                                                          aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                       <  1:1157935/123‑1 (MQ=255)                                                                          aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                       >  2:1157935/1‑123 (MQ=255)                                                                          aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc                                                       >  2:682233/1‑139 (MQ=255)                                                                            ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCa                                                     >  2:1026454/1‑139 (MQ=255)                                                                            ggTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATAAACCCa                                                     >  2:24920/1‑139 (MQ=255)                                                                                     gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGC‑AAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                       <  2:535452/111‑1 (MQ=255)                                                                                     gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGC‑AAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                       >  1:535452/1‑111 (MQ=255)                                                                                       aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                             <  2:119533/56‑1 (MQ=255)                                                                                       aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                             >  1:119533/1‑56 (MQ=255)                                                                                       aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                       >  2:326118/1‑110 (MQ=255)                                                                                       aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                       <  1:326118/110‑1 (MQ=255)                                                                                            ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                                     >  1:300290/1‑139 (MQ=255)                                                                                               tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                                    >  2:1090579/1‑89 (MQ=255)                                                                                               tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                                    <  1:1090579/89‑1 (MQ=255)                                                                                               tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGc                                  >  1:1112051/1‑139 (MQ=255)                                                                                               tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATGAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAACGAGCGATCATGCATATAACCAAAACTCGataatatacctc                                  >  2:1284477/1‑128 (MQ=255)                                                                                                      gagaGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                                   <  2:418804/99‑1 (MQ=255)                                                                                                      gagaGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                                   >  1:418804/1‑99 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                            >  1:864566/1‑37 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                            <  2:864566/37‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                             >  1:1262317/1‑36 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                             <  2:1262317/36‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                                                                     >  1:798914/1‑60 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                                                                     <  2:798914/60‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                <  1:1134307/65‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                >  2:1134307/1‑65 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                                   >  2:1273425/1‑94 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                                   <  1:1273425/94‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                                                                         <  2:299021/88‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                                                                         >  1:299021/1‑88 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                                                                         >  1:622041/1‑88 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcg                                                                         <  2:622041/88‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                                 <  2:622527/96‑1 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                                 >  1:622527/1‑96 (MQ=255)                                                                                                           ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                      >  1:820636/1‑139 (MQ=255)                                                                                                                   aCGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat              <  2:300290/139‑1 (MQ=255)                                                                                                                   aCGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat              <  1:682233/139‑1 (MQ=255)                                                                                                                      aCTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                                   <  1:843173/83‑1 (MQ=255)                                                                                                                      aCTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                                   >  2:843173/1‑83 (MQ=255)                                                                                                                        tctcATAGATTTTATTAAGAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATGAGaa         >  1:887005/1‑139 (MQ=255)                                                                                                                               gATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAg  <  2:1112051/139‑1 (MQ=255)                                                                                                                                  tttATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat              >  2:505674/1‑124 (MQ=255)                                                                                                                                  tttATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat              <  1:505674/124‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    TGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAG  >  NC_000913/3171746‑3172010

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