Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,788,944	T→C	F282L (TTT→CTT)	yibQ →	putative polysaccharide deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,788,944	0	T	C	100.0%	142.7 / NA	37	F282L (TTT→CTT)	yibQ	putative polysaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/19);  total (18/19)

GCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGA  >  NC_000913/3788811‑3789088                                                                                                                                      |                                                                                                                                                 gCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGc                                                                                                                                    <  1:194668/148‑1 (MQ=255) gCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGa                                                                                                                                  <  2:652923/150‑1 (MQ=255)     aatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCa                                                                                                                             <  1:197001/151‑1 (MQ=255)      atgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCAt                                                                                                                            >  1:318974/1‑151 (MQ=255)        gaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATc                                                                                                                           <  1:309634/150‑1 (MQ=255)               aGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCag                                                                                                                     >  1:458416/1‑149 (MQ=255)                     aTACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCt                                                                                                             <  2:125502/151‑1 (MQ=255)                       aCTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGa                                                                                                           <  2:350233/151‑1 (MQ=255)                                ccgccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACt                                                                                                  >  1:104687/1‑151 (MQ=255)                                 cgccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACtt                                                                                                 >  1:575509/1‑151 (MQ=255)                                  gccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTc                                                                                                >  2:338590/1‑151 (MQ=255)                                   ccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTcc                                                                                               >  1:213953/1‑151 (MQ=255)                                         aaCGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAACAt                                                                                          <  1:282664/150‑1 (MQ=255)                                         aaCGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAACAt                                                                                          >  2:282663/1‑150 (MQ=255)                                               gtgCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTg                                                                                     <  1:520590/149‑1 (MQ=255)                                                                tAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAAc                                                                    >  2:322621/1‑149 (MQ=255)                                                                 aaTCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGc                                                                 <  2:380056/151‑1 (MQ=255)                                                                      gttccgtGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGa                                                             >  1:627331/1‑150 (MQ=255)                                                                        tccgtGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAg                                                           <  2:434999/150‑1 (MQ=255)                                                                          cgtGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAgcgg                                                        >  1:649153/1‑151 (MQ=255)                                                                             ggCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGa                                                     <  2:86218/151‑1 (MQ=255)                                                                                   aaGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTa                                               <  1:185198/151‑1 (MQ=255)                                                                                         tgCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCt                                          >  2:599990/1‑150 (MQ=255)                                                                                            aaaCCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAgcg                                      >  2:220462/1‑151 (MQ=255)                                                                                              aCCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCa                                    <  1:19941/151‑1 (MQ=255)                                                                                               ccGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCaa                                   <  1:351327/151‑1 (MQ=255)                                                                                                  aagaaACCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAAtt                                >  2:187204/1‑151 (MQ=255)                                                                                                     aaaCCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTaa                              >  2:53722/1‑150 (MQ=255)                                                                                                       aCCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTaa                              <  1:605702/148‑1 (MQ=255)                                                                                                       aCCGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTaa                              >  2:605708/1‑148 (MQ=255)                                                                                                        ccGATAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGgc                          <  2:506632/151‑1 (MQ=255)                                                                                                           aTAGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCt                       <  2:347802/151‑1 (MQ=255)                                                                                                             aGAACCTGTTTATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTAc                     >  1:653716/1‑151 (MQ=255)                                                                                                                       tATGCTAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCg           >  1:23494/1‑151 (MQ=255)                                                                                                                            tAATCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGt      >  1:614013/1‑151 (MQ=255)                                                                                                                             aaTCGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATgg       <  2:121004/149‑1 (MQ=255)                                                                                                                                cGCTTCCTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa  <  2:207318/151‑1 (MQ=255)                                                                                                                                      |                                                                                                                                                 GCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGA  >  NC_000913/3788811‑3789088

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