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

GGGCAGCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCA  >  NC_000913/2910621‑2910885                                                                                                                                        |                                                                                                                                  gggCAGCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAc                                                                                                                                <  2:987918/139‑1 (MQ=255) gggCAGCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAc                                                                                                                                >  2:77571/1‑139 (MQ=255)                    cgcgCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAAttt                                                                                                             <  1:772221/139‑1 (MQ=255)                                         tGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                          <  2:309243/137‑1 (MQ=255)                                         tGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                          >  1:309243/1‑137 (MQ=255)                                         tGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAAtt                                                                                        >  2:742106/1‑139 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                  <  2:52407/127‑1 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                  >  1:979238/1‑127 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                  >  1:52407/1‑127 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                  <  2:979238/127‑1 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGt                                                                                            <  1:242085/133‑1 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGt                                                                                            >  2:242085/1‑133 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                      >  2:907786/1‑139 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                      >  2:419369/1‑139 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                      >  2:309036/1‑139 (MQ=255)                                           aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                      >  1:558525/1‑139 (MQ=255)                                                      gATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAAc                                                                           >  2:516201/1‑139 (MQ=255)                                                              tCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTa                                                                                                                                 <  1:642610/77‑1 (MQ=255)                                                                         tCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGct                                                        >  1:921616/1‑139 (MQ=255)                                                                          ccTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTggg                                                                                >  1:507143/1‑114 (MQ=255)                                                                          ccTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTggg                                                                                <  2:507143/114‑1 (MQ=255)                                                                                   tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAAt                                                                                                   >  1:744757/1‑86 (MQ=255)                                                                                   tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAAt                                                                                                   <  2:744757/86‑1 (MQ=255)                                                                                   tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGccc                                              >  2:504770/1‑139 (MQ=255)                                                                                        gTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCa                                                                                                                             <  1:36308/55‑1 (MQ=255)                                                                                        gTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCa                                                                                                                             >  2:36308/1‑55 (MQ=255)                                                                                           tGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTg                                                                         >  2:432293/1‑104 (MQ=255)                                                                                           tGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTg                                                                         <  1:432293/104‑1 (MQ=255)                                                                                                   ggAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACtttt                              >  2:375004/1‑139 (MQ=255)                                                                                                       aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGt                                                            >  2:486171/1‑105 (MQ=255)                                                                                                       aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGt                                                            <  1:486171/105‑1 (MQ=255)                                                                                                       aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                     >  2:186582/1‑128 (MQ=255)                                                                                                       aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                     <  1:186582/128‑1 (MQ=255)                                                                                                       aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAcc                          >  2:204418/1‑139 (MQ=255)                                                                                                          tgtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGa                       <  1:375004/139‑1 (MQ=255)                                                                                                           gtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGAt                      >  1:848261/1‑139 (MQ=255)                                                                                                            tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATc                     >  2:607952/1‑139 (MQ=255)                                                                                                                   aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAAttt                                                                                                             <  2:446990/44‑1 (MQ=255)                                                                                                                   aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAAttt                                                                                                             >  1:446990/1‑44 (MQ=255)                                                                                                                        aaGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                     <  2:341294/111‑1 (MQ=255)                                                                                                                        aaGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGa                                     >  1:341294/1‑111 (MQ=255)                                                                                                                          gTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGc       <  1:419369/139‑1 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGa                                                                                                 <  2:369902/46‑1 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGa                                                                                                 >  1:369902/1‑46 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTgg                                                                                 >  2:999971/1‑62 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTgg                                                                                 <  1:999971/62‑1 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAgg                                        <  1:819296/103‑1 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAgg                                        >  2:819296/1‑103 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  2:383773/1‑139 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  2:713176/1‑139 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  1:855145/1‑139 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  1:301860/1‑139 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  1:129385/1‑139 (MQ=255)                                                                                                                             aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  1:108504/1‑139 (MQ=255)                                                                                                                               ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa  <  1:204418/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  GGGCAGCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCA  >  NC_000913/2910621‑2910885

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