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

CGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTA  >  NC_000913/3385729‑3385984                                                                                                                                 |                                                                                                                                cgcAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                       >  2:46092/1‑139 (MQ=255)       ggCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaa                                                                                                                 <  1:46092/139‑1 (MQ=255)                   tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCa                                                                                                           >  2:750651/1‑133 (MQ=255)                   tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCa                                                                                                           <  1:750651/133‑1 (MQ=255)                             cGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctc                                                                                           >  2:455483/1‑139 (MQ=255)                               tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                >  2:330832/1‑132 (MQ=255)                               tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                <  1:330832/132‑1 (MQ=255)                               tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCtt                                                                                         >  1:1302020/1‑139 (MQ=255)                                     tGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaaa                                                                                                               >  1:1359617/1‑111 (MQ=255)                                     tGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaaa                                                                                                               <  2:1359617/111‑1 (MQ=255)                                           cTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCtcat                                                                             >  1:1481025/1‑139 (MQ=255)                                             ggCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCtcatca                                                                           <  1:455483/139‑1 (MQ=255)                                                  gCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg                                                                      <  1:51608/139‑1 (MQ=255)                                                   cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACg                                                                                                                         <  2:1637950/87‑1 (MQ=255)                                                   cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACg                                                                                                                         >  1:1637950/1‑87 (MQ=255)                                                   cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGc                                                                     >  1:718639/1‑139 (MQ=255)                                                            aaCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                        >  1:113565/1‑127 (MQ=255)                                                            aaCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                        <  2:113565/127‑1 (MQ=255)                                                             aCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAAcg                                                           <  1:1545849/139‑1 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                       >  1:1212111/1‑65 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                       <  2:1212111/65‑1 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCct                                                                                      >  1:178335/1‑98 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCct                                                                                      <  2:178335/98‑1 (MQ=255)                                                                           cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                             >  2:54789/1‑139 (MQ=255)                                                                               cATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGAcgcg                                         >  2:453055/1‑139 (MQ=255)                                                                                   gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                                     <  2:718639/139‑1 (MQ=255)                                                                                   gTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                                     <  2:1481025/139‑1 (MQ=255)                                                                                     cTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTcc                                                                                       <  1:1249967/87‑1 (MQ=255)                                                                                     cTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTcc                                                                                       >  2:1249967/1‑87 (MQ=255)                                                                                         tCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTcc                               <  1:453055/139‑1 (MQ=255)                                                                                              cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtctct                                 >  1:1437706/1‑132 (MQ=255)                                                                                              cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtctct                                 <  2:1437706/132‑1 (MQ=255)                                                                                              cgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTCCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGAt                          >  2:1143074/1‑139 (MQ=255)                                                                                                  atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg                                                                      >  1:1110001/1‑91 (MQ=255)                                                                                                  atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAg                                                                      <  2:1110001/91‑1 (MQ=255)                                                                                                   ttGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACa                     >  1:1593357/1‑139 (MQ=255)                                                                                                          cGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGc                                                                                                            >  1:375811/1‑45 (MQ=255)                                                                                                          cGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGc                                                                                                            <  2:375811/45‑1 (MQ=255)                                                                                                           ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                <  2:273724/56‑1 (MQ=255)                                                                                                           ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                >  1:273724/1‑56 (MQ=255)                                                                                                           ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGct                                                                                            >  1:1380501/1‑60 (MQ=255)                                                                                                           ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGct                                                                                            <  2:1380501/60‑1 (MQ=255)                                                                                                             ttCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGc                                                                                             <  2:712407/57‑1 (MQ=255)                                                                                                             ttCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGc                                                                                             >  1:712407/1‑57 (MQ=255)                                                                                                               cGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                <  2:1657312/52‑1 (MQ=255)                                                                                                               cGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                >  1:1657312/1‑52 (MQ=255)                                                                                                                gggggCTTTTTTTTGCGTTAATGACGAACATTaaaa                                                                                                               <  2:503504/36‑1 (MQ=255)                                                                                                                gggggCTTTTTTTTGCGTTAATGACGAACATTaaaa                                                                                                               >  1:503504/1‑36 (MQ=255)                                                                                                                     cTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGt   >  1:1649405/1‑139 (MQ=255)                                                                                                                      ttttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTa  >  1:236339/1‑139 (MQ=255)                                                                                                                       tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCtt                                                                                                  >  2:1082765/1‑42 (MQ=255)                                                                                                                       tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCtt                                                                                                  <  1:1082765/42‑1 (MQ=255)                                                                                                                       tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                        <  1:1311338/68‑1 (MQ=255)                                                                                                                       tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                        >  2:1311338/1‑68 (MQ=255)                                                                                                                                 |                                                                                                                                CGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTA  >  NC_000913/3385729‑3385984

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