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%	153.2 / 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)

TTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAAC  >  NC_000913/3131212‑3131473                                                                                                                                  |                                                                                                                                     ttAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAgg                                                                                                                             >  2:1042244/1‑139 (MQ=17)  tAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACAGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGt                                                                                                                            >  2:1229942/1‑139 (MQ=18)         cGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGc                                                                                                                     >  2:206121/1‑139 (MQ=255)          gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    <  1:380006/139‑1 (MQ=255)          gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    >  1:338674/1‑139 (MQ=255)          gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    <  2:24679/139‑1 (MQ=255)          gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    <  1:806617/139‑1 (MQ=255)          gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    >  2:806617/1‑139 (MQ=255)                      tAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCtt                                                                                                                   >  1:289921/1‑128 (MQ=255)                      tAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCtt                                                                                                                   <  2:289921/128‑1 (MQ=255)                             gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    <  1:1019802/120‑1 (MQ=255)                             gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    >  2:1019802/1‑120 (MQ=255)                             gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaa                                                                                                 >  1:1285392/1‑139 (MQ=255)                             gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACAGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATa                                                                                                                        <  2:73248/116‑1 (MQ=255)                             gCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACAGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATa                                                                                                                        >  1:73248/1‑116 (MQ=255)                                  ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                            >  2:269143/1‑139 (MQ=255)                                  ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                            <  1:1229942/139‑1 (MQ=255)                                  ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCACACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                                                            >  2:1422784/1‑139 (MQ=255)                                            ttACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTaa                                                                                                       >  2:1283980/1‑118 (MQ=255)                                            ttACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTaa                                                                                                       <  1:1283980/118‑1 (MQ=255)                                             tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    >  2:479534/1‑104 (MQ=255)                                             tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCt                                                                                                                    <  1:479534/104‑1 (MQ=255)                                                 ggAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGc                                                                             <  1:378524/139‑1 (MQ=255)                                                  gagaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTGAACGCCAAACACCAGCGAAATGAGCGGCt                                                                            >  2:1090496/1‑139 (MQ=255)                                                     aaaaaaTCGGCTCAAACATGAAGAAATGAAATGCCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCAGCTTCCTTAAATCAGGTCATACGCTTCgag                                                                                                               <  1:942268/101‑1 (MQ=255)                                                     aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCgag                                                                                                               >  2:942268/1‑101 (MQ=255)                                                       aaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                                       <  2:832829/139‑1 (MQ=255)                                                           tCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcgacga                                                                   <  2:1285392/139‑1 (MQ=255)                                                                           gaaatgaaatGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAgcgc                                                   <  2:1079173/139‑1 (MQ=255)                                                                           gaaatgaaatGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCACACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAgcgc                                                   <  1:1422784/139‑1 (MQ=255)                                                                            aaatgaaatgACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGcc                                                  <  1:206121/139‑1 (MQ=255)                                                                             aatgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGcca                                                 <  2:1139020/139‑1 (MQ=255)                                                                                      aCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATAc                                        >  1:1297455/1‑139 (MQ=255)                                                                                       cTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATAcc                                       <  2:1297455/139‑1 (MQ=255)                                                                                         gAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGa                                     <  1:269143/139‑1 (MQ=255)                                                                                                ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCa                                                                                                                          <  1:299330/47‑1 (MQ=255)                                                                                                ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCa                                                                                                                          >  2:299330/1‑47 (MQ=255)                                                                                                ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACCCCGCGAGCGCCACCAGAATACCGATTccgcc                              <  1:1042244/139‑1 (MQ=255)                                                                                                             ccccGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACCCCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCg                 <  1:738684/139‑1 (MQ=255)                                                                                                             ccccGCTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCg                 <  2:1013788/139‑1 (MQ=255)                                                                                                                  cTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                                       >  1:547806/1‑80 (MQ=255)                                                                                                                  cTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCcg                                                                       <  2:547806/80‑1 (MQ=255)                                                                                                                  cTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccac                                                <  2:1327702/103‑1 (MQ=255)                                                                                                                  cTTATTCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccac                                                >  1:1327702/1‑103 (MQ=255)                                                                                                                      ttCGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACCCCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTcaca        <  1:1422261/139‑1 (MQ=255)                                                                                                                        cGCACCTTCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACAcc      >  2:113477/1‑139 (MQ=255)                                                                                                                            ccttccttAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTaac  <  1:113477/139‑1 (MQ=255)                                                                                                                               tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGa                                     >  2:1096508/1‑101 (MQ=255)                                                                                                                               tCCTTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGa                                     <  1:1096508/101‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     TTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAAC  >  NC_000913/3131212‑3131473

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