Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,091,800	Δ1 bp	coding (567/1452 nt)	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,797	0	G	.	100.0%	185.7 / NA	44	coding (570/1452 nt)	frvB	putative PTS enzyme IIBC component FrvB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (24/20);  total (24/20)

CGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATAT  >  NC_000913/4091719‑4091873                                                                               |                                                                             cGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATgg                                                                             >  2:828901/1‑80 (MQ=255) cGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATgg                                                                             >  1:826828/1‑80 (MQ=255)  gCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATggga                                                                           <  2:826760/81‑2 (MQ=255)    cAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATgggat                                                                          >  2:826804/1‑78 (MQ=255)     aGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATgggatca                                                                        >  1:826833/1‑77 (MQ=255)       aTAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCAt                                                                      >  2:826813/1‑80 (MQ=255)         aGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATa                                                                     >  2:826847/1‑79 (MQ=255)         aGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATa                                                                     >  2:826814/1‑79 (MQ=255)          gCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAg                                                                  <  1:826761/81‑1 (MQ=255)            aaaCCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGt                                                                <  1:826762/81‑1 (MQ=255)             aaCCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTg                                                               <  2:826763/81‑1 (MQ=255)               ccAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTg                                                               >  1:826826/1‑79 (MQ=255)                  gAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGcc                                                          <  1:826764/81‑1 (MQ=255)                   aaaCGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGcc                                                          >  1:826839/1‑80 (MQ=255)                       gCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAgg                                                      >  1:826832/1‑80 (MQ=255)                          ggAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAAc                                                  <  1:826765/81‑1 (MQ=255)                           gAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAAcc                                                 >  1:826830/1‑81 (MQ=255)                           gAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAAcc                                                 <  1:826766/81‑1 (MQ=255)                           gAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAAcc                                                 <  1:826767/81‑1 (MQ=255)                            aGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAAcc                                                 >  2:826848/1‑80 (MQ=255)                            aGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCg                                                <  1:826768/81‑1 (MQ=255)                             gCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGa                                               >  1:826851/1‑81 (MQ=255)                              cGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGAt                                              >  2:826850/1‑81 (MQ=255)                                   gCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGAtt                                         <  2:826769/81‑1 (MQ=255)                                      ggTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCaa                                      >  2:826840/1‑81 (MQ=255)                                            tCGGCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTAcc                                 >  2:826854/1‑80 (MQ=255)                                              ggCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAg                              >  1:826834/1‑81 (MQ=255)                                               gCAATCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAgg                             >  2:826852/1‑81 (MQ=255)                                                 aaTCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGt                            >  1:826837/1‑80 (MQ=255)                                                 aaTCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTa                           <  1:826770/81‑1 (MQ=255)                                                  aTCGACGAGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTa                           >  1:826856/1‑80 (MQ=255)                                                    cgacgaGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAac                        >  2:826855/1‑81 (MQ=255)                                                     gacgaGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAac                        >  1:826831/1‑80 (MQ=255)                                                     gacgaGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAac                        >  2:826841/1‑80 (MQ=255)                                                      acgaGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAacga                      >  2:826838/1‑81 (MQ=255)                                                        gaGGCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAacgacg                    >  1:826861/1‑81 (MQ=255)                                                          ggCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGgc                  <  1:826771/81‑1 (MQ=255)                                                           gCAATATAAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGgcg                 <  2:826772/81‑1 (MQ=255)                                                               tataAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGcc             <  2:826774/81‑1 (MQ=255)                                                               tataAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGcc             <  1:826773/81‑1 (MQ=255)                                                                ataAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCt            >  2:826865/1‑81 (MQ=255)                                                                 taAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCtt           >  1:826853/1‑81 (MQ=255)                                                                 taAGCGCCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCt            <  1:826775/80‑1 (MQ=255)                                                                    gcgcCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTc         >  2:826842/1‑80 (MQ=255)                                                                    gcgcCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCg        <  1:826776/81‑1 (MQ=255)                                                                     cgcCCATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGa       <  1:826777/81‑1 (MQ=255)                                                                       cccATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACa     <  2:826778/81‑1 (MQ=255)                                                                        ccATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACat    <  2:826779/81‑1 (MQ=255)                                                                         cATGAT‑GGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACata   <  1:826780/81‑1 (MQ=255)                                                                          atgat‑gGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACatat  >  1:826857/1‑81 (MQ=255)                                                                               |                                                                             CGCCAGATAGCAAACCAGAAACGCCGGAGCGAAAGCCGGTTTATCGGCAATCGACGAGGCAATATAAGCGCCCATGATGGGGATCATAAAGGTGAAGCCAAGGTAACCGATGGATTCAACTACCCAGGTAAACGACGGCGCGCCTTTCGACATAT  >  NC_000913/4091719‑4091873

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