Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,501,912 T→A I117I (ATT→ATA tehB → tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,501,9120TA92.2% 209.2 / 7.4 64I117I (ATT→ATAtehBtellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base A (30/29);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTT  >  NC_000913/1501814‑1502007
                                                                                                  |                                                                                               
gTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTg                                                                                                <  2:715826/100‑1 (MQ=255)
  tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                              <  2:1286290/100‑1 (MQ=255)
    aTCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCaa                                                                                            >  1:776256/1‑100 (MQ=255)
    aTCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCaa                                                                                            <  2:1370336/100‑1 (MQ=255)
     tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAat                                                                                           <  1:212235/100‑1 (MQ=255)
     tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                           <  2:1344119/100‑1 (MQ=255)
     tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                           <  2:144164/100‑1 (MQ=255)
        gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                        <  1:1212080/100‑1 (MQ=255)
        gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                        <  2:375084/100‑1 (MQ=255)
          aTAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCa                                                                                      >  2:1532829/1‑100 (MQ=255)
               cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                 >  1:1718321/1‑100 (MQ=255)
               cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                 >  1:435238/1‑100 (MQ=255)
               cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTACCTTGAGGCGAAAACGATCCCCGGATTGATAGCCAATATGCAACGtt                                                                                 >  1:1312164/1‑100 (MQ=255)
                 cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                               <  2:986896/100‑1 (MQ=255)
                 cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                               >  1:467907/1‑100 (MQ=255)
                  tcaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                               <  1:1183748/98‑1 (MQ=255)
                  acaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCa                                                                              <  1:624528/100‑1 (MQ=255)
                   caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                             >  1:1898460/1‑100 (MQ=255)
                       tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaaa                                                                         <  1:1444344/100‑1 (MQ=255)
                         aTAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAcc                                                                       <  1:1682177/100‑1 (MQ=255)
                         aTAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAcc                                                                       >  2:1090438/1‑100 (MQ=255)
                            gACTGTTCGATTTTATTCTTTCGACTTTGGTTCTGATGTTCCTTGAGGCTATAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtg                                                                     <  1:1908151/99‑1 (MQ=255)
                            gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                    <  2:1394012/100‑1 (MQ=255)
                                gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                <  2:1718332/100‑1 (MQ=255)
                                   cGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTAc                                                             >  1:675464/1‑100 (MQ=255)
                                    gATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACa                                                            <  2:627453/100‑1 (MQ=255)
                                     aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                                           >  1:94247/1‑100 (MQ=255)
                                           ttCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAt                                                     >  2:162831/1‑100 (MQ=255)
                                           ttCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAt                                                     <  1:1013722/100‑1 (MQ=255)
                                            tCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAtt                                                    <  2:610644/100‑1 (MQ=255)
                                            tCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAtt                                                    >  2:572909/1‑100 (MQ=255)
                                                  gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                              >  2:1068835/1‑100 (MQ=255)
                                                  gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                              <  2:1401411/100‑1 (MQ=255)
                                                    cTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                            <  1:640544/100‑1 (MQ=255)
                                                     tgtgGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGgcggcg                                           >  1:770423/1‑100 (MQ=255)
                                                       tggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCCATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGAt                                         >  1:242008/1‑100 (MQ=255)
                                                         gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                       <  2:1253367/100‑1 (MQ=255)
                                                            cTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATa                                    >  1:750731/1‑100 (MQ=255)
                                                               aTGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCg                                 >  2:938872/1‑100 (MQ=255)
                                                                tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                <  1:472684/100‑1 (MQ=255)
                                                                 gTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAcgg                                 >  1:1896302/1‑96 (MQ=255)
                                                                 gTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCt                               >  2:1653516/1‑100 (MQ=255)
                                                                  ttCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                              >  2:746276/1‑100 (MQ=255)
                                                                    ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                            <  1:1163553/100‑1 (MQ=255)
                                                                      ttGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAtt                           >  2:1286813/1‑99 (MQ=255)
                                                                       gcagcctaccaCCGTCCCCTGGGTCATAGCCAACATGCAACGTTGTACTAAACCCGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTAt                         <  2:1750711/90‑1 (MQ=255)
                                                                        gAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATc                        >  2:19589/1‑100 (MQ=255)
                                                                         aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                       <  2:1510217/100‑1 (MQ=255)
                                                                         aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                       <  2:340260/100‑1 (MQ=255)
                                                                           gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                     <  2:435240/100‑1 (MQ=255)
                                                                           gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                     >  2:732382/1‑100 (MQ=255)
                                                                            cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATg                    >  2:1906986/1‑100 (MQ=255)
                                                                             tAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                   >  1:982964/1‑100 (MQ=255)
                                                                              aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                   <  1:1286806/99‑1 (MQ=255)
                                                                               aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                 >  1:1691484/1‑100 (MQ=255)
                                                                               aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                 >  1:346234/1‑100 (MQ=255)
                                                                                  ccATCCCCGGGTTGATAGCAAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGt              >  2:1135619/1‑100 (MQ=255)
                                                                                      ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc          >  1:718572/1‑100 (MQ=255)
                                                                                      ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc          >  2:810882/1‑100 (MQ=255)
                                                                                       cccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCt         >  2:1334578/1‑100 (MQ=255)
                                                                                         cGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTc       <  2:280336/100‑1 (MQ=255)
                                                                                            gTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCg    <  2:718573/100‑1 (MQ=255)
                                                                                             ttGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGt   <  2:999609/100‑1 (MQ=255)
                                                                                               gATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACTACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGtt  <  1:938871/99‑1 (MQ=255)
                                                                                                  |                                                                                               
GTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTT  >  NC_000913/1501814‑1502007

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