Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	704,008	C→A	P22Q (CCG→CAG)	nagE →	N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	704,008	0	C	A	100.0%	100.5 / NA	36	P22Q (CCG→CAG)	nagE	N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/24);  total (12/24)

AGTACAATTTGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGC  >  NC_000913/703874‑704147                                                                                                                                       |                                                                                                                                             aGTACAATTTGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTggcgg                                                                                                                                      >  2:102700/1‑143 (MQ=255)  gTACAATTTGCAGCAAAATAAAAATACGGCTTGCAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTggcggc                                                                                                                                     <  2:207271/143‑1 (MQ=255)    aCAATTTGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCAc                                                                                                                                   <  2:185195/143‑1 (MQ=255)     cAATTTGCAGCAAAGTAGAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGATTATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACt                                                                                                                                  <  1:280062/143‑1 (MQ=255)       aTTTGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGt                                                                                                                                <  2:86882/143‑1 (MQ=255)        tttGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGAAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACGGtt                                                                                                                               >  1:368932/1‑143 (MQ=255)                aaaataaaaataCGGCTTGAAACGAGCCAAATAGGGTTCTCGGAGGGGGAATAAGATGAATAATTTAGGATTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCCCGGTGCTGCAGGTGGCGGACCCGGTGCTGCGAt                                                                                                                       >  1:598397/1‑143 (MQ=255)                           agggCCTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCa                                                                                                           <  1:38590/141‑1 (MQ=255)                                  gAAACGAGCCAAATCGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTGCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTAc                                                                                                    <  1:45888/143‑1 (MQ=255)                                  gAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTAc                                                                                                    <  2:606893/143‑1 (MQ=255)                                              aTAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGt                                                                                        <  2:158451/143‑1 (MQ=255)                                              aTAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTGCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGt                                                                                        <  1:641123/143‑1 (MQ=255)                                              aTAGGGTGCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGt                                                                                        <  1:644111/143‑1 (MQ=255)                                                    ttCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGCCGGGACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAAGGTTGCGTTTAtta                                                                                  >  1:437982/1‑142 (MQ=255)                                                         gTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGTTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAg                                                                             >  2:619521/1‑143 (MQ=255)                                                                             atTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGa                                                         <  2:356611/143‑1 (MQ=255)                                                                               tttAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTACCGTTACGTTTATTGCCCAGGCGGGCGGTGCGTTTTTTGAtg                                                       >  1:248684/1‑142 (MQ=255)                                                                               tttAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATa                                                       <  1:612119/143‑1 (MQ=255)                                                                                  aGGTTCTGTCCAGCGCCGCTGTAGGGCGGTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAAcc                                                    <  1:595232/143‑1 (MQ=255)                                                                                     ttttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCTGTGCTGCAGGTGGCGGCCCTGTTG‑TTCGATAGCACTCGGCTTGAGGGGCCCAGTGGTGCCTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCg                                                 <  2:170444/143‑1 (MQ=255)                                                                                     ttttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAG‑TCGGTCCGCCAGATTTACTTAACGTAGCGTTGCTTGCCCCGGCGGGCGGCGCGGTTCTGGATCGCCTCg                                                 >  1:719301/1‑143 (MQ=255)                                                                                       ttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCa                                               <  1:585654/143‑1 (MQ=255)                                                                                        tttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCAt                                              <  2:502808/143‑1 (MQ=255)                                                                                            cAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAAt                                          <  2:527754/143‑1 (MQ=255)                                                                                            cAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTCCGTTTATTGCCACGGCCGGGGGTGCGATTGTTGATAACCTCGCATTaag                                          >  1:492695/1‑142 (MQ=255)                                                                                              gCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCACGCGGCCGGTGCGATTTTTGATAACCTCGCATTAATCt                                        >  2:23417/1‑143 (MQ=255)                                                                                               cggcTCGGCAGGGCGTTACGGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATcgt                                       <  1:298365/140‑3 (MQ=255)                                                                                                gattgggTAGGGCGTTACAGCTCCCTATCGTGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTc                                      <  1:318312/138‑1 (MQ=255)                                                                                                gagtCGGGAGGGCGTTACAGCTCCCTCTCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTc                                      <  1:471471/140‑1 (MQ=255)                                                                                                  tggggTAGGGCGTTACCGCCCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGc                                    <  1:661196/140‑1 (MQ=255)                                                                                                         gggCGTTACTGCTCCCTATCGCGGGGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGgt                             <  1:694461/143‑1 (MQ=255)                                                                                                                aCAGCTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTGACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCAt                      >  1:173309/1‑143 (MQ=255)                                                                                                                    cTCCCTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGAGTTACTTAACGTTGCGTTGATTGCCCAGGCGGGCGGTGCGGCATTTGATAACATCGCGCTACTCTGCGCCATCTGGGTGGCATcctg                  >  1:410997/1‑141 (MQ=255)                                                                                                                     gcccTATCGCGGTGCTGCAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGc                 <  1:193472/142‑1 (MQ=255)                                                                                                                             gcgGTGCTGCATGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGaa         <  1:77828/143‑1 (MQ=255)                                                                                                                                    tgcAGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAgc  >  2:414152/1‑143 (MQ=255)                                                                                                                                       |                                                                                                                                             AGTACAATTTGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGAT‑TCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGC  >  NC_000913/703874‑704147

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