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

CGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTA  >  NC_000913/3385729‑3385984                                                                                                                                 |                                                                                                                                 cgcAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGATGaa                                                                                                                        <  2:696755/139‑1 (MQ=255)   cAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACa                                                                                                                      >  1:863518/1‑139 (MQ=255)                 aTTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATg                                                                                                       <  1:671494/139‑1 (MQ=255)                   tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACAt                                                                                                                     >  1:323099/1‑124 (MQ=255)                   tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACAt                                                                                                                     <  2:323099/124‑1 (MQ=255)                       gAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTc                                                                                                 >  2:532607/1‑139 (MQ=255)                           cTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGc                                                                                             >  1:917967/1‑139 (MQ=255)                               tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCtt                                                                                         >  2:66741/1‑139 (MQ=255)                                     tGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTg                                                                                   <  1:66741/139‑1 (MQ=255)                                     tGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTg                                                                                   <  2:863518/139‑1 (MQ=255)                                       acacCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                                 >  2:812971/1‑139 (MQ=255)                                              gCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaa                                                                          >  1:184566/1‑139 (MQ=255)                                               cACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaaa                                                                         <  2:184566/139‑1 (MQ=255)                                                aCGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                        <  1:888353/139‑1 (MQ=255)                                                aCGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                        <  2:642656/139‑1 (MQ=255)                                                   cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGc                                                                                                      >  1:855394/1‑106 (MQ=255)                                                   cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGc                                                                                                      <  2:855394/106‑1 (MQ=255)                                                            aaCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAAc                                                            >  1:550556/1‑139 (MQ=255)                                                                gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACCGCAAATGCCCTTCCAGct                                                                                            >  2:199762/1‑104 (MQ=255)                                                                gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACCGCAAATGCCCTTCCAGct                                                                                            <  1:199762/104‑1 (MQ=255)                                                                 tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTcc                                                                                                <  1:1060644/98‑1 (MQ=255)                                                                 tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTcc                                                                                                >  2:1060644/1‑98 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaaa                                                                                                                >  2:357689/1‑73 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaaa                                                                                                                <  1:357689/73‑1 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                             >  1:648450/1‑139 (MQ=255)                                                                                aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                                 >  2:420534/1‑98 (MQ=255)                                                                                aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                                 <  1:420534/98‑1 (MQ=255)                                                                                   gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                                     <  1:479599/139‑1 (MQ=255)                                                                                    tCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtc                                    >  2:724333/1‑139 (MQ=255)                                                                                     cTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtct                                   >  2:156020/1‑139 (MQ=255)                                                                                           cgacgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCa                             <  1:536322/139‑1 (MQ=255)                                                                                             ccgaTATTGCCCCCGTTTGGGGGCTTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGa                           <  2:183605/138‑1 (MQ=255)                                                                                             acgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGa                           <  2:550556/139‑1 (MQ=255)                                                                                              cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGAt                          >  1:809009/1‑139 (MQ=255)                                                                                              cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCCTTTATGACGAAACTTTAAA‑CGCAAATGCCCTTCCAGCCCTTCCTCTGCCTCATCAAACAAGAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGAt                          >  1:183605/1‑139 (MQ=255)                                                                                                   ttGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGa                           <  2:917996/133‑1 (MQ=255)                                                                                                   ttGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGa                           >  1:917996/1‑133 (MQ=255)                                                                                                              tCGGGGGC‑TTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCa         <  1:765765/139‑1 (MQ=255)                                                                                                                gggggCTTTTTTTTGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCaa        >  1:954755/1‑139 (MQ=255)                                                                                                                      ttttttttGCGTTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTa  >  1:171497/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                 CGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAA‑CGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTA  >  NC_000913/3385729‑3385984

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