Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,385,853	(T)8→9	intergenic (+52/+4)	yhcN → / ← yhcO	DUF1471 domain‑containing stress‑induced protein YhcN/putative barnase inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,385,845	1	.	T	100.0%	80.4 / NA	23	intergenic (+44/+12)	yhcN/yhcO	DUF1471 domain‑containing stress‑induced protein YhcN/putative barnase inhibitor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (11/12);  total (11/12)

GACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑TTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAA  >  NC_000913/3385766‑3385919                                                                                 |                                                                           gACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑t                                                                           >  1:576636/1‑81 (MQ=255)         gCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑ttttttttt                                                                   >  1:576643/1‑80 (MQ=255)         gCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑ttttttt                                                                     <  2:576595/79‑1 (MQ=255)           aCGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑tttttttttgc                                                                 >  1:576631/1‑78 (MQ=255)           aCGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑tttttttttgc                                                                 <  2:576596/81‑4 (MQ=255)            cGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑tttttttttgcg                                                                >  2:576633/1‑77 (MQ=255)                  ggCTGAACTGTACAATTAAACCCTCACCGTCTTTTCCGACGATATTGCCTCCGGTTCGGGGGCTTTTTTTTTGCGCTAATg                                                           >  2:576643/1‑81 (MQ=255)                     tGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATTAc                                                         >  2:576634/1‑80 (MQ=255)                       agcTGTACAAGTAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCAGTTCGGGGGCTTTTTTTTTGCGCTGATGACGaa                                                      <  2:576597/79‑1 (MQ=255)                            tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTa                                                 <  1:576599/81‑1 (MQ=255)                              cAAATAAACCCTCATCGTCTTGTACGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACGTTaa                                                <  2:576599/80‑1 (MQ=255)                               aaaTAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCTCTAATGACGAACATTaaaa                                              >  2:576631/1‑81 (MQ=255)                                   aaaCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGc                                           >  2:576639/1‑80 (MQ=255)                                    aaCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCaa                                         <  1:576601/81‑1 (MQ=255)                                           aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCt                                  >  1:576627/1‑81 (MQ=255)                                             cGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTc                                >  1:576633/1‑81 (MQ=255)                                              gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTcc                               <  2:576602/81‑1 (MQ=255)                                              gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTcc                               <  1:576602/81‑1 (MQ=255)                                               tCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCa                              <  2:576603/81‑1 (MQ=255)                                                 ttGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGc                            >  1:576629/1‑81 (MQ=255)                                                       gacgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTcc                      >  1:576646/1‑81 (MQ=255)                                                         cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCctc                    <  2:576604/81‑1 (MQ=255)                                                          gaTATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCctct                   <  1:576606/81‑1 (MQ=255)                                                              ttGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGcc                <  1:576608/80‑1 (MQ=255)                                                              ttGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCt               <  2:576606/81‑1 (MQ=255)                                                                       gTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaaa        >  2:576644/1‑79 (MQ=255)                                                                         tCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg     >  2:576638/1‑80 (MQ=255)                                                                           gggggCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCaa  >  2:576656/1‑81 (MQ=255)                                                                           gggggCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCa   <  2:576608/80‑1 (MQ=255)                                                                                 |                                                                           GACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGC‑TTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAA  >  NC_000913/3385766‑3385919

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