Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,854 A→T 15.4% intergenic (‑141/‑142) yefM ← / → hisL antitoxin of the YoeB‑YefM toxin‑antitoxin system/his operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,089,8540AT15.4% 69.5 / 7.2 26intergenic (‑141/‑142)yefM/hisLantitoxin of the YoeB‑YefM toxin‑antitoxin system/his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (11/11);  new base T (2/2);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

TCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTA  >  NC_000913/2089737‑2089953
                                                                                                                     |                                                                                                   
tCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                             <  2:804635/126‑1 (MQ=255)
tCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                             >  1:804635/1‑126 (MQ=255)
         aaTGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                             <  4:1047664/117‑1 (MQ=255)
         aaTGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                             >  3:1047664/1‑117 (MQ=255)
                        ttaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                    <  3:1101251/95‑1 (MQ=255)
                        ttaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                    >  4:1101251/1‑95 (MQ=255)
                                      attaGTTTGTAGACCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCctct                                                                                              <  4:448081/87‑1 (MQ=255)
                                      attaGTTTGTAGACCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCctct                                                                                              >  3:448081/1‑87 (MQ=255)
                                                      aGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                          <  2:26732/107‑1 (MQ=255)
                                                      aGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                          >  1:26732/1‑107 (MQ=255)
                                                               cAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatata                                       >  2:732548/1‑117 (MQ=255)
                                                               cAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatata                                       <  1:732548/117‑1 (MQ=255)
                                                                                aGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           <  4:673862/128‑1 (MQ=255)
                                                                                aGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           >  3:673862/1‑128 (MQ=255)
                                                                                 ggCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                         <  1:831535/113‑1 (MQ=255)
                                                                                 ggCTTCGCTATTTTGGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                         >  2:831535/1‑113 (MQ=255)
                                                                                        cTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg      >  1:981238/1‑125 (MQ=255)
                                                                                        cTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg      <  2:981238/125‑1 (MQ=255)
                                                                                          aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaa                                            >  1:423124/1‑85 (MQ=255)
                                                                                          aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaa                                            <  2:423124/85‑1 (MQ=255)
                                                                                            ttttGGTGCCATCAGCTAAGAGGACATTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                   >  3:934991/1‑60 (MQ=255)
                                                                                            ttttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                   <  4:934991/60‑1 (MQ=255)
                                                                                                        cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTa  <  4:1111096/113‑1 (MQ=255)
                                                                                                        cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTa  >  3:1111096/1‑113 (MQ=255)
                                                                                                             aaGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                         >  3:925797/1‑85 (MQ=255)
                                                                                                             aaGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                         <  4:925797/85‑1 (MQ=255)
                                                                                                                     |                                                                                                   
TCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTA  >  NC_000913/2089737‑2089953

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