Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,197,963	A→G	intergenic (+665/‑2316)	yehF → / → yehI	putative protein YehF/DUF4132 domain‑containing protein YehI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,197,963	0	A	G	100.0%	106.6 / NA	30	intergenic (+665/‑2316)	yehF/yehI	putative protein YehF/DUF4132 domain‑containing protein YehI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/14);  total (16/14)

AGTTGGCGCATTATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAA  >  NC_000913/2197833‑2198096                                                                                                                                   |                                                                                                                                      aGTTGGCGCATTATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAg                                                                                                                              <  1:36675/140‑1 (MQ=255) aGTTGGCGCATTATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAg                                                                                                                              <  1:403385/140‑1 (MQ=255)   ttGGCGCATTATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCGGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAgcg                                                                                                                            <  2:268239/140‑1 (MQ=255)          aTTATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCt                                                                                                                     >  1:417060/1‑140 (MQ=255)          aTTATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCt                                                                                                                     >  1:338071/1‑140 (MQ=255)           ttATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTa                                                                                                                    >  1:318519/1‑140 (MQ=255)             aTGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACt                                                                                                                  >  2:291725/1‑140 (MQ=255)                   aTGGCAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGc                                                                                                            >  2:151117/1‑140 (MQ=255)                               aTCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCgta                                                                                                        >  1:86257/1‑132 (MQ=255)                               aTCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCgta                                                                                                        <  2:86257/132‑1 (MQ=255)                                tctcACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGTCAGCCGTAGTAGTctcc                                                                                               <  1:81740/140‑1 (MQ=255)                                    acaACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTctcccctc                                                                                           >  2:246063/1‑140 (MQ=255)                                    acaACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTctcccctc                                                                                           >  2:333088/1‑140 (MQ=255)                                           gCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGtctc                                                                                                >  1:129516/1‑128 (MQ=255)                                           gCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGtctc                                                                                                <  2:129516/128‑1 (MQ=255)                                                                   tGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGAt                                                            >  1:279076/1‑140 (MQ=255)                                                                    ggCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGGTCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTGGAAAAAGAAAAAATCGCCGAtt                                                           >  2:528428/1‑140 (MQ=255)                                                                     gCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTc                                                          <  2:417060/140‑1 (MQ=255)                                                                     gCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTc                                                          <  1:210333/140‑1 (MQ=255)                                                                     gCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTc                                                          <  2:318519/140‑1 (MQ=255)                                                                       tctcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCg                                                        <  2:394003/140‑1 (MQ=255)                                                                        ctcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCgg                                                       >  2:358961/1‑140 (MQ=255)                                                                         tcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGt                                                      >  1:62955/1‑140 (MQ=255)                                                                                   cAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTa                                            >  2:317833/1‑140 (MQ=255)                                                                                          cGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAg                                           >  1:232155/1‑134 (MQ=255)                                                                                          cGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAg                                           <  2:232155/134‑1 (MQ=255)                                                                                           gTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCAtt                                    <  1:528428/140‑1 (MQ=255)                                                                                                     cATGCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAg                          <  1:333088/140‑1 (MQ=255)                                                                                                        gCCAGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCa                       <  2:338071/140‑1 (MQ=255)                                                                                                           aGCAACAACTGACACAGCCCTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATc                    >  1:209737/1‑140 (MQ=255)                                                                                                                             ccTCAAGCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAgaagaa  >  1:447530/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                      AGTTGGCGCATTATGCTAATGACAATGCTTATCTCACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAA  >  NC_000913/2197833‑2198096

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