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%	113.7 / NA	34	N178S (AAC→AGC)	gltJ	glutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/18);  total (16/18)

TGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGTTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATCGCGGCGCGTACC  >  NC_000913/685991‑686283                                                                                                                                                   |                                                                                                                                                   tGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGGAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGGTGGCGGAGGTTTTCACCCGGCtc                                                                                                                                                  >  1:662181/1‑149 (MQ=255)                    cacaTACGCCAGCGTAATGGCGGAAAACGCGTCCCGGGCACTGGCTGTGGGATCGTGCAAGTTAACCCCCGGCGCCGTCTTATTCACCAGACCGATAAGCCCGGGGAAGGCGGGGGTTTTGACCAGGCTCATCATTTTTGAggga                                                                                                                                   >  2:390675/1‑143 (MQ=255)                       aTACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGGTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGGGGGGTTTTTCACCAGGCTCATCATCTCTGAGGGCATCGGCggg                                                                                                                          >  2:301942/1‑151 (MQ=255)                           gCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCCTATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCggg                                                                                                                          >  1:271968/1‑147 (MQ=255)                              aGCGTAATGGCGGTAAACGACCCCCCGGCATGGGCTGAGTAATCGGGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGGTGGCGGGGTTTTTCACCCGGCTCATCATCTCTTAGGTCATCCGCggggagac                                                                                                                     >  1:176261/1‑144 (MQ=255)                                           tAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGCTAACACGATAAGCAtt                                                                                                      >  2:61347/1‑151 (MQ=255)                                                       aGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGAGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAAAGCGATAAGCATTAGGCAGCAGTa                                                                                           >  2:601607/1‑150 (MQ=255)                                                                   aGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAAAGCGATAAACATTAGGCAGCAGTACGTATCGATaa                                                                                >  2:100039/1‑149 (MQ=255)                                                                      aagggaaaatttcccccGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATaa                                                                                <  1:51076/132‑1 (MQ=255)                                                                       accgacattttcccccGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGc                                                                          <  1:753069/138‑1 (MQ=255)                                                                        tCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAAAGCGGTAAACAATTGGCAGCAGTACGTAACGCAAAACCTGc                                                                          >  1:20592/1‑150 (MQ=255)                                                                          gAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGGCGGTAACGCGGTAAGCCTTAGGCAGCAGAACGTATCGATAAGCCTGCgg                                                                        >  2:529390/1‑150 (MQ=255)                                                                           agcatttcccccGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCg                                                                         <  1:40049/139‑1 (MQ=255)                                                                           agaattttCCCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGc                                                                          <  1:373209/142‑1 (MQ=255)                                                                               aTTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGGTAACGCGATAAGCAATAGGGAGGAGGACGTATCGATAAAACTGCGGGAGCg                                                                   >  1:557366/1‑150 (MQ=255)                                                                                tcccccccccTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGc                                                                    <  1:142435/144‑1 (MQ=255)                                                                                    cccGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGc                                                             <  2:485485/151‑1 (MQ=255)                                                                                      cGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGccc                                                           <  2:55725/151‑1 (MQ=255)                                                                                            cgccgcCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCg                                                       <  1:630090/149‑1 (MQ=255)                                                                                                   atatCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCgg                                              <  2:677391/151‑1 (MQ=255)                                                                                                           caGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCAtttttt                                      <  2:288470/151‑1 (MQ=255)                                                                                                            aGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCAttttttt                                     <  2:168336/151‑1 (MQ=255)                                                                                                                   tAGTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGAGGGTAAAGCGATAAGCATTATGCAGCAGTACGGATCGATAAGCCTGCGGGAGCGTCAGGCCCATCCCCCGCGCGGCAtttttttttcccc                               >  1:631822/1‑144 (MQ=255)                                                                                                                     gTCGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTcgcg                            <  2:44681/151‑1 (MQ=255)                                                                                                                       cGAGGCGATGGCGGAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCgg                           <  1:356267/150‑1 (MQ=255)                                                                                                                         aGGCGATGGCGGAGCTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGGCGATAACGCGATAAGCATTAAGCAGCAGTAAGTATCGATAAGCCTGCGGGAGCGTTAGGGCCCTCGCCCGCGCGGGCTTTTTTTTccccc                               >  1:603843/1‑142 (MQ=255)                                                                                                                                   ggAGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATc              <  1:154229/151‑1 (MQ=255)                                                                                                                                     aGTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATcgc            <  2:11670/151‑1 (MQ=255)                                                                                                                                      gTTTTTCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGGCGATAAAGCGCTAAGAAATAAGCAGCAGTACGGATCGATAAGCCTGCGGGAGCGTCAGGCCCATCGCCAGCGCGGCAATTTTTTTCCCTCGCgggagggaacgaa                >  1:527805/1‑141 (MQ=255)                                                                                                                                       tttttCACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATcgcg           <  2:457850/150‑1 (MQ=255)                                                                                                                                            cACCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATCGCGGcgcg      <  1:682569/150‑1 (MQ=255)                                                                                                                                             aCCAGGCTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATCGCGGCGCGTa    <  1:641455/151‑1 (MQ=255)                                                                                                                                              ccAGGCTCATCATCTCTGAGGTCATCGGCGGGACGCTAAAGCGATAAGCAATAGGCAGCAGGACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCGTTTTTTTTCCCTCGCGGCGGCGACTGAAACGCGGCGCGGAc   >  1:467154/1‑151 (MQ=255)                                                                                                                                                aGGCTCATCATCTCTGAGGTCATCGGCGGGGAGATAACCCGATAAACCATAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCCGCGCGGCATTTTTTTTCCCTCGCGGGAGCGACTGAATCGCGGCGCGGAcc  >  1:684299/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAATCGAGCAATTTACCCGCCTGCGCCGCCATATCCACCAGACCGATAGTCGAGGCGATGGCGGAGTTTTTCACCAGGTTCATCATCTCTGAGGTCATCGGCGGGACGATAACGCGATAAGCATTAGGCAGCAGTACGTATCGATAAGCCTGCGGTAGCGTCAGGCCCATCGCCAGCGCGGCATTTTTTTGCCCTCGCGGCAGCGACTGAATCGCGGCGCGTACC  >  NC_000913/685991‑686283

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