Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	686,137	T→C	N178S (AAC→AGC)	gltJ ←	glutamate/aspartate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	686,137	0	T	C	100.0%	99.4 / NA	29	N178S (AAC→AGC)	gltJ	glutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/13);  total (16/13)

AGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGTTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATCGCG  >  NC_000913/685998‑686274                                                                                                                                            |                                                                                                                                           aGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAAAGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGGTTTTCACCAGGCTCATCAtctc                                                                                                                                  >  2:274124/1‑150 (MQ=255)        aTCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCAGAGTAATCGAGCAATTTACCCCCCTGCGCCGCCCTATCCACCAGACCGATAGTCGAGGCGATGGCGGAGGTTTTCACCAGGCTCATCATCTCTGAgg                                                                                                                             >  1:404442/1‑148 (MQ=255)            acacATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATcgg                                                                                                                      >  2:179010/1‑151 (MQ=255)               caTACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATcggcgg                                                                                                                   >  2:296684/1‑151 (MQ=255)                     ccAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCggg                                                                                                                  >  2:187255/1‑146 (MQ=255)                      cAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGGTTTTCACCAGGCTCATCATCTCTGAGGGCATCGGCgggg                                                                                                                 >  1:692340/1‑145 (MQ=255)                              gggggaaaacccccccccAGGCAGGGCCGAAATAATCGACCATTTCCCCCGCCGGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATa                                                                                                    <  1:756144/136‑1 (MQ=255)                                     aaaCGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGGGAACGCGGTaa                                                                                                   >  1:359142/1‑145 (MQ=255)                                         gACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAAAGCGATTAggaat                                                                                              >  2:474110/1‑142 (MQ=255)                                             cccAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGCCGAGAACGCGATAAGCAATA‑Ggaagc                                                                                      >  1:610079/1‑149 (MQ=255)                                                 ggCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGAGCGATAACCCGGTTAGCAATaa                                                                                            >  2:849351/1‑139 (MQ=255)                                                              tAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGGAAACGCGATAAGCATTA‑GGCAGCAGGACGTATCGATAAGcc                                                                    >  1:319321/1‑150 (MQ=255)                                                                  gaacaaattgccccTCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCg                                                                <  2:628846/139‑1 (MQ=255)                                                                            cccccccGGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCa                                                       <  2:311388/148‑1 (MQ=255)                                                                            ccccGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGt                                                         <  1:343503/146‑1 (MQ=255)                                                                             cccGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGGAGCAGGACGTATCGATAAGCCTGCGGTAGCGTCAgg                                                     >  2:417395/1‑150 (MQ=255)                                                                             cccGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCa                                                       <  2:273416/148‑1 (MQ=255)                                                                                gCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGGCGGAAAAGCGATAAGGAATAGGGCAG‑AGTACGGATAGATAAGccgtcgggagggtta                                                       >  2:598459/1‑132 (MQ=255)                                                                                                          gATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAAAGCGATAAGAATTA‑GGCAGCAGGACGTATCGATAAACCCGCGGGAGCGGCAGGCCCATCCCCAGCGCGGCATTTTTTTTccc                        >  2:262935/1‑150 (MQ=255)                                                                                                              gTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTcgcg                   <  2:439240/151‑1 (MQ=255)                                                                                                              gTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTcgcg                   <  1:233860/151‑1 (MQ=255)                                                                                                                 gAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCa                <  1:435605/151‑1 (MQ=255)                                                                                                                   ggCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGc              <  2:727731/151‑1 (MQ=255)                                                                                                                     cGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGGCGGTAACGCGCTAAACATTA‑GGCAGCAGTACGTATCGAAAAACCTGAGGGAGCGGCAGGCCCCTCGCCCGCGCGGGATTTTTTTTCCCTcgccg                  >  1:622228/1‑143 (MQ=255)                                                                                                                         ggCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCATGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGa        <  1:788083/151‑1 (MQ=255)                                                                                                                            ggAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATc     <  1:376903/151‑1 (MQ=255)                                                                                                                              aGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATcgc   <  1:105577/151‑1 (MQ=255)                                                                                                                               gTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAATGCGATAAGCATTA‑AGCAGCAGGACGGATCGATAAACCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTTCCCTCGCGGCAGCGACTGAATcgc   >  2:511232/1‑150 (MQ=255)                                                                                                                               gTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATcgcg  <  2:245034/151‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           AGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGTTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTA‑GGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATCGCG  >  NC_000913/685998‑686274

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