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%	179.7 / NA	56	*112* (TAG→TAA)	mazF	mRNA interferase toxin, antitoxin is MazE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/26);  total (30/26)

CGGGCAGCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAT  >  NC_000913/2910620‑2910883                                                                                                                                         |                                                                                                                                cGGGCAGCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTa                                                                                                                               <  1:1564828/139‑1 (MQ=255)  gggCAGCCGTTTTCCGGATCGCGCAGGCGCTGTATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAc                                                                                                                              >  2:665353/1‑139 (MQ=255)        ccGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATc                                                                                                                        <  1:1626661/139‑1 (MQ=255)        ccGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATc                                                                                                                        >  2:1626661/1‑139 (MQ=255)                                        aGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                        <  1:665353/139‑1 (MQ=255)                                            aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCa                                                                                                                           <  1:361343/100‑1 (MQ=255)                                            aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCa                                                                                                                           >  2:361343/1‑100 (MQ=255)                                            aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                    >  2:1622725/1‑139 (MQ=255)                                            aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                    >  2:175227/1‑139 (MQ=255)                                            aGCAAACGGACGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                    >  2:1120149/1‑139 (MQ=255)                                                             aTTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTcc                                                                   >  2:1391156/1‑139 (MQ=255)                                                                        tGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTGCCCTTCTTCGTTg                                                        <  1:1120149/139‑1 (MQ=255)                                                                         gTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGc                                                       <  1:175227/139‑1 (MQ=255)                                                                          tCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAc                                                                                                                              >  2:857287/1‑67 (MQ=255)                                                                          tCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAc                                                                                                                              <  1:857287/67‑1 (MQ=255)                                                                            cTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGctcc                                                    >  1:1153013/1‑139 (MQ=255)                                                                            cTGAAAATTGCGGGTCTGTAAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGctcc                                                    >  1:1157247/1‑139 (MQ=255)                                                                               aaaaTTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGTATAGATGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCctct                                                 >  2:1118201/1‑139 (MQ=255)                                                                               aaaaTTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTACTTTCTTCGTTGCtcttct                                                 >  1:1177179/1‑139 (MQ=255)                                                                                    tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCaa                                                                          <  2:705496/109‑1 (MQ=255)                                                                                    tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCaa                                                                          >  1:705496/1‑109 (MQ=255)                                                                                    tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTttctt                                                             <  2:116487/122‑1 (MQ=255)                                                                                    tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTttctt                                                             >  1:116487/1‑122 (MQ=255)                                                                                      cGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGt                                                                                                                     >  2:284265/1‑64 (MQ=255)                                                                                      cGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGt                                                                                                                     <  1:284265/64‑1 (MQ=255)                                                                                             gTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                   <  1:1622725/139‑1 (MQ=255)                                                                                              tCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCaa                                                                          >  1:1281833/1‑99 (MQ=255)                                                                                              tCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCaa                                                                          <  2:1281833/99‑1 (MQ=255)                                                                                                   tGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTa                                                                                                   >  2:596253/1‑69 (MQ=255)                                                                                                   tGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTa                                                                                                   <  1:596253/69‑1 (MQ=255)                                                                                                   tGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGt                                                                                          <  1:1262476/78‑1 (MQ=255)                                                                                                   tGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGt                                                                                          >  2:1262476/1‑78 (MQ=255)                                                                                                     gAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTg                                                        >  2:605565/1‑110 (MQ=255)                                                                                                     gAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTg                                                        <  1:605565/110‑1 (MQ=255)                                                                                                      aaaCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGa                                                                                               >  1:659239/1‑70 (MQ=255)                                                                                                      aaaCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGa                                                                                               <  2:659239/70‑1 (MQ=255)                                                                                                        aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACg                                                                                                                  <  2:358654/49‑1 (MQ=255)                                                                                                        aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACg                                                                                                                  >  1:358654/1‑49 (MQ=255)                                                                                                        aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAcc                        >  2:798389/1‑139 (MQ=255)                                                                                                           tgtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCaa                                                                          <  1:178927/86‑1 (MQ=255)                                                                                                           tgtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCaa                                                                          >  2:178927/1‑86 (MQ=255)                                                                                                           tgtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGa                     <  2:1177179/139‑1 (MQ=255)                                                                                                            gtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGTTACTTTTTACCTGAt                    <  1:1118201/139‑1 (MQ=255)                                                                                                             tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATCCtt                              <  1:147326/128‑1 (MQ=255)                                                                                                             tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACtt                              >  2:147326/1‑128 (MQ=255)                                                                                                             tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATc                   >  1:665157/1‑139 (MQ=255)                                                                                                                ccAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGc                <  2:1157247/139‑1 (MQ=255)                                                                                                                    aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGtt                                                                                            <  1:910040/59‑1 (MQ=255)                                                                                                                    aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGtt                                                                                            >  2:910040/1‑59 (MQ=255)                                                                                                                        gAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTa        <  1:798389/139‑1 (MQ=255)                                                                                                                              aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTGCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTg                                               <  1:248321/94‑1 (MQ=255)                                                                                                                              aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTg                                               >  2:248321/1‑94 (MQ=255)                                                                                                                              aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                   <  2:1165017/106‑1 (MQ=255)                                                                                                                              aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                   >  1:1165017/1‑106 (MQ=255)                                                                                                                              aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt  >  2:1602934/1‑139 (MQ=255)                                                                                                                              aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt  >  2:1540247/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                CGGGCAGCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAT  >  NC_000913/2910620‑2910883

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