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%	139.1 / NA	44	*91* (TAG→TAA)	yhcO	putative barnase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/22);  total (22/22)

GAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAA  >  NC_000913/3385723‑3385978                                                                                                                                       |                                                                                                                          gaAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAAt                                                                                                                       >  1:289596/1‑139 (MQ=255)              gCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaa                                                                                                          <  2:1228324/139‑1 (MQ=255)                       aTTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                 <  2:289596/139‑1 (MQ=255)                       aTTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                 <  1:1073873/139‑1 (MQ=255)                         tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGcc                                                                                               >  2:962101/1‑139 (MQ=255)                             gAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTc                                                                                           <  1:962101/139‑1 (MQ=255)                                gCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGATCGGGGGCTTTTTTTTGCGTTAATGACGAAAATTAAAACGCAAATGCCCTTCCac                                                                                        >  2:1072821/1‑138 (MQ=255)                                 cTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGc                                                                                       >  2:149213/1‑139 (MQ=255)                                   cGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctc                                                                                     >  1:1184931/1‑139 (MQ=255)                                   cGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctc                                                                                     >  1:920508/1‑139 (MQ=255)                                          gTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCctct                                                                              >  1:290907/1‑139 (MQ=255)                                                        gCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg                                                                <  2:1184931/139‑1 (MQ=255)                                                         cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTC‑GGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCa                                                              >  2:403581/1‑139 (MQ=255)                                                                      gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                           <  1:374168/114‑1 (MQ=255)                                                                      gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                           >  2:374168/1‑114 (MQ=255)                                                                       tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                                    >  2:265749/1‑104 (MQ=255)                                                                       tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                                    <  1:265749/104‑1 (MQ=255)                                                                              aaaCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaaa                                                                   >  2:620134/1‑114 (MQ=255)                                                                              aaaCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaaa                                                                   <  1:620134/114‑1 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                 >  2:595772/1‑65 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                 <  1:595772/65‑1 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCaaa                                                                                                   <  2:655990/79‑1 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCaaa                                                                                                   >  1:655990/1‑79 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                 <  1:826122/81‑1 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                 >  2:826122/1‑81 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGc                                                                                                <  2:548466/82‑1 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGc                                                                                                >  1:548466/1‑82 (MQ=255)                                                                                 cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                       <  1:149213/139‑1 (MQ=255)                                                                                      aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAACTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGt                                  >  2:501086/1‑139 (MQ=255)                                                                                         gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                               <  2:920508/139‑1 (MQ=255)                                                                                                   acgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGa                     >  2:219916/1‑139 (MQ=255)                                                                                                    cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTc                                                                                           >  2:21710/1‑68 (MQ=255)                                                                                                    cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTc                                                                                           <  1:21710/68‑1 (MQ=255)                                                                                                    cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTc                                                                                  <  2:156779/77‑1 (MQ=255)                                                                                                    cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTc                                                                                  >  1:156779/1‑77 (MQ=255)                                                                                                     gaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATg                   <  1:677426/139‑1 (MQ=255)                                                                                                               ccGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCa         <  1:219916/139‑1 (MQ=255)                                                                                                                 ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCtcatca                                                                     <  1:95333/77‑1 (MQ=255)                                                                                                                 ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCtcatca                                                                     >  2:95333/1‑77 (MQ=255)                                                                                                                      gggggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCa                                                                                         <  1:1200107/52‑1 (MQ=255)                                                                                                                      gggggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCa                                                                                         >  2:1200107/1‑52 (MQ=255)                                                                                                                      gggggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCaa  <  2:1057695/139‑1 (MQ=255)                                                                                                                             tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACa                                                                 >  1:864648/1‑69 (MQ=255)                                                                                                                             tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACa                                                                 <  2:864648/69‑1 (MQ=255)                                                                                                                                       |                                                                                                                          GAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAA  >  NC_000913/3385723‑3385978

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