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

TCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGG  >  NC_000913/3171748‑3172013                                                                                                                                    |                                                                                                                                       tCCATCGAATCGGCAAATTGCATTGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                                 <  1:359565/139‑1 (MQ=255)  ccAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                                <  2:628696/139‑1 (MQ=255)   cAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                               <  2:647521/139‑1 (MQ=255)           cGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATaaa                                                                                                                       >  1:715582/1‑139 (MQ=255)                  ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTccc                                                                                                                >  2:677909/1‑139 (MQ=255)                       tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCaaaa                                                                                                           >  1:857779/1‑139 (MQ=255)                               aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                   <  1:677909/139‑1 (MQ=255)                               aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                   <  2:1190986/139‑1 (MQ=255)                                  ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                                                                                >  1:1080920/1‑139 (MQ=255)                                   tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                                                                               <  2:1080920/139‑1 (MQ=255)                                      aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                                                                            <  2:577307/139‑1 (MQ=255)                                         tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                         <  1:1125909/139‑1 (MQ=255)                                          gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                                                                                        >  1:100060/1‑139 (MQ=255)                                           ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                                 <  2:1012352/97‑1 (MQ=255)                                           ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                                 >  1:1012352/1‑97 (MQ=255)                                                gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                                <  2:880294/93‑1 (MQ=255)                                                gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                                >  1:880294/1‑93 (MQ=255)                                                       aaGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTGTTAATAGCaa                                                                                                                               >  1:1056093/1‑87 (MQ=255)                                                       aaGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCaa                                                                                                                               <  2:1056093/87‑1 (MQ=255)                                                         ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                         >  2:1374094/1‑139 (MQ=255)                                                           tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                                                                            >  2:463702/1‑118 (MQ=255)                                                           tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                                                                            <  1:463702/118‑1 (MQ=255)                                                             gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                                                                     <  2:100060/139‑1 (MQ=255)                                                                ttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa                                                                  <  2:715582/139‑1 (MQ=255)                                                                  gAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                   >  2:292330/1‑104 (MQ=255)                                                                  gAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                   <  1:292330/104‑1 (MQ=255)                                                                     gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                                             <  2:965071/107‑1 (MQ=255)                                                                     gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                                             >  1:965071/1‑107 (MQ=255)                                                                       gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                                    <  1:814803/130‑1 (MQ=255)                                                                       gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                                    >  2:814803/1‑130 (MQ=255)                                                                                   gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                               >  2:778498/1‑139 (MQ=255)                                                                                    tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                                 <  2:1128704/56‑1 (MQ=255)                                                                                    tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                                 >  1:1128704/1‑56 (MQ=255)                                                                                             tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGc                                     >  2:225988/1‑139 (MQ=255)                                                                                                  aTGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                         >  2:746634/1‑98 (MQ=255)                                                                                                  aTGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                         <  1:746634/98‑1 (MQ=255)                                                                                                  aTGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                                   <  1:988235/120‑1 (MQ=255)                                                                                                  aTGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                                   >  2:988235/1‑120 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                                >  2:56637/1‑36 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                                <  1:56637/36‑1 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                         <  2:553857/75‑1 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                         >  1:553857/1‑75 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                          >  2:967380/1‑90 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                          <  1:967380/90‑1 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                         >  2:399022/1‑139 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                         >  2:1168396/1‑139 (MQ=255)                                                                                                         ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACAGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                         >  2:678548/1‑139 (MQ=255)                                                                                                           gTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                       <  1:225988/139‑1 (MQ=255)                                                                                                            ttAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                                                                            >  2:973978/1‑69 (MQ=255)                                                                                                            ttAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                                                                            <  1:973978/69‑1 (MQ=255)                                                                                                             tAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                          <  1:149336/86‑1 (MQ=255)                                                                                                             tAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                          >  2:149336/1‑86 (MQ=255)                                                                                                                 aCGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat                 <  1:778498/139‑1 (MQ=255)                                                                                                                      tctcATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaa            >  1:1192700/1‑139 (MQ=255)                                                                                                                                tttATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                          >  1:1186815/1‑67 (MQ=255)                                                                                                                                tttATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                          <  2:1186815/67‑1 (MQ=255)                                                                                                                                tttATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataa               >  1:1474968/1‑126 (MQ=255)                                                                                                                                tttATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataa               <  2:1474968/126‑1 (MQ=255)                                                                                                                                tttATTAATAGCAAAATATAAACCGTCACCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTgg  >  2:998815/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       TCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGG  >  NC_000913/3171748‑3172013

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