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_0009132,910,7560CT100.0% 105.1 / NA 33*112* (TAG→TAAmazFmRNA interferase toxin, antitoxin is MazE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/15);  total (18/15)

GCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCA  >  NC_000913/2910626‑2910885
                                                                                                                                  |                                                                                                                                 
gCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAAt                                                                                                                           <  1:470989/139‑1 (MQ=255)
gCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAAt                                                                                                                           <  2:248946/139‑1 (MQ=255)
                                     aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                      >  1:250971/1‑139 (MQ=255)
                                     aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                      >  2:315039/1‑139 (MQ=255)
                                     aGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTcc                                                                                      >  1:431275/1‑139 (MQ=255)
                                             gTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGc                                                                              >  1:298427/1‑139 (MQ=255)
                                             gTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGc                                                                              >  2:26707/1‑139 (MQ=255)
                                                                             tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAg                                                                                                                        >  1:295416/1‑65 (MQ=255)
                                                                             tGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAg                                                                                                                        <  2:295416/65‑1 (MQ=255)
                                                                                           gTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACtt                                <  2:250971/139‑1 (MQ=255)
                                                                                               aaaCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGa                                                                                                 <  2:649237/70‑1 (MQ=255)
                                                                                               aaaCCTGTGACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGa                                                                                                 >  1:649237/1‑70 (MQ=255)
                                                                                                     gtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                          >  1:169024/1‑71 (MQ=255)
                                                                                                     gtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                          <  2:169024/71‑1 (MQ=255)
                                                                                                     gtgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGAt                      <  2:278602/139‑1 (MQ=255)
                                                                                                     gtgACCAGAATAGAAGTGAATTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGAt                      <  1:94041/139‑1 (MQ=255)
                                                                                                      tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAc                           >  1:415093/1‑133 (MQ=255)
                                                                                                      tgACCAGAATAGAAGTGAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTAc                           <  2:415093/133‑1 (MQ=255)
                                                                                                                      gAGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCa     <  1:26707/139‑1 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAAt                                                                                                   <  2:141951/44‑1 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAAt                                                                                                   >  1:141951/1‑44 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                          <  1:260517/53‑1 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTaa                                                                                          >  2:260517/1‑53 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGc                                                                              >  1:428509/1‑65 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGc                                                                              <  2:428509/65‑1 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGc                  <  1:449840/125‑1 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGc                  >  2:449840/1‑125 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  2:795212/1‑139 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  2:528006/1‑139 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  2:335686/1‑139 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  1:401832/1‑139 (MQ=255)
                                                                                                                       aGTTAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCAt    >  1:774450/1‑139 (MQ=255)
                                                                                                                         ttAGTAACATTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCa  <  1:528006/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                 
GCCGTTTTCCGGATCGCGCAGGCGCTGCATAATAGTGAGCAAACGGTCGATTTGATTCATTGAATTGTCCTGAAAATTGCGGGTCTGTCAGGTGGAAACCTGTGACCAGAATAGAAGTGAGTTAGTAACACTACCCAATCAGTACGTTAATTTTGGCTTTAATGAGTTGTAATTCCTCTGGGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCA  >  NC_000913/2910626‑2910885

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: