Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,529,039	C→T	G402D (GGT→GAT)	sgcC ←	putative PTS enzyme IIC component SgcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,529,039	0	C	T	100.0%	114.0 / NA	35	G402D (GGT→GAT)	sgcC	putative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/13);  total (22/13)

GCTTTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGG  >  NC_000913/4528793‑4529283                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     gCTTTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGa                                                                                                                                                                                                                                                   >  2:54582/1‑250 (MQ=255) gCTTTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGa                                                                                                                                                                                                                                                   >  1:42288/1‑250 (MQ=255)     tGTCGATAAACCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGAACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACat                                                                                                                                                                                                                                               >  2:306820/1‑250 (MQ=255)                     cATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCCCGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATCCCGATAATTCACGTTTACTTAAGAACAGCACGACGCTGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACg                                                                                                                                                                                                                              <  1:277797/251‑1 (MQ=255)                        ccccacttGCGCATCAAAACTGGGTTCGCGGGGCAATGCAGGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCTTTATTCCGATAATTCACGTTTACGTAATAACAACACGACGCTGGCGACAAGTCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATGCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTg                                                                                                                                                                                                                           <  1:142472/244‑1 (MQ=255)                                   cataaaaccacGGGTCGCCGGGCAATGCACATAGATGGCACAGGGCAATGACCGTTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACCCCTTATTCCGAGAATTCACGTTTACGTAAGACCAGCACGACGCTGGCAACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCa                                                                                                                                                                                                                <  2:161016/241‑1 (MQ=255)                                       aaaaCTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACttt                                                                                                                                                                                                            >  1:226943/1‑251 (MQ=255)                                                        gCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCAAATACCCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGaaa                                                                                                                                                                                           <  1:248278/251‑1 (MQ=255)                                                              cacttctctttcccatagCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCGAATCATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATATCCCCCTGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCt                                                                                                                                                                                     >  1:119565/18‑251 (MQ=255)                                                                                           ttttCCGTCACAATAACCTCTTTCAGCCAAC‑TCATTGACAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTACGgcgc                                                                                                                                                        >  2:43273/1‑251 (MQ=255)                                                                                           ttttCAGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGCAGTAAGtcgc                                                                                                                                                        >  2:235173/1‑251 (MQ=255)                                                                                                   cccAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTggg                                                                                                                                                >  1:299211/1‑251 (MQ=255)                                                                                                    ccAATACCCTCTTTCA‑CCAACCTCAGTAAAAGACCTCCTTATCCCGAAAATTCGCTTTTACTAAAGAACAGCACAACCCTGGCGAAATTCCCGCCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTTAAGTAAGGCGCGAACTGGGt                                                                                                                                               <  2:305075/251‑1 (MQ=255)                                                                                                       aTAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGc                                                                                                                                            >  2:299996/1‑251 (MQ=255)                                                                                                              ctTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATcagc                                                                                                                                     >  1:139890/1‑251 (MQ=255)                                                                                                                       aaC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCt                                                                                                                                                                                     >  2:109587/1‑194 (MQ=255)                                                                                                                            cATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAatc                                                                                                                         >  2:5031/1‑250 (MQ=255)                                                                                                                                 aaaGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACg                                                                                                                   >  1:256838/1‑251 (MQ=255)                                                                                                                                   aGTACTCCTTATTCCGATAATTCACGTTTACGTAACAACAGCACGACGCTGGCGCAAATACAGACAGCCACGACAACGCCACTCATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACACCGCTGAGCTTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGCACTGGGTGGCGATCAGCAGGACGGTAATCATCACGAt                                                                                                                 >  2:276194/1‑251 (MQ=255)                                                                                                                                     tACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATg                                                                                                                >  1:238750/1‑250 (MQ=255)                                                                                                                                              ttCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCg                                                                                                                                                                                                       >  2:248525/1‑154 (MQ=255)                                                                                                                                              ttCCGATAATTCACGTTTACGTAAGAACAGCACGACGCGGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCg                                                                                                                                                                                                       <  1:248525/154‑1 (MQ=255)                                                                                                                                                            gTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGAcc                                                                                        >  2:163069/1‑251 (MQ=255)                                                                                                                                                                acgaaaGAACAGCACGACGCTGGCGAAAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAgg                                                                                     <  1:110607/246‑1 (MQ=255)                                                                                                                                                                              cgaggcggggcaaaaccccgcccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATa                                                                      <  2:139890/230‑1 (MQ=255)                                                                                                                                                                                   ggggCGACAACCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGt                                                                 <  2:164234/249‑1 (MQ=255)                                                                                                                                                                                         aCAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGAcc                                                                                        >  2:285659/1‑222 (MQ=255)                                                                                                                                                                                             ccccgaccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAg                                                       <  2:226943/249‑1 (MQ=255)                                                                                                                                                                                                 gaccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGCTAACCGTAGCCATACAGATa                                                    >  1:228017/1‑250 (MQ=255)                                                                                                                                                                                                        aCGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCTAACTTATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTACTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAActt                                            <  2:7031/251‑3 (MQ=255)                                                                                                                                                                                                             aaCGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATAATGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTAAGTGAAGTAATGCGCGAACTGGGTGGCGATCAGCATGACGGTAATCATCACGATGACCACGCTTATCAACGTCCTTACCAGGTCTCCGCGGTGGATAACCGTATCCATACAGATCAAGAACGGttct                                        >  2:100947/1‑247 (MQ=255)                                                                                                                                                                                                               cGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTg                                                                                                                                                                                                                    >  2:117102/1‑76 (MQ=255)                                                                                                                                                                                                               cGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTAcc                                     >  2:28871/1‑251 (MQ=255)                                                                                                                                                                                                                    agaatCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACAGATGTCCCACGGACAACGCGGAGATTTGCGCATTTTCACCGGAAAAATTGAAACCGCCTTTGAGCGCCATTTCAGTTAAGTAAGGCGCTAACTGGGTGGCGATCAGCATGACGGTAATCATCACGATGACCCCTCTTATCAACGTTCTGACCAGGCCTCCGCGGTGGATAACCGTGGCCATACATATAAAGAACGGTGCTACCGGCAg                                <  2:192022/249‑1 (MQ=255)                                                                                                                                                                                                                                                  gACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCgg  <  1:68475/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     GCTTTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAAC‑TCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGG  >  NC_000913/4528793‑4529283

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