Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,089,854	A→T	10.0%	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_000913	2,089,854	0	A	T	10.0%	121.7 / 6.4	40	intergenic (‑141/‑142)	yefM/hisL	antitoxin of the YoeB‑YefM toxin‑antitoxin system/his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base T (2/2);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-02

GTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGCTGTAGCTAATTGTACGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAA  >  NC_000913/2089590‑2090120                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                           gTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGCTGTAGCTAATTGTACGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTCTTAATTAGTTTGTAGATCAAGGTATTGTCCGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGAAAGTCCTCTTa                                                                                                                                                                                                                                                                   >  2:286264/1‑274 (MQ=255)                                     taataaGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGCTGTAGCTAATTGTACGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATcc                                                                                                                                                                                                                              >  2:486057/1‑275 (MQ=255)                                                                                                                                                    tCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                                                                                                                                                                                                    <  6:804635/126‑1 (MQ=255)                                                                                                                                                    tCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                                                                                                                                                                                                    >  5:804635/1‑126 (MQ=255)                                                                                                                                                             aaTGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                                                                                                                                                                                                    <  8:1047664/117‑1 (MQ=255)                                                                                                                                                             aaTGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCtt                                                                                                                                                                                                                                                                    >  7:1047664/1‑117 (MQ=255)                                                                                                                                                               tgatgaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGCCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCATCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAACGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACCTCAGTTGAATAAACATTCACAGAGACTTTTATGTCACGCGTTCAATTTAAAcaccac                                                                                                    >  2:368985/1‑275 (MQ=255)                                                                                                                                                               tgatgaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAACAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAacacacc                                                                                                    >  4:525892/1‑272 (MQ=255)                                                                                                                                                                            ttaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                                                                                                                                                                                           >  8:1101251/1‑95 (MQ=255)                                                                                                                                                                            ttaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                                                                                                                                                                                           <  7:1101251/95‑1 (MQ=255)                                                                                                                                                                                cTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATcc                                                                                   >  4:138141/1‑275 (MQ=255)                                                                                                                                                                                          attaGTTTGTAGACCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCctct                                                                                                                                                                                                                                                                     <  8:448081/87‑1 (MQ=255)                                                                                                                                                                                          attaGTTTGTAGACCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCctct                                                                                                                                                                                                                                                                     >  7:448081/1‑87 (MQ=255)                                                                                                                                                                                                          aGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                                                                                                                                                                                                 >  5:26732/1‑107 (MQ=255)                                                                                                                                                                                                          aGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                                                                                                                                                                                                 <  6:26732/107‑1 (MQ=255)                                                                                                                                                                                                                tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGaa                                                    >  2:21064/1‑274 (MQ=255)                                                                                                                                                                                                                   cAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatata                                                                                                                                                                                                              >  6:732548/1‑117 (MQ=255)                                                                                                                                                                                                                   cAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatata                                                                                                                                                                                                              <  5:732548/117‑1 (MQ=255)                                                                                                                                                                                                                      tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGGGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTc                                             >  4:23439/1‑275 (MQ=255)                                                                                                                                                                                                                          aCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAg                                                                                                                                                                                                                                                                           <  2:68893/49‑1 (MQ=255)                                                                                                                                                                                                                                    aGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa                                                                                                                                                                                  <  8:673862/128‑1 (MQ=255)                                                                                                                                                                                                                                    aGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa                                                                                                                                                                                  >  7:673862/1‑128 (MQ=255)                                                                                                                                                                                                                                     ggCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                                                                                                                                                                                                <  5:831535/113‑1 (MQ=255)                                                                                                                                                                                                                                     ggCTTCGCTATTTTGGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                                                                                                                                                                                                >  6:831535/1‑113 (MQ=255)                                                                                                                                                                                                                                          cGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAANAACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTAAATTTAAACACCACCATCATCACCATCNACATGACTAGTTTTTCAGGCGAGGTTGTCTNNAAGAAATTAAAAATCTCAAGTGGTGCATg                         >  3:70301/1‑275 (MQ=255)                                                                                                                                                                                                                                          cGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACAcg                                                                                                                       >  1:80731/1‑181 (MQ=255)                                                                                                                                                                                                                                            cTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg                                                                                                                                                                             >  5:981238/1‑125 (MQ=255)                                                                                                                                                                                                                                            cTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg                                                                                                                                                                             <  6:981238/125‑1 (MQ=255)                                                                                                                                                                                                                                              aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaa                                                                                                                                                                                                                   <  6:423124/85‑1 (MQ=255)                                                                                                                                                                                                                                              aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaa                                                                                                                                                                                                                   >  5:423124/1‑85 (MQ=255)                                                                                                                                                                                                                                              aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGtt                                                                                                                   <  3:80731/181‑1 (MQ=255)                                                                                                                                                                                                                                                ttttGGTGCCATCAGCTAAGAGGACATTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                                                                                                                                                                                          >  7:934991/1‑60 (MQ=255)                                                                                                                                                                                                                                                ttttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                                                                                                                                                                                          <  8:934991/60‑1 (MQ=255)                                                                                                                                                                                                                                                    ggTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATa                                                                                                                                                                                                                     >  1:25572/1‑77 (MQ=255)                                                                                                                                                                                                                                                        ccATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAat                                                                                                                                                                                                                 <  3:25572/77‑1 (MQ=255)                                                                                                                                                                                                                                                            cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTa                                                                                                                                                                         >  7:1111096/1‑113 (MQ=255)                                                                                                                                                                                                                                                            cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTa                                                                                                                                                                         <  8:1111096/113‑1 (MQ=255)                                                                                                                                                                                                                                                                 aaGAGTACAGTCCTCTTAGCCACCTCCTTTCCCCGCTCAGTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTGGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGGTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGaa  <  4:486056/275‑1 (MQ=255)                                                                                                                                                                                                                                                                 aaGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                                                                                                                                                                                                >  7:925797/1‑85 (MQ=255)                                                                                                                                                                                                                                                                 aaGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                                                                                                                                                                                                <  8:925797/85‑1 (MQ=255)                                                                                                                                                                                                                                                                         aGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAg                                           >  4:372336/1‑226 (MQ=255)                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                           GTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGCTGTAGCTAATTGTACGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAA  >  NC_000913/2089590‑2090120

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