Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,385,857	C→T	*91* (TAG→TAA)	yhcO ←	putative barnase inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,385,857	0	C	T	100.0%	115.1 / NA	37	*91* (TAG→TAA)	yhcO	putative barnase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/18);  total (19/18)

TACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCAG  >  NC_000913/3385747‑3385991                                                                                                               |                                                                                                                                       tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGAc                                                                                                                                 <  2:986824/118‑1 (MQ=255) tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGAc                                                                                                                                 >  1:986824/1‑118 (MQ=255)   cTGAATCTCGTAGCGGTGACACCTGGCACGCTAGGGCTGACCTGTACAAATAAACCCTCATCGTCTTGCCCGACGATATTGCCCCCGGATCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCt                                                                                                          <  2:557177/139‑1 (MQ=255)      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATCTTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGAc                                                                                                                                 >  1:1075397/1‑113 (MQ=255)      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGAc                                                                                                                                 <  2:1075397/113‑1 (MQ=255)      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCa                                                                                                                  >  2:1112004/1‑128 (MQ=255)      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCa                                                                                                                  <  1:1112004/128‑1 (MQ=255)      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                       >  2:412269/1‑139 (MQ=255)      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                       <  1:201573/139‑1 (MQ=255)              aGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCATCCAGCTCTTc                                                                                               >  1:721256/1‑139 (MQ=255)                                gCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg                                                                             <  1:412269/139‑1 (MQ=255)                                 cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACg                                                                                                                                >  1:632745/1‑87 (MQ=255)                                 cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACg                                                                                                                                <  2:632745/87‑1 (MQ=255)                                 cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGc                                                                            <  1:223461/139‑1 (MQ=255)                                         gAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg                                                                             >  2:981927/1‑130 (MQ=255)                                         gAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg                                                                             <  1:981927/130‑1 (MQ=255)                                          aaCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                       >  2:52991/1‑103 (MQ=255)                                          aaCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                       <  1:52991/103‑1 (MQ=255)                                            cTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaa                                                                                                                       <  1:854612/85‑1 (MQ=255)                                            cTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaa                                                                                                                       >  2:854612/1‑85 (MQ=255)                                               tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCACCCGGTTCGGAGGCTTTTTTTTGCGTTAATAACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAAAGCGcc                                                              >  2:1032189/1‑139 (MQ=255)                                                      aaaCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACgg                                                       >  2:408763/1‑139 (MQ=255)                                                         cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCa                                                                                                                  >  1:280379/1‑77 (MQ=255)                                                         cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCa                                                                                                                  <  2:280379/77‑1 (MQ=255)                                                         cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                    <  1:408763/139‑1 (MQ=255)                                                         cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                    >  1:27065/1‑139 (MQ=255)                                                                       tCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTcc                                      >  1:634070/1‑139 (MQ=255)                                                                                atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGc                                                      <  2:998627/114‑1 (MQ=255)                                                                                atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGc                                                      >  1:998627/1‑114 (MQ=255)                                                                                         ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAg                                                                                                     <  2:1028945/58‑1 (MQ=255)                                                                                         ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAg                                                                                                     >  1:1028945/1‑58 (MQ=255)                                                                                         ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTc                                                                                               >  2:43368/1‑64 (MQ=255)                                                                                         ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTc                                                                                               <  1:43368/64‑1 (MQ=255)                                                                                                     tttttttGCGTTAATGACGAACATTAAAACGCAAATGc                                                                                                             >  2:377486/1‑38 (MQ=255)                                                                                                     tttttttGCGTTAATGACGAACATTAAAACGCAAATGc                                                                                                             <  1:377486/38‑1 (MQ=255)                                                                                                     tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAg        >  2:214983/1‑139 (MQ=255)                                                                                                           tGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCAg  <  1:214983/139‑1 (MQ=255)                                                                                                               |                                                                                                                                       TACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCAG  >  NC_000913/3385747‑3385991

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