Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,227,995	C→T	D105D (GAC→GAT)	ycgL →	PF05166 family protein YcgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,227,995	0	C	T	100.0%	98.8 / NA	29	D105D (GAC→GAT)	ycgL	PF05166 family protein YcgL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/15);  total (14/15)

GAAGAAACTGGTCAATGCCGATATTGAAAAGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGACACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAG  >  NC_000913/1227860‑1228128                                                                                                                                        |                                                                                                                                      gaagaaACTGGTCAATGCCGATATTGAAAAGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTa                                                                                                                                   <  1:811854/140‑1 (MQ=255)   agaaACTTGTCAATGCCGATATTGAAAAGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAAc                                                                                                                                 >  2:524255/1‑140 (MQ=255)                tGCCGATATTGAAAAGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGatat                                                                                                                    >  1:691378/1‑140 (MQ=255)                tGCCGATATTGAAAAGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGatat                                                                                                                    >  1:725730/1‑140 (MQ=255)                tGCCGATATTGAAAAGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGatat                                                                                                                    >  1:844073/1‑140 (MQ=255)                            aaaGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCa                                                                                                        <  2:648476/140‑1 (MQ=255)                                          cATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGgc                                                                                          >  1:745233/1‑140 (MQ=255)                                          cATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGgc                                                                                          >  2:454859/1‑140 (MQ=255)                                                   aGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGtttt                                                                                 >  1:504540/1‑140 (MQ=255)                                                   aGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGtttt                                                                                 >  1:105368/1‑140 (MQ=255)                                                     cAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCa                                                                               <  2:691378/140‑1 (MQ=255)                                                        ggTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCatat                                                                            <  1:196821/140‑1 (MQ=255)                                                          ttACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATcc                                                                          <  2:380300/140‑1 (MQ=255)                                                           tACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCa                                                                         >  1:504809/1‑140 (MQ=255)                                                              tATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTc                                                                                >  2:86637/1‑130 (MQ=255)                                                              tATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTc                                                                                <  1:86637/130‑1 (MQ=255)                                                                  tGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAg                                                                  >  2:336091/1‑140 (MQ=255)                                                                            gCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGaa                                                        <  2:504809/140‑1 (MQ=255)                                                                              caccacCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACg                                                      <  2:906759/140‑1 (MQ=255)                                                                               accaccCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGt                                                     <  2:844073/140‑1 (MQ=255)                                                                                            ttGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACAt                                        >  2:685042/1‑140 (MQ=255)                                                                                                      aaCATCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGc                              <  2:725730/140‑1 (MQ=255)                                                                                                        aaTCTTTCCGTCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCaa                            <  2:363498/139‑1 (MQ=255)                                                                                                                  tCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTg                        <  1:442735/134‑1 (MQ=255)                                                                                                                  tCATGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTg                        >  2:442735/1‑134 (MQ=255)                                                                                                                    aTGGGGCAGAAAACAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGctgctg                >  2:824462/1‑140 (MQ=255)                                                                                                                             aaaaCAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCg       <  2:542751/140‑1 (MQ=255)                                                                                                                             aaaaCAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCg       <  2:46174/140‑1 (MQ=255)                                                                                                                               aaCAGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGt     <  1:824462/140‑1 (MQ=255)                                                                                                                                  aGACGATACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAg  <  2:105368/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      GAAGAAACTGGTCAATGCCGATATTGAAAAGGTTAAACAGGCATTAACCGAGCAAGGTTACTATTTGCAGTTACCGCCACCACCCGAAGATTTGCTGAAGCAACATCTTTCCGTCATGGGGCAGAAAACAGACGACACTAACAAATAACCGATATCCGGCGGTGGCATTATCTTTGTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAG  >  NC_000913/1227860‑1228128

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