Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,838,055	Δ1 bp	coding (509/1458 nt)	ascF →	PTS system enzyme IIBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,055	0	G	.	100.0%	389.0 / NA	86	coding (509/1458 nt)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (38/48);  total (38/48)

TTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATT  >  NC_000913/2837985‑2838125                                                                        |                                                                       ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCa                                                                   <  1:1296700/76‑1 (MQ=255)  tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaa                                                                  >  1:999962/1‑76 (MQ=255)  tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaa                                                                  <  1:2981009/76‑1 (MQ=255)  tAACCATTCTGAATGGGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaa                                                                  <  2:874067/76‑1 (MQ=255)   aaCCATTCTGAATGTGATTGGTGACTGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                 <  2:3207135/76‑1 (MQ=255)    aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAt                                                                >  1:3793513/1‑76 (MQ=255)     ccATTCTGAATGTGATTGGGGACGGGGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAtt                                                               >  1:3484427/1‑76 (MQ=255)      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAttt                                                              <  1:4455715/76‑1 (MQ=255)      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAttt                                                              <  1:3706657/76‑1 (MQ=255)       aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTa                                                             <  2:2454093/76‑1 (MQ=255)        ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaa                                                            <  1:3768307/76‑1 (MQ=255)        ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaa                                                            >  1:613431/1‑76 (MQ=255)        ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaa                                                            <  1:3036212/76‑1 (MQ=255)        ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaa                                                            <  2:3878653/76‑1 (MQ=255)         tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaa                                                           <  1:217899/76‑1 (MQ=255)          cTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaaa                                                          <  1:2638000/76‑1 (MQ=255)          cTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaaa                                                          >  1:777589/1‑76 (MQ=255)           tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAc                                                         <  2:3637707/76‑1 (MQ=255)            gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAcc                                                        >  2:3543338/1‑76 (MQ=255)            gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAcc                                                        >  1:1699902/1‑76 (MQ=255)             aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                       <  2:2523000/76‑1 (MQ=255)              aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCaa                                                      >  2:461312/1‑76 (MQ=255)               tgtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                     <  1:3077426/76‑1 (MQ=255)                gtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                    <  2:3072764/76‑1 (MQ=255)                  ggttGGTGCCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                  <  2:1763100/74‑1 (MQ=255)                  gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                  >  2:3547640/1‑76 (MQ=255)                   gttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGt                                                 <  2:3251922/75‑1 (MQ=255)                     tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                               <  2:4420961/76‑1 (MQ=255)                     tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                               <  2:579995/76‑1 (MQ=255)                      ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                              <  1:978791/76‑1 (MQ=255)                      ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                              >  2:3111124/1‑76 (MQ=255)                       gTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                             <  2:999961/76‑1 (MQ=255)                        tGACGGTGTTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                            <  2:2304784/76‑1 (MQ=255)                        tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                            <  1:2127254/76‑1 (MQ=255)                        tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                            >  1:2740332/1‑76 (MQ=255)                        tGACGGTGCTTTCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                            <  2:4332591/76‑1 (MQ=255)                        tGACGATGCTTTCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAGCCAATAGGTCGCTg                                            <  2:3685766/76‑1 (MQ=255)                         gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                           <  1:2118851/76‑1 (MQ=255)                          aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgtc                                          >  1:3076626/1‑74 (MQ=255)                           cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCg                                         <  1:626662/76‑1 (MQ=255)                               gCTTTCTTCTTCCGGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                     <  2:4096696/76‑1 (MQ=255)                                cTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGc                                    >  1:4324058/1‑76 (MQ=255)                                cTTTCTTCTGCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGc                                    >  2:2841257/1‑76 (MQ=255)                                  ttcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                  >  1:1204281/1‑76 (MQ=255)                                   tcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                 <  1:266810/76‑1 (MQ=255)                                   tcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                 <  1:3513670/76‑1 (MQ=255)                                     ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                               >  2:1239943/1‑76 (MQ=255)                                       cttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGCTTGCCATctcg                             >  2:910995/1‑76 (MQ=255)                                       cttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                             >  1:3241408/1‑76 (MQ=255)                                       cttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                             <  1:3823850/76‑1 (MQ=255)                                       ctcccTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                             <  1:139371/73‑1 (MQ=255)                                           cTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                         <  1:2841257/76‑1 (MQ=255)                                           cTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                         >  1:2746657/1‑76 (MQ=255)                                            tGCCGCTGATGGTGGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAGCGCGGgtgt                        <  1:1186401/76‑1 (MQ=255)                                            tGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgg                        >  1:2299922/1‑75 (MQ=255)                                            cgCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAACCCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                        <  1:1985620/75‑1 (MQ=255)                                              ccGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCCTCGCGGGGGTGc                      >  2:204089/1‑76 (MQ=255)                                               cGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTCAAACCAATATGTCGCTGGCGATTGCCCTCGCGGGGGTGCt                     >  2:2855938/1‑76 (MQ=255)                                                gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                    >  2:4402161/1‑76 (MQ=255)                                                gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                    <  2:4416321/76‑1 (MQ=255)                                                gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                    >  1:3886689/1‑76 (MQ=255)                                                  tGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                  <  1:1288219/76‑1 (MQ=255)                                                  tGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                  <  1:3850817/76‑1 (MQ=255)                                                  tGATGGGCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                  <  2:278859/76‑1 (MQ=255)                                                   gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                 >  2:2541400/1‑76 (MQ=255)                                                   gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                 >  2:2076758/1‑76 (MQ=255)                                                   gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGGGtggtgggg                 >  2:751686/1‑74 (MQ=255)                                                    agggTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                <  2:225647/74‑1 (MQ=255)                                                     tGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa               <  1:2731360/76‑1 (MQ=255)                                                     tGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa               >  1:4073616/1‑76 (MQ=255)                                                       gTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc             >  2:914573/1‑76 (MQ=255)                                                       gTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGGGCTGGTACCTc             >  2:76469/1‑76 (MQ=255)                                                           cggcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg         >  2:1639114/3‑76 (MQ=255)                                                            ggcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt        <  2:3679926/75‑1 (MQ=255)                                                            ggcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt        >  1:4300060/2‑76 (MQ=255)                                                             gcATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt       <  1:3457680/76‑1 (MQ=255)                                                              cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttt      <  1:1133104/76‑1 (MQ=255)                                                              cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttt      <  1:4343745/76‑1 (MQ=255)                                                               aTCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtttt     >  2:2979605/1‑76 (MQ=255)                                                               aTCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtttt     >  1:3999815/1‑76 (MQ=255)                                                               aTCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGGACATCCGAGtttt     >  2:3560744/1‑76 (MQ=255)                                                                tCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa    <  1:1926612/76‑1 (MQ=255)                                                                 cTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt   <  2:3552892/76‑1 (MQ=255)                                                                  tGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt  >  2:4209183/1‑76 (MQ=255)                                                                  tGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt  >  1:4205091/1‑76 (MQ=255)                                                                  tGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGAGGGTGTGCTGGTACATCCGAGTTTTAtt  >  1:1870185/1‑76 (MQ=255)                                                                        |                                                                       TTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATT  >  NC_000913/2837985‑2838125

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