Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,911,083 G→A S3S (AGC→AGT 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,911,0830GA82.1% 120.4 / 22.6 56S3S (AGC→AGTmazFmRNA interferase toxin, antitoxin is MazE
Reads supporting (aligned to +/- strand):  ref base G (6/4);  new base A (23/23);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAAT  >  NC_000913/2910946‑2911217
                                                                                                                                         |                                                                                                                                      
acacacAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGc                                                                                                                                       <  2:456137/139‑1 (MQ=255)
  acacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACtt                                                                                                                                     >  1:630510/1‑139 (MQ=255)
  acacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACtt                                                                                                                                     >  1:422468/1‑139 (MQ=255)
             tACCTGTTTTGTTGTTTTACATGAAAGGACTCAGGAAAACAGCTGGACGATGTCAAGCTTGCTCGTTACCTTTGGTCGGGTAAAAATCAACCCAAATAAGATCGCCCATTTCGGGTACGTATCGACttaccattacaga                                                                                                                          <  1:292977/139‑4 (MQ=255)
             tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAg                                                                                                                          >  1:188046/1‑139 (MQ=255)
             tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAg                                                                                                                          <  2:188046/139‑1 (MQ=255)
              aCCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGa                                                                                                                         >  2:642554/1‑139 (MQ=255)
                  gTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTc                                                                                                                     >  2:378385/1‑139 (MQ=255)
                    tttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCt                                                                                                                   >  1:577333/1‑139 (MQ=255)
                           ttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCttt                                                                                                            <  1:161848/139‑1 (MQ=255)
                            tgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCt                                                                                                                   >  2:711617/1‑131 (MQ=255)
                            tgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCt                                                                                                                   <  1:711617/131‑1 (MQ=255)
                                  tGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGctct                                                                                                     >  2:107023/1‑139 (MQ=255)
                                   gAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCtctc                                                                                                    >  1:517234/1‑139 (MQ=255)
                                    aaaGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                                                                                   <  1:305140/139‑1 (MQ=255)
                                           tCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGGCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCg                                                                                            <  2:630510/139‑1 (MQ=255)
                                             aGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                                                                                          >  1:330989/1‑139 (MQ=255)
                                             aGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                                                                                          >  2:650270/1‑139 (MQ=255)
                                                    caGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCg                                                                                   >  1:298249/1‑139 (MQ=255)
                                                         ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTccc                                                                                                  >  1:89773/1‑119 (MQ=255)
                                                         ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAg                                                                              >  2:119588/1‑139 (MQ=255)
                                                         ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACACAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTccc                                                                                                  <  2:89773/119‑1 (MQ=255)
                                                               tGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTccc                                                                                                  <  2:586213/113‑1 (MQ=255)
                                                               tGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTccc                                                                                                  >  1:586213/1‑113 (MQ=255)
                                                                       ttGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGa                                                                <  1:107023/139‑1 (MQ=255)
                                                                           tCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTc                                                            >  1:307756/1‑139 (MQ=255)
                                                                           tCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTc                                                            >  2:14658/1‑139 (MQ=255)
                                                                            cGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCg                                                           <  2:557403/139‑1 (MQ=255)
                                                                              cTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGtt                                                         <  2:422468/139‑1 (MQ=255)
                                                                                aCCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGa                                                       <  1:650270/139‑1 (MQ=255)
                                                                                aCCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGa                                                       <  2:330989/139‑1 (MQ=255)
                                                                                      tGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGc                                                                                                                                       <  1:313232/53‑1 (MQ=255)
                                                                                      tGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGc                                                                                                                                       >  2:313232/1‑53 (MQ=255)
                                                                                          gggTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGc                                             <  2:517234/139‑1 (MQ=255)
                                                                                           ggTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGc                                                                    >  1:423841/1‑115 (MQ=255)
                                                                                           ggTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGc                                                                    <  2:423841/115‑1 (MQ=255)
                                                                                            gTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCaa                                           >  1:439876/1‑139 (MQ=255)
                                                                                            gTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCaa                                           <  1:642554/139‑1 (MQ=255)
                                                                                                          aaaTCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGc                             <  1:378385/139‑1 (MQ=255)
                                                                                                          aaaTCAGATCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGc                             <  2:33817/139‑1 (MQ=255)
                                                                                                                 aTCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCCCTCCCCAGTCGATATTCTCGTGGAGGGTTTCCGGCGTGATGTCGTTGACCAGTTCAGCATGCGTAACTACGGGTTCTTTAc                      >  2:354187/1‑139 (MQ=255)
                                                                                                                 aTCGCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTAc                      >  2:655937/1‑139 (MQ=255)
                                                                                                                    gCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcccc                                                                                                 <  1:337179/61‑1 (MQ=255)
                                                                                                                    gCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcccc                                                                                                 >  2:337179/1‑61 (MQ=255)
                                                                                                                    gCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCa                   <  2:307756/139‑1 (MQ=255)
                                                                                                                    gCCCATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCa                   <  2:298249/139‑1 (MQ=255)
                                                                                                                      ccATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACggg                              <  2:479380/126‑1 (MQ=255)
                                                                                                                      ccATATCGGGTACGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACggg                              >  1:479380/1‑126 (MQ=255)
                                                                                                                                  cGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGt                                                                                              <  2:611379/50‑1 (MQ=255)
                                                                                                                                  cGTATCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGt                                                                                              >  1:611379/1‑50 (MQ=255)
                                                                                                                                     aTCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAAtaattaat  >  1:641996/1‑139 (MQ=255)
                                                                                                                                     aTCGACTTACCATTACCAGACTTCGTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACgg                               >  2:484124/1‑110 (MQ=255)
                                                                                                                                     aTCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACgg                               <  1:484124/110‑1 (MQ=255)
                                                                                                                                     aTCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGc                             <  2:164704/112‑1 (MQ=255)
                                                                                                                                     aTCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGc                             >  1:164704/1‑112 (MQ=255)
                                                                                                                                     aTCGACTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAAtaattaat  >  2:92105/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                      
ACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAAT  >  NC_000913/2910946‑2911217

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