Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,910,756	C→T	*112* (TAG→TAA)	mazF ←	mRNA interferase toxin, antitoxin is MazE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,910,756	0	C	T	100.0%	189.4 / NA	58	*112* (TAG→TAA)	mazF	mRNA interferase toxin, antitoxin is MazE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/26);  total (32/26)

ATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGT  >  NC_000913/2910657‑2910888                                                                                                    |                                                                                                                                     ataGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGt                                                                                               >  2:302830/1‑139 (MQ=255) ataGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGt                                                                                               >  1:299490/1‑139 (MQ=255)     tGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAAtt                                                                                           >  1:309001/1‑139 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAc                                                                                                                                   <  2:420809/97‑1 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAc                                                                                                                                   >  1:420809/1‑97 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAAt                                                                                                                              <  2:227683/102‑1 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAAt                                                                                                                              >  1:227683/1‑102 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTa                                                                                                        >  2:697195/1‑124 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTa                                                                                                        <  1:697195/124‑1 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAAt                                                                                                      <  2:93139/126‑1 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAAt                                                                                                      >  1:93139/1‑126 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                         >  2:496799/1‑139 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                         >  1:67295/1‑139 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                         >  2:771431/1‑139 (MQ=255)       aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                         >  2:778380/1‑139 (MQ=255)         cAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGt                                                                                                  <  1:450357/128‑1 (MQ=255)         cAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGt                                                                                                  >  2:450357/1‑128 (MQ=255)            aCGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACTCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                     <  2:408723/122‑1 (MQ=255)            aCGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACTCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                     >  1:408723/1‑122 (MQ=255)               gTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGc                                                                                 <  2:309001/139‑1 (MQ=255)                                          aaaaTTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCctct                                                      >  2:85648/1‑139 (MQ=255)                                                gCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCg                                                <  1:771431/139‑1 (MQ=255)                                                 cGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGc                                               <  1:149946/139‑1 (MQ=255)                                                  gggTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGcc                                              <  1:996183/139‑1 (MQ=255)                                                   ggTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCa                                             <  1:365441/139‑1 (MQ=255)                                                        gTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                        <  1:302830/139‑1 (MQ=255)                                                         tCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                        <  1:191402/138‑1 (MQ=255)                                                         tCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                        >  2:191402/1‑138 (MQ=255)                                                                   aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAACCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAcc                             >  2:361492/1‑139 (MQ=255)                                                                       gtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGAt                         >  2:1210512/1‑139 (MQ=255)                                                                        tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATc                        >  1:1129030/1‑139 (MQ=255)                                                                        tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATc                        >  1:679587/1‑139 (MQ=255)                                                                            cAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCt                    <  1:85648/139‑1 (MQ=255)                                                                            cAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCt                    <  1:639282/139‑1 (MQ=255)                                                                               aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAAc                 >  2:194770/1‑139 (MQ=255)                                                                               aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAAc                 >  1:511793/1‑139 (MQ=255)                                                                               aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAAc                 >  2:516924/1‑139 (MQ=255)                                                                               aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAAc                 >  2:538157/1‑139 (MQ=255)                                                                                   gAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTa             <  1:1210512/139‑1 (MQ=255)                                                                                      gTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGc          <  2:299490/139‑1 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                     >  1:16042/1‑45 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                     <  2:16042/45‑1 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCg                                                                >  2:1093929/1‑82 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCg                                                                <  1:1093929/82‑1 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGt                                                               <  2:572796/83‑1 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGt                                                               >  1:572796/1‑83 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                        >  1:645623/1‑106 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                        <  2:645623/106‑1 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGa                          >  1:127166/1‑120 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGa                          <  2:127166/120‑1 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt       >  1:578540/1‑139 (MQ=255)                                                                                         aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt       >  1:937518/1‑139 (MQ=255)                                                                                           ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa     <  2:937518/139‑1 (MQ=255)                                                                                           ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa     <  2:679587/139‑1 (MQ=255)                                                                                           ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa     <  2:578540/139‑1 (MQ=255)                                                                                           ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa     >  2:1141562/1‑139 (MQ=255)                                                                                           ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa     <  2:1129030/139‑1 (MQ=255)                                                                                              gTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGCTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGt  >  2:751286/1‑139 (MQ=255)                                                                                                    |                                                                                                                                     ATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGT  >  NC_000913/2910657‑2910888

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