Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,507,334	T→C	L187L (TTA→TTG)	fryC ←	putative PTS enzyme IIC component FryC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,507,334	0	T	C	100.0%	158.1 / NA	45	L187L (TTA→TTG)	fryC	putative PTS enzyme IIC component FryC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/24);  total (21/24)

TACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGTAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGT  >  NC_000913/2507210‑2507458                                                                                                                             |                                                                                                                             tACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTgg                                                                                                               >  1:530895/1‑140 (MQ=255)         ccgccgTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCa                                                                                                       <  2:177317/140‑1 (MQ=255)         ccgccgTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTCTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCa                                                                                                       >  2:131852/1‑140 (MQ=255)           gccgTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATcc                                                                                                        >  2:478131/1‑137 (MQ=255)           gccgTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATcc                                                                                                        <  1:478131/137‑1 (MQ=255)           gccgTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACt                                                                                                     >  1:500066/1‑140 (MQ=255)                tGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTcaca                                                                                                <  2:530895/140‑1 (MQ=255)                       gCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAAc                                                                                         >  2:556560/1‑140 (MQ=255)                            caAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTcac                                                                                                 >  1:491088/1‑127 (MQ=255)                            caAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTcac                                                                                                 <  2:491088/127‑1 (MQ=255)                            caAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCa                                                                                    >  2:428336/1‑140 (MQ=255)                              aaGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCCGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCAGCAAc                                                                                  <  2:736909/140‑1 (MQ=255)                                gggggTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGa                                                                                <  1:428336/140‑1 (MQ=255)                                 ggggTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTaa                                                                                          <  1:396236/129‑1 (MQ=255)                                 ggggTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTaa                                                                                          >  2:396236/1‑129 (MQ=255)                                  gggTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGa                                                                                >  2:628131/1‑138 (MQ=255)                                  gggTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGa                                                                                <  1:628131/138‑1 (MQ=255)                                      gatgaCGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGaaa                                                                          >  1:368004/1‑140 (MQ=255)                                             tagtaCATCGCCAGCATGACAAGAATGGCGGGAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAg                                                                   <  2:321301/140‑1 (MQ=255)                                                taCATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTTGTTCATCCGCTTTACAGGGTAACCGGCAACGATAGAAATAATCAgcgc                                                                >  2:263745/1‑140 (MQ=255)                                                  cATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCAc                                                              >  1:364394/1‑140 (MQ=255)                                                  cATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCAc                                                              <  2:500066/140‑1 (MQ=255)                                                     cGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCaa                                                           >  2:619549/1‑140 (MQ=255)                                                     cGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCaa                                                           >  2:518183/1‑140 (MQ=255)                                                       ccAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATTCACTTCACAAGGTAACGGGCAACGAAAGAAATAATCAGCGCACCAATg                                                         >  2:291550/1‑140 (MQ=255)                                                          gCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAg                                                      <  1:575208/140‑1 (MQ=255)                                                          gCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAg                                                      >  1:438457/1‑140 (MQ=255)                                                                     aTGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCAc                                           >  1:720/1‑140 (MQ=255)                                                                         cGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGc                                       >  1:36046/1‑140 (MQ=255)                                                                                aaTCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAgg                                <  1:556560/140‑1 (MQ=255)                                                                                       aGCAAAAATGTGGTTTTGAACGTTAACAGTAAATCCGCCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCa                         <  1:263745/140‑1 (MQ=255)                                                                                       aGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCa                         >  2:394811/1‑140 (MQ=255)                                                                                       aGCAAAAATGTGGTTATGAACGGTAGCAGCAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCa                         <  1:310781/140‑1 (MQ=255)                                                                                          aaaaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTc                              >  1:416612/1‑132 (MQ=255)                                                                                          aaaaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTc                              <  2:416612/132‑1 (MQ=255)                                                                                           aaaaTGTGGTTTTGAACGCTAATAGCAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGc                     <  1:291550/140‑1 (MQ=255)                                                                                            aaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCa                    <  2:438457/140‑1 (MQ=255)                                                                                            aaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCa                    <  1:670925/140‑1 (MQ=255)                                                                                            aaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCa                    <  2:541257/140‑1 (MQ=255)                                                                                            aaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCa                    <  1:619549/140‑1 (MQ=255)                                                                                            aaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCa                    <  1:518183/140‑1 (MQ=255)                                                                                            aaaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCa                    <  1:394811/140‑1 (MQ=255)                                                                                             aaTGTGGTTTTGAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCAt                   <  1:675934/140‑1 (MQ=255)                                                                                                        gAACGCTAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGAt        <  2:364394/140‑1 (MQ=255)                                                                                                              tAACAGGAAATCCGGCAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGt  >  1:3328/1‑140 (MQ=255)                                                                                                                             |                                                                                                                             TACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAAAATGTGGTTTTGAACGCTAACAGGAAATCCGGTAACTGGATTTTCTGGTTCATCCACTTCACAAGGTAACCGGCAACGATAGAAATAATCAGCGCACCAATGAAGGTGGACGGCACCGGCGAAGAGGTCGCCCATTGCATAGTGCTGGGATCGAAGT  >  NC_000913/2507210‑2507458

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