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

CAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCA  >  NC_000913/3385731‑3385990                                                                                                                               |                                                                                                                                      caatggCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACa                                                                                                                           >  1:764257/5‑139 (MQ=255)  aaCGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACAt                                                                                                                          >  2:1159961/1‑139 (MQ=255)                     gAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTc                                                                                                       >  1:882576/1‑139 (MQ=255)                     gAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTc                                                                                                       >  1:377740/1‑139 (MQ=255)                      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                             <  2:206662/116‑1 (MQ=255)                      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                             >  1:206662/1‑116 (MQ=255)                      aaGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                      <  1:217446/139‑1 (MQ=255)                              aGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTc                                                                                              <  1:335679/139‑1 (MQ=255)                                            gCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaa                                                                                <  2:1258966/139‑1 (MQ=255)                                              aCGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                              <  2:882576/139‑1 (MQ=255)                                                 cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaa                                                                                                                       <  2:871187/95‑1 (MQ=255)                                                 cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaa                                                                                                                       >  1:871187/1‑95 (MQ=255)                                                  tACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGc                                                                                                            >  1:468655/1‑105 (MQ=255)                                                  tACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGc                                                                                                            <  2:468655/105‑1 (MQ=255)                                                              gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                                       <  2:704196/114‑1 (MQ=255)                                                              gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                                       >  1:704196/1‑114 (MQ=255)                                                               tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGcc                                                             >  1:448820/1‑139 (MQ=255)                                                                         cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             <  2:1002422/81‑1 (MQ=255)                                                                         cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             >  1:1002422/1‑81 (MQ=255)                                                                         cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                   >  2:771375/1‑139 (MQ=255)                                                                             cATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             >  2:443965/1‑77 (MQ=255)                                                                             cATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             <  1:443965/77‑1 (MQ=255)                                                                              aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             >  1:1149063/1‑76 (MQ=255)                                                                              aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             <  2:1149063/76‑1 (MQ=255)                                                                              aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCt                                                                                      >  1:600164/1‑99 (MQ=255)                                                                              aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCt                                                                                      <  2:600164/99‑1 (MQ=255)                                                                                   cTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                             >  1:1043742/1‑55 (MQ=255)                                                                                   cTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                             <  2:1043742/55‑1 (MQ=255)                                                                                       tCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                              <  1:607237/98‑1 (MQ=255)                                                                                       tCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                              >  2:607237/1‑98 (MQ=255)                                                                                                atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTACCTCATCAAACAGc                                                                           <  2:413610/92‑1 (MQ=255)                                                                                                atTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTACCTCATCAAACAGc                                                                           >  1:413610/1‑92 (MQ=255)                                                                                                    cccccGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             <  1:603369/54‑1 (MQ=255)                                                                                                    cccccGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                             >  2:603369/1‑54 (MQ=255)                                                                                                         ggTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAAt                   <  1:1159961/139‑1 (MQ=255)                                                                                                          gTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctc                                                                                                 >  1:595577/1‑60 (MQ=255)                                                                                                          gTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctc                                                                                                 <  2:595577/60‑1 (MQ=255)                                                                                                                gggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGc                                                     <  1:939829/98‑1 (MQ=255)                                                                                                                gggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGc                                                     >  2:939829/1‑98 (MQ=255)                                                                                                                   cTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                   >  2:407221/1‑97 (MQ=255)                                                                                                                   cTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                   <  1:407221/97‑1 (MQ=255)                                                                                                                     tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                   <  2:941152/95‑1 (MQ=255)                                                                                                                     tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                   >  1:941152/1‑95 (MQ=255)                                                                                                                     tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAg       >  1:283284/1‑139 (MQ=255)                                                                                                                          ttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCa  <  1:470097/139‑1 (MQ=255)                                                                                                                            gcgtTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                                           <  2:569287/96‑1 (MQ=255)                                                                                                                            gcgtTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGtttt                                           >  1:569287/1‑96 (MQ=255)                                                                                                                               |                                                                                                                                      CAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCA  >  NC_000913/3385731‑3385990

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