Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,847,724 (C)6→7 coding (285/1044 nt) selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,847,7181.C94.1% 56.8 / ‑2.9 17coding (291/1044 nt)selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (9/7);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAA  >  NC_000913/1847670‑1847761
                                                 |                                            
aaGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTg                                               >  1:1369050/1‑49 (MQ=255)
aaGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTg                                               >  1:1939014/1‑49 (MQ=255)
 aGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGc                                             >  1:922489/1‑49 (MQ=255)
  gTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGcc                                           >  1:291745/1‑49 (MQ=255)
    ttGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGcccc                                         >  1:70578/1‑49 (MQ=255)
       ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGcccccc                                       >  1:590878/1‑48 (MQ=255)
       ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGcccccc                                       >  1:1028095/1‑48 (MQ=255)
          aTCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCACCCCCATc                                    <  1:2758705/49‑1 (MQ=38)
             ggCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATcgc                                  >  1:555904/1‑48 (MQ=38)
                cAGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATcgcg                                 >  1:336015/1‑46 (MQ=38)
                cAGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAg                              <  1:1329658/49‑1 (MQ=38)
                    cGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGatat                          <  1:627596/49‑1 (MQ=38)
                    cGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGatat                          >  1:243254/1‑49 (MQ=38)
                        gATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCt                      >  1:2815144/1‑49 (MQ=38)
                        gATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCt                      >  1:954721/1‑49 (MQ=38)
                            gCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATg                  >  1:2442025/1‑49 (MQ=38)
                              aaTCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGc                >  1:1170958/1‑49 (MQ=39)
                                  gCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTa            <  1:548881/49‑1 (MQ=39)
                                    cATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGt          >  1:827230/1‑49 (MQ=39)
                                       aaTCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgc       >  1:2220251/1‑49 (MQ=38)
                                        aTCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgcc      <  1:1583369/49‑1 (MQ=38)
                                        aTCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgc       <  1:146506/48‑1 (MQ=38)
                                           ggTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCa   <  1:583114/49‑1 (MQ=38)
                                            gTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCaa  <  1:2511002/49‑1 (MQ=38)
                                                 |                                            
AAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAA  >  NC_000913/1847670‑1847761

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