Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,629,247 C→T I11I (ATC→ATT)  ydfZ → putative selenoprotein YdfZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,629,2470CT96.0% 67.7 / ‑4.1 25I11I (ATC→ATT) ydfZputative selenoprotein YdfZ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (12/12);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

GACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTT  >  NC_000913/1629018‑1629470
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
gACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGTTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACGCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACGGTAACGCAATTACCGCTGGCAGCCGTGTTATg                                                                                                                                                                                                            <  1:310330/251‑1 (MQ=255)
    ggTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGtt                                                                                                                                                                                                         >  1:248459/1‑250 (MQ=255)
     gTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGtt                                                                                                                                                                                                         >  1:204265/1‑249 (MQ=255)
     gTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAg                                                                                                                                                                                                       <  1:19212/251‑1 (MQ=255)
            gcggagTAATTCATTTCCTCTCATCCCATCCGGGGTGAGCGTCTTTTCACCCTCCTTATGTCTCATGCCTGCATCAAAAAAGATGTGAGCTTGAACAAAAACAAAAAATATTTCACTCGAAAGGAGGATTTATATTGCGCCCGTTACGAGGGCTTCGAATGTAAATCAGAAAGTAGAAAACACCTATGACTACCTACGATAGAAACCGTAACGCAATTACCATTGGCAGCCGTGTTATGGTTAGCGGcccc                                                                                                                                                                                                >  2:298638/4‑248 (MQ=255)
            gAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCAcc                                                                                                                                                                                                >  1:239732/1‑251 (MQ=255)
             aaGATTAACTCATTTCCGCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTGATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATGTCCCTCGACAGGAGTATTTATATTGCGCCCGGTACGTGGGCTTCGACTGTAAATCAGAAAGGATAAAACACCTATGACGAGCTACGATCGTAACCGTAACGCAATTACCACTGGTAGCCGTGTTATGGTTAGCGGCACCg                                                                                                                                                                                               <  2:280344/251‑1 (MQ=255)
                                                             cccGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCGTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGGGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTg                                                                                                                                               <  1:268446/251‑1 (MQ=255)
                                                                     aTGGCTCATGCATGCATCAAACAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGCCAAGATCCTGTCGATTGATACTGCAGGTCTGACCGCTGAGCAAAt                                                                                                                                        >  2:75584/1‑250 (MQ=255)
                                                                        gCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATccgc                                                                                                                                    >  1:62480/1‑251 (MQ=255)
                                                                                              aTGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGttgttg                                                                                                              >  1:265771/1‑251 (MQ=255)
                                                                                               tgtgAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGCGAAAAAACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTAGGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGATCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGttgttg                                                                                                              >  1:355590/1‑250 (MQ=255)
                                                                                                      ttGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGCCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTtgt                                                                                                      >  1:307149/1‑251 (MQ=255)
                                                                                                            aaaaacaaaaaATATTTACCTCGAAAGGAGTATTTAGATTTCCCCTGTTACGTGGGCTGCGCCTGTAAATCAGAAAGGAGAAAACTCCTATGACGGCCTACGATCGTAACCGTAACGCAATTACCACTGTCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTTGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCGGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAagag                                                                                                <  2:175866/251‑1 (MQ=255)
                                                                                                                   aaaaaTATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTgg                                                                                         >  1:294157/1‑251 (MQ=255)
                                                                                                                     aaaTATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGa                                                                                                                                                             <  2:199473/181‑1 (MQ=255)
                                                                                                                     aaaTATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGa                                                                                                                                                             >  1:199473/1‑181 (MQ=255)
                                                                                                                                       gAGTATTTAAATTGCGCCCGTTACGTTGGATTCGACGGAAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGt                                                                     <  2:204265/251‑1 (MQ=255)
                                                                                                                                           gtttAAATTTCGCCCGTTACGTGGGATTCGAATGTAAATCAGAAAGGAGAAAAAACCTATGACTACCTACGATCGTAACCGTAACGCCATTTCACCTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACGTGCAAGATCCTGGCAATTGATACTGAAGGTCTGACCGCTGAGAAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGAGTCGtctc                                                                  <  2:215903/249‑1 (MQ=255)
                                                                                                                                                       cccGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGAGCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAg                                                     <  1:208165/251‑1 (MQ=255)
                                                                                                                                                               gTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACATACGATCGTAACCGAAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGCTCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCAGATTCGTCTCGGCATGAACTACGCGTGTGac                                             >  1:356207/1‑250 (MQ=255)
                                                                                                                                                                       cGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAATGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGccgccg                                     <  1:248441/251‑1 (MQ=255)
                                                                                                                                                                                  cAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGAGAATGCCGCCGATGGCGGCAtt                          >  1:339077/1‑251 (MQ=255)
                                                                                                                                                                                           gaAAACACCTAGGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGTTGAGCAAATCCGCCGCGGAAAAACCTTAgttg                                                                                                                 <  1:277677/155‑1 (MQ=255)
                                                                                                                                                                                                          gacCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAttt  <  1:146437/251‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
GACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTT  >  NC_000913/1629018‑1629470

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