Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,995,912	C→T	H110H (CAC→CAT)	yigA →	conserved protein YigA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,912	0	C	T	100.0%	138.2 / NA	39	H110H (CAC→CAT)	yigA	conserved protein YigA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/16);  total (23/16)

CTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGC  >  NC_000913/3995784‑3996045                                                                                                                                 |                                                                                                                                      cTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGgcgcg                                                                                                                            <  1:584060/140‑1 (MQ=255)        agaACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGAtctc                                                                                                                     >  1:68872/1‑140 (MQ=255)               ggCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTgg                                                                                                              >  1:167570/1‑140 (MQ=255)                gcgcTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGc                                                                                                             <  1:76469/140‑1 (MQ=255)                    tGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGt                                                                                                         >  2:417765/1‑140 (MQ=255)                             aaCAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTg                                                                                                <  2:167570/140‑1 (MQ=255)                                aGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTgcgc                                                                                             >  2:155206/1‑140 (MQ=255)                                 ggCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGcc                                                                                            <  1:319748/140‑1 (MQ=255)                                  gCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTTTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCt                                                                                           >  1:472612/1‑140 (MQ=255)                                    tATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGt                                                                                         >  1:497900/1‑140 (MQ=255)                                      tCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTgcgc                                                                                             >  2:237171/1‑134 (MQ=255)                                      tCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTgcgc                                                                                             <  1:237171/134‑1 (MQ=255)                                                            tCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTg                                                                 >  2:144762/1‑140 (MQ=255)                                                               ccTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTgcgc                                                              <  1:386013/140‑1 (MQ=255)                                                               ccTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTgcgc                                                              >  2:386013/1‑140 (MQ=255)                                                               ccTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTgc                                                                >  2:326912/1‑138 (MQ=255)                                                               ccTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTgc                                                                <  1:326912/138‑1 (MQ=255)                                                                       aCCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAAc                                                      <  2:68872/140‑1 (MQ=255)                                                                             aGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGTAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACt                                                >  1:531471/1‑140 (MQ=255)                                                                             aGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACt                                                >  1:327934/1‑140 (MQ=255)                                                                              gcgcAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACcaca                                                  >  2:216807/1‑137 (MQ=255)                                                                              gcgcAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACcaca                                                  <  1:216807/137‑1 (MQ=255)                                                                              gcgcAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTc                                               >  2:470672/1‑140 (MQ=255)                                                                              gcgcAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTc                                               >  1:108128/1‑140 (MQ=255)                                                                                    tctcACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTgg                                         >  1:103779/1‑140 (MQ=255)                                                                                    tctcACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTgg                                         >  2:478513/1‑140 (MQ=255)                                                                                         ccgccgCCAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTa                                    <  2:497900/140‑1 (MQ=255)                                                                                              gccAGCAGTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCg                               >  1:294157/1‑140 (MQ=255)                                                                                                 agcagTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCa                            <  2:327934/140‑1 (MQ=255)                                                                                                 agcagTCTCGACGATATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCa                            <  1:478513/140‑1 (MQ=255)                                                                                                             gaTATGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAAcc                <  2:294157/140‑1 (MQ=255)                                                                                                                 tGCTGATGCGCTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTg            <  2:441467/140‑1 (MQ=255)                                                                                                                         cgcTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCa    >  1:346439/1‑140 (MQ=255)                                                                                                                         cgcTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCa    >  2:569580/1‑140 (MQ=255)                                                                                                                           cTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGaa                  >  2:24131/1‑124 (MQ=255)                                                                                                                           cTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGaa                  <  1:24131/124‑1 (MQ=255)                                                                                                                           cTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCg               >  2:306665/1‑127 (MQ=255)                                                                                                                           cTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCg               <  1:306665/127‑1 (MQ=255)                                                                                                                           cTTTCATCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGc  >  1:588670/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                      CTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGC  >  NC_000913/3995784‑3996045

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