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

CGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGTTCC  >  NC_000913/2910642‑2910891                                                                                                                   |                                                                                                                                        cgcAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAg                                                                                                                              >  1:373454/1‑126 (MQ=255) cgcAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAg                                                                                                                              <  2:373454/126‑1 (MQ=255)          tGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                        <  1:441870/139‑1 (MQ=255)                    tGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGt                                                                                                                         >  2:897145/1‑112 (MQ=255)                    tGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGt                                                                                                                         <  1:897145/112‑1 (MQ=255)                    tGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTACtt                                                                                              >  1:289795/1‑139 (MQ=255)                      aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAcc                                                                                                                                     <  2:892116/98‑1 (MQ=255)                      aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTAcc                                                                                                                                     >  1:892116/1‑98 (MQ=255)                      aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCa                                                                                                                               <  2:874231/104‑1 (MQ=255)                      aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCa                                                                                                                               >  1:874231/1‑104 (MQ=255)                      aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                                >  2:1071947/1‑135 (MQ=255)                      aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                                <  1:1071947/135‑1 (MQ=255)                           aCGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACg                                                                                                                          >  1:83842/1‑104 (MQ=255)                           aCGGTCGATTTGATTCAATGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACg                                                                                                                          <  2:83842/104‑1 (MQ=255)                                cGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGGATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCaa                                                                                  <  2:811848/139‑1 (MQ=255)                                                    tCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCa                                                                                                                                   >  2:557385/1‑70 (MQ=255)                                                    tCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCa                                                                                                                                   <  1:557385/70‑1 (MQ=255)                                                      cTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                        >  1:182401/1‑95 (MQ=255)                                                      cTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATg                                                                                                        <  2:182401/95‑1 (MQ=255)                                                                                  aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTgg                                                                                       >  2:1063144/1‑84 (MQ=255)                                                                                  aCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTgg                                                                                       <  1:1063144/84‑1 (MQ=255)                                                                                  aCCTGTGACCAGAATAGAAGTGAGTCAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAcc                                >  1:912423/1‑139 (MQ=255)                                                                                   ccTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCt                               <  1:749057/139‑1 (MQ=255)                                                                                              aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTg                                                                               <  2:955224/80‑1 (MQ=255)                                                                                              aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTg                                                                               >  1:955224/1‑80 (MQ=255)                                                                                              aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTaa                     >  1:628866/1‑138 (MQ=255)                                                                                              aaTAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTaa                     <  2:628866/138‑1 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTgggg                                                                                     <  2:178468/64‑1 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTgggg                                                                                     >  1:178468/1‑64 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGc                                                  >  1:564278/1‑99 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGc                                                  <  2:56136/99‑1 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGc                                                  <  2:564278/99‑1 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGc                                                  >  1:56136/1‑99 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAc                                 <  2:446166/116‑1 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAc                                 >  1:446166/1‑116 (MQ=255)                                                                                                        aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt          >  1:1069733/1‑139 (MQ=255)                                                                                                          ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTg                                                                <  2:545347/83‑1 (MQ=255)                                                                                                          ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTg                                                                >  1:545347/1‑83 (MQ=255)                                                                                                          ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa        >  2:495177/1‑139 (MQ=255)                                                                                                                acatTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGTTcc  >  1:253194/1‑139 (MQ=255)                                                                                                                   |                                                                                                                                        CGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGTTCC  >  NC_000913/2910642‑2910891

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