Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,863,090	T→C	T172A (ACG→GCG)	glvC ←	putative PTS enzyme II component GlvC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,863,090	0	T	C	100.0%	230.5 / NA	61	T172A (ACG→GCG)	glvC	putative PTS enzyme II component GlvC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/31);  total (30/31)

GACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAG  >  NC_000913/3862956‑3863225                                                                                                                                       |                                                                                                                                        gACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCGATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCt                                                                                                                                    >  2:1182399/1‑140 (MQ=255) gACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCt                                                                                                                                    >  1:1223002/1‑140 (MQ=255) gACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCt                                                                                                                                    >  2:131822/1‑140 (MQ=255)        aGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAAt                                                                                                                             >  2:391090/1‑140 (MQ=255)                  ccgAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGaaaa                                                                                                                   >  2:299887/1‑140 (MQ=255)                    gAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAAc                                                                                                                 >  2:469766/1‑140 (MQ=255)                       cGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCgg                                                                                                              <  1:377642/140‑1 (MQ=255)                         cAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCa                                                                                                            <  2:419260/140‑1 (MQ=255)                               cGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGttttt                                                                                                      >  1:342288/1‑140 (MQ=255)                                  ttGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTAt                                                                                                   <  2:1188688/140‑1 (MQ=255)                                    gCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCg                                                                                                 <  1:500044/140‑1 (MQ=255)                                      agagATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGaa                                                                                               <  2:1223002/140‑1 (MQ=255)                                       gagaTTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAAc                                                                                              <  1:773511/140‑1 (MQ=255)                                       gagaTTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAAc                                                                                              <  1:512668/140‑1 (MQ=255)                                        agaTTCAATCCCCGTTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACa                                                                                             <  1:873988/140‑1 (MQ=255)                                        agaTTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACa                                                                                             <  1:335339/140‑1 (MQ=255)                                        agaTTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAACCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACa                                                                                             <  1:806943/140‑1 (MQ=255)                                          aTTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACGTAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGa                                                                                           >  2:342242/1‑140 (MQ=255)                                                    cATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCa                                                                                 >  1:69957/1‑140 (MQ=255)                                                     aTTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAg                                                                                <  2:378586/140‑1 (MQ=255)                                                          tACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCg                                                                           >  2:184540/1‑140 (MQ=255)                                                           aCTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGt                                                                          >  2:98465/1‑140 (MQ=255)                                                            cTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTc                                                                         >  1:60937/1‑140 (MQ=255)                                                             ttttGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTca                                                                        <  1:486235/140‑1 (MQ=255)                                                                   ccAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGc                                                                  <  1:718130/140‑1 (MQ=255)                                                                    cAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGcc                                                                 >  1:737956/1‑140 (MQ=255)                                                                       ccGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGa                                                              >  2:554676/1‑140 (MQ=255)                                                                       ccGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGa                                                              <  1:563811/140‑1 (MQ=255)                                                                                 tCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAAtt                                                    >  1:1065074/1‑140 (MQ=255)                                                                                 tCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAAtt                                                    >  1:922634/1‑140 (MQ=255)                                                                                       aGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCg                                              <  2:471939/140‑1 (MQ=255)                                                                                         gCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgcgc                                                <  1:696431/136‑1 (MQ=255)                                                                                         gCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgcgc                                                >  2:696431/1‑136 (MQ=255)                                                                                         gCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGAt                                            <  1:131822/140‑1 (MQ=255)                                                                                         gCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGAt                                            <  1:554676/140‑1 (MQ=255)                                                                                         gCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGAt                                            <  1:523989/140‑1 (MQ=255)                                                                                          cacaGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATa                                           >  2:1255665/1‑140 (MQ=255)                                                                                               gggATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgc                                                  >  2:121444/1‑128 (MQ=255)                                                                                               gggATCATCACCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgc                                                  <  1:121444/128‑1 (MQ=255)                                                                                                   tcatcaCCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTa                                  >  1:489391/1‑140 (MQ=255)                                                                                                       caCCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGa                              <  2:1065074/140‑1 (MQ=255)                                                                                                       caCCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGa                              <  1:98465/140‑1 (MQ=255)                                                                                                        aCCAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAg                             <  2:60937/140‑1 (MQ=255)                                                                                                          cAGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAggg                           >  2:890512/1‑140 (MQ=255)                                                                                                           aGGAAGGCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGt                          >  2:890341/1‑140 (MQ=255)                                                                                                                 gCGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAAt                    >  1:1182281/1‑140 (MQ=255)                                                                                                                  cGATAATCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATc                   >  1:224037/1‑140 (MQ=255)                                                                                                                      aaTCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCgg               <  2:342288/140‑1 (MQ=255)                                                                                                                      aaTCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCgg               <  2:409395/140‑1 (MQ=255)                                                                                                                      aaTCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCgg               <  1:890512/140‑1 (MQ=255)                                                                                                                      aaTCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCgg               <  1:299887/140‑1 (MQ=255)                                                                                                                       aTCACCACATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGc              <  1:1255665/140‑1 (MQ=255)                                                                                                                          accacATAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCcat           <  2:922634/140‑1 (MQ=255)                                                                                                                               aTAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCgg               <  2:13417/131‑1 (MQ=255)                                                                                                                               aTAAGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCgg               >  1:13417/1‑131 (MQ=255)                                                                                                                                 aaGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTc    <  2:299023/140‑1 (MQ=255)                                                                                                                                 aaGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTc    >  1:94245/1‑140 (MQ=255)                                                                                                                                 aaGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTc    >  1:391787/1‑140 (MQ=255)                                                                                                                                  aGACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCa   >  2:502988/1‑140 (MQ=255)                                                                                                                                   gACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAg  >  2:939514/1‑140 (MQ=255)                                                                                                                                   gACGCCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGCTTTAATCCCGGCCATCATTGTCAg  >  2:384827/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                        GACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAG  >  NC_000913/3862956‑3863225

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