Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,091,682	C→A	G229* (GGA→TGA)	frvB ←	putative PTS enzyme IIBC component FrvB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,091,682	0	C	A	100.0%	86.4 / NA	27	G229* (GGA→TGA)	frvB	putative PTS enzyme IIBC component FrvB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/12);  total (15/12)

CCGACCAGAAACGCTGCCGCAAACTTCATTGATGGAGGAATGGTATTCAGGAAGTGGAGCAGCCCGCCCATGAGGTCGGACATCACCGGCCCGATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATAC  >  NC_000913/4091432‑4091927                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      ccGACCAGAAACGCTGCCGCAAACGTCATTGATGGAGGAATGGTATTCAGGAAGTGGAGCAGCCCGCCCATGAGGTCGGACATCACCGGCCCGATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATcc                                                                                                                                                                                                                                                       <  2:33129/251‑1 (MQ=255)           aCGCTGCCGCAAACTTCATTGATGGAGGAATGGTATTCAGGAAGTGGAGCAGCCCGCCCATGAGGTCGGACATCACCGGCCCGATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGAc                                                                                                                                                                                                                                             >  1:168620/1‑251 (MQ=255)                                     ggaATGGTATTCAGGAAGTGGAGCAGCCCGCCCATGAGGTCGGACATCACCGGCCCGATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGtt                                                                                                                                                                                                                   >  2:57867/1‑251 (MQ=255)                                      gaATGGTATTCAGGAAGTGGAGCAGCCCGCCCATGAGGTCGGACATCACCGGGCCGATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTc                                                                                                                                                                                                                  <  2:161397/251‑1 (MQ=255)                                       aaTGGTATTCAGGAAGTGGAGCAGCCCGCCCATGAGGTCGGACATCACCGGCCCGATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCg                                                                                                                                                                                                                 >  2:27339/1‑251 (MQ=255)                                                                           gTCGGACATCACCGGCCCGATAACGGAGTAAGTGAGCACGCCGAAAACCAGCCGCGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGACGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGACGGAGCGAAAGCCg                                                                                                                                                                             >  1:144368/1‑251 (MQ=255)                                                                                          ccagATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTGTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATcg                                                                                                                                                              <  2:91244/248‑1 (MQ=255)                                                                                                aaCGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAgg                                                                                                                                                        >  1:203806/1‑251 (MQ=255)                                                                                                 aCGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGc                                                                                                                                                       >  2:93766/1‑251 (MQ=255)                                                                                                         aGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGTACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAAtata                                                                                                                                                 >  1:91866/1‑249 (MQ=255)                                                                                                                        aaCCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGa                                                                                                                                <  2:186203/251‑1 (MQ=255)                                                                                                                                gagTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTGGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATTGGGATCATAAAg                                                                                                                        <  2:146862/251‑1 (MQ=255)                                                                                                                                    aacgaacgGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGa                                                                                                                    <  1:70408/251‑1 (MQ=255)                                                                                                                                       gatcgaggtcagaatggAACCGCGCAGCGGTTGTAATTATTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCCCCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGGCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGTAGc                                                                                                                 <  2:1277/235‑1 (MQ=255)                                                                                                                                         acgGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCa                                                                                                               <  1:41498/251‑1 (MQ=255)                                                                                                                                              aTCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAGGACGAAATAGCCAACCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGTGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTa                                                                                                          <  2:167744/251‑1 (MQ=255)                                                                                                                                                              gcagcGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACt                                                                                          >  1:20553/1‑251 (MQ=255)                                                                                                                                                               cagcGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACt                                                                                          >  1:129674/1‑250 (MQ=255)                                                                                                                                                                  cGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGACACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGCTCATAAAGGTGAAGCCAAGGTAACCGATGGATTCCACTAccc                                                                                      >  2:136084/1‑251 (MQ=255)                                                                                                                                                                   ggTTGTAATGCTTTGCCAAGGCGCCCTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGACGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCa                                                                                     <  2:54065/251‑1 (MQ=255)                                                                                                                                                                                   cAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGgcgcg                                                                     >  1:38230/1‑251 (MQ=255)                                                                                                                                                                                                 cGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATc                                                       >  1:182181/1‑251 (MQ=255)                                                                                                                                                                                                          aagaCGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACgg                                              >  1:195293/1‑251 (MQ=255)                                                                                                                                                                                                           agaCGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGc                                             >  1:5672/1‑251 (MQ=255)                                                                                                                                                                                                              cGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGTTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCaaa                                          <  1:100845/251‑1 (MQ=255)                                                                                                                                                                                                              cGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCaa                                           >  1:234517/1‑250 (MQ=255)                                                                                                                                                                                                                           cGCCAGCCCCAGCACCACCGCGCCGAGGAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCAcc                              <  2:183847/250‑1 (MQ=255)                                                                                                                                                                                                                                                      gAATCACGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACGGCACCAGCATGTTGGCAAGCGCCACCAGGATac  <  1:219222/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      CCGACCAGAAACGCTGCCGCAAACTTCATTGATGGAGGAATGGTATTCAGGAAGTGGAGCAGCCCGCCCATGAGGTCGGACATCACCGGCCCGATAACGTAGTAAGTGAGCACGCCGAAAACCAGCAGAGTAACGAACGGGATCAGCATTGAACCGAGCAGCGGTTGTAATGCTTTGCCAAGGCGCACTTTGCGGAACCAGAAGACGAAATAGCCAATCGCCAGCCCCAGCACCACCGCGCCGAGGAATCCCGCACCCGACTGGGTGCCGAGCAGCGCTTTGTCGTTCGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATATCGGTATACGGCAAACCGAACTGCACCAGCATGTTGGCAAGCGCCACCAGGATAC  >  NC_000913/4091432‑4091927

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