Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,507,885	(T)6→5	coding (10/1248 nt)	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,880	0	T	.	100.0%	119.0 / NA	29	coding (15/1248 nt)	fryC	putative PTS enzyme IIC component FryC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (9/20);  total (9/20)

CAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGA  >  NC_000913/2507800‑2507956                                                                                 |                                                                             cAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGt                                                                              <  1:445297/81‑1 (MQ=255)      gCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGttttta                                                                         <  2:445298/81‑2 (MQ=255)       cTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtttttaa                                                                        <  1:445299/81‑3 (MQ=255)          gTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtttttaat                                                                       <  1:445301/79‑4 (MQ=255)          gTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATgg                                                                    <  1:445300/81‑1 (MQ=255)              tAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑GTTTTAATGGCCa                                                                 >  2:445347/1‑80 (MQ=255)               agagGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCAt                                                                <  1:445302/80‑1 (MQ=255)               agagGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCAt                                                                <  1:445303/80‑1 (MQ=255)               agagGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCAt                                                                <  1:445304/80‑1 (MQ=255)                 agGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATaa                                                              >  2:483076/1‑80 (MQ=255)                 agGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAAc                                                             <  2:445305/81‑1 (MQ=255)                         ggATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATcc                                                     <  2:445306/81‑1 (MQ=255)                                 aaCTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTc                                                <  1:445308/78‑1 (MQ=255)                                 aaCTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGa                                              <  1:445307/80‑1 (MQ=255)                                      ccgctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAatat                                          <  1:445310/79‑1 (MQ=255)                                      ccgctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                        <  2:445309/81‑1 (MQ=255)                                          tgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgt                                    >  2:445349/1‑81 (MQ=255)                                            ccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgt                                    <  2:445311/79‑1 (MQ=255)                                              gcACCGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTgttgtt                                <  2:445312/81‑1 (MQ=255)                                                  cGGATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCaga                            <  2:445313/81‑1 (MQ=255)                                                    gATGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAgg                          <  1:445314/81‑1 (MQ=255)                                                     aTGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGc                         <  2:445315/81‑1 (MQ=255)                                                      tGCGCCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCa                        <  2:445316/81‑1 (MQ=255)                                                        cgcCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCaa                       <  2:445317/80‑1 (MQ=255)                                                        cgcCAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAat                      >  1:445350/1‑81 (MQ=255)                                                          ccAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAatca                    <  2:445318/81‑1 (MQ=255)                                                          ccAGGCACAACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAatca                    >  2:445348/1‑81 (MQ=255)                                                               cacaACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATctct               >  2:445358/1‑81 (MQ=255)                                                               cacaACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATctct               >  2:445362/1‑81 (MQ=255)                                                                 caACGGTTGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTc             <  2:445319/81‑1 (MQ=255)                                                                       ttGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTc       >  2:445342/1‑81 (MQ=255)                                                                       ttGCACTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTc       <  1:445320/81‑1 (MQ=255)                                                                            cTGCG‑TTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtga  >  1:445344/1‑81 (MQ=255)                                                                                 |                                                                             CAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGA  >  NC_000913/2507800‑2507956

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