Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,435 Δ1 bp 14.3% intergenic (+17/+18) ydfZ → / ← ydfI selenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,629,4350G.14.3% 84.9 / 10.4 28intergenic (+17/+18)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base . (2/2);  total (13/15)
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.81e-01

TGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCA  >  NC_000913/1629307‑1629558
                                                                                                                                |                                                                                                                            
tGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                         <  4:218219/149‑1 (MQ=255)
           gTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATAGCGGCATTGCTTTTTTACTTCACGGAat                                                                                              <  3:1256683/149‑1 (MQ=255)
                aCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAt                                                                                                                   <  1:615002/123‑1 (MQ=255)
                aCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAt                                                                                                                   >  2:615002/1‑123 (MQ=255)
                                     ggAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCa                                                                    >  4:692402/1‑149 (MQ=255)
                                                 gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaa                                                        <  3:692402/149‑1 (MQ=255)
                                                            ttGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCta                                                             >  4:873484/1‑133 (MQ=255)
                                                            ttGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCta                                                             <  3:873484/133‑1 (MQ=255)
                                                            ttGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGt                                             >  3:1114448/1‑149 (MQ=255)
                                                                gAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACgg                                                                                  >  1:1244000/1‑108 (MQ=255)
                                                                gAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACgg                                                                                  <  2:1244000/108‑1 (MQ=255)
                                                                   gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGaa                                                                                               >  1:1065874/1‑92 (MQ=255)
                                                                   gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGaa                                                                                               <  2:1065874/92‑1 (MQ=255)
                                                                           ggCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGc                              <  4:1114448/149‑1 (MQ=255)
                                                                                   tGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACaa                                                   >  4:349222/1‑120 (MQ=255)
                                                                                   tGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACaa                                                   <  3:349222/120‑1 (MQ=255)
                                                                                               gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGa                                                      >  2:532150/1‑105 (MQ=255)
                                                                                               gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGa                                                      <  1:532150/105‑1 (MQ=255)
                                                                                                      cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                      <  4:634111/114‑1 (MQ=255)
                                                                                                      cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                      >  3:634111/1‑114 (MQ=255)
                                                                                                           aCTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgc                                                                        >  3:864245/1‑75 (MQ=255)
                                                                                                           aCTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgc                                                                        <  4:864245/75‑1 (MQ=255)
                                                                                                                  gtgtGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                      <  1:603433/102‑1 (MQ=255)
                                                                                                                  gtgtGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                      >  2:603433/1‑102 (MQ=255)
                                                                                                                     tgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaa                                                            <  4:755030/77‑1 (MQ=255)
                                                                                                                     tgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaa                                                            >  3:755030/1‑77 (MQ=255)
                                                                                                                       aaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGa                                                      <  1:1295283/81‑1 (MQ=255)
                                                                                                                       aaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGa                                                      >  2:1295283/1‑81 (MQ=255)
                                                                                                                            cgccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca  >  1:557259/1‑128 (MQ=255)
                                                                                                                            cgccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca  <  2:557259/128‑1 (MQ=255)
                                                                                                                                |                                                                                                                            
TGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCA  >  NC_000913/1629307‑1629558

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: