Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,131,341	C→T	*356* (TAG→TAA)	yghQ ←	putative inner membrane polysaccharide flippase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,131,341	0	C	T	100.0%	78.7 / NA	26	*356* (TAG→TAA)	yghQ	putative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/18);  total (8/18)

TGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAA  >  NC_000913/3131217‑3131472                                                                                                                             |                                                                                                                                    tGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCAGCTTCCTTAAATCAGGTCATa                                                                                                                       <  1:715311/139‑1 (MQ=255)  ggCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATAc                                                                                                                      <  2:1041518/139‑1 (MQ=37)                 tAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTa                                                                                                       <  2:426443/139‑1 (MQ=255)                        gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                   >  1:101661/1‑120 (MQ=255)                        gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                   <  2:101661/120‑1 (MQ=255)                        gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaa                                                                                                >  2:315425/1‑139 (MQ=255)                         cTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaaa                                                                                               >  1:338857/1‑139 (MQ=255)                                   gCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCCAGTAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCACGTCATACGCTTCGAGATACTTAACGCCAAACACAAGCGaa                                                                                     >  2:359540/1‑139 (MQ=18)                                        tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgag                                                                                                              <  2:163596/109‑1 (MQ=255)                                        tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgag                                                                                                              >  1:163596/1‑109 (MQ=255)                                        tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAg                                                                                <  1:291075/139‑1 (MQ=255)                                                aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGTGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGc                                                                        <  2:673036/139‑1 (MQ=255)                                                         gCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                                      <  2:502141/132‑1 (MQ=255)                                                         gCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                                      >  1:502141/1‑132 (MQ=255)                                                                         atgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccac                                               <  2:399700/139‑1 (MQ=255)                                                                                           ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgcc                             <  1:939249/139‑1 (MQ=255)                                                                                           ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgcc                             <  1:475242/139‑1 (MQ=255)                                                                                            cGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCGTCGAGATACTTAACGCCAAACACTAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATAccgccc                            >  1:16532/1‑139 (MQ=255)                                                                                                aagaaTTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCa                        <  1:315425/139‑1 (MQ=255)                                                                                                   aaTTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTa                     <  1:449044/139‑1 (MQ=255)                                                                                                        ccccGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACCCCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCg                <  2:680171/139‑1 (MQ=255)                                                                                                             cTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACtt           >  1:867145/1‑139 (MQ=255)                                                                                                                      aCCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACCCCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTaa  <  2:867145/139‑1 (MQ=255)                                                                                                                          tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAAcc                 <  1:1055345/120‑1 (MQ=255)                                                                                                                          tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAAc                  <  1:724185/119‑1 (MQ=255)                                                                                                                          tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAAc                  <  1:646028/119‑1 (MQ=255)                                                                                                                             |                                                                                                                                    TGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAA  >  NC_000913/3131217‑3131472

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