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%	152.0 / NA	48	*356* (TAG→TAA)	yghQ	putative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/27);  total (21/27)

CTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCA  >  NC_000913/3131211‑3131465                                                                                                                                   |                                                                                                                             cTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAg                                                                                                                      >  1:858371/1‑139 (MQ=14)       tGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATa                                                                                                                <  2:1552343/139‑1 (MQ=25)         gCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACg                                                                                                              >  1:38120/1‑139 (MQ=38)         gCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACg                                                                                                              >  2:549876/1‑139 (MQ=38)          cGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCCCAAACATGAAGAAATGAAATGAATGAGTCAGCCGAGAAGAATTTCCCCGCTTATTAGCACCGTCCTTACATCAGGTCATTCGc                                                                                                             >  2:1232210/1‑139 (MQ=11)           gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                            <  2:272791/139‑1 (MQ=255)                 aCGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgaga                                                                                                      <  2:858371/139‑1 (MQ=255)                    gTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATAc                                                                                                   >  2:175061/1‑139 (MQ=255)                     tCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACt                                                                                                  <  1:175061/139‑1 (MQ=255)                        aaaTAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTaa                                                                                               <  2:38120/139‑1 (MQ=255)                         aaTAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAAc                                                                                              >  2:1481340/1‑139 (MQ=255)                               cTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaaa                                                                                        <  1:93820/139‑1 (MQ=255)                                 gATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAaca                                                                                      <  1:1481340/139‑1 (MQ=255)                                   ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                    >  1:1506660/1‑139 (MQ=255)                                            tttACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATa                                                                                                                >  2:772557/1‑102 (MQ=255)                                            tttACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATa                                                                                                                <  1:772557/102‑1 (MQ=255)                                              tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAg                                                                         <  1:1109306/139‑1 (MQ=255)                                                      aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGc                                                                 <  2:1675727/139‑1 (MQ=255)                                                      aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGc                                                                 >  1:1315466/1‑139 (MQ=255)                                                        aaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTGTTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                               <  2:1214457/139‑1 (MQ=255)                                                        aaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                               <  2:1315466/139‑1 (MQ=255)                                                          aaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcgac                                                             <  1:1232210/139‑1 (MQ=255)                                                                cTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAg                                                                         >  1:888808/1‑121 (MQ=255)                                                                cTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAg                                                                         <  2:888808/121‑1 (MQ=255)                                                                         gaagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCtt                                                                                                           >  2:679440/1‑78 (MQ=255)                                                                         gaagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCtt                                                                                                           <  1:679440/78‑1 (MQ=255)                                                                          aagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTaa                                                                                               >  2:237897/1‑89 (MQ=255)                                                                          aagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTaa                                                                                               <  1:237897/89‑1 (MQ=255)                                                                          aagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACcg                                                  >  2:1573957/1‑134 (MQ=255)                                                                          aagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACcg                                                  <  1:1573957/134‑1 (MQ=255)                                                                               atgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccac                                        >  2:954465/1‑139 (MQ=255)                                                                                       aCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgaga                                                                                                      >  1:376682/1‑69 (MQ=255)                                                                                       aCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgaga                                                                                                      <  2:376682/69‑1 (MQ=255)                                                                                           aGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAg                                              >  1:1018009/1‑121 (MQ=255)                                                                                           aGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAg                                              <  2:1018009/121‑1 (MQ=255)                                                                                           aGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGa                                    >  2:1544243/1‑131 (MQ=255)                                                                                           aGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGa                                    <  1:1544243/131‑1 (MQ=255)                                                                                               aGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccg                        <  1:549876/139‑1 (MQ=255)                                                                                                 ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaaa                                                                                        <  1:983479/73‑1 (MQ=255)                                                                                                 ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACCCCGCGAGCGCCACCAGAATACCGATTccgcc                      <  2:1506660/139‑1 (MQ=255)                                                                                                 ccGAGAAGAATGTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaaa                                                                                        >  2:983479/1‑73 (MQ=255)                                                                                                           tttCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGcca                                         <  1:1473761/110‑1 (MQ=255)                                                                                                           tttCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGcca                                         >  2:1473761/1‑110 (MQ=255)                                                                                                                ccGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGa       >  1:567805/1‑139 (MQ=255)                                                                                                                ccGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGa       <  1:954465/139‑1 (MQ=255)                                                                                                                     tATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTca  <  2:567805/139‑1 (MQ=255)                                                                                                                                tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc                         >  2:865829/1‑105 (MQ=255)                                                                                                                                tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc                         <  1:865829/105‑1 (MQ=255)                                                                                                                                   |                                                                                                                             CTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCA  >  NC_000913/3131211‑3131465

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