Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,658,592 C→T A175V (GCC→GTC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,5920CT98.0% 161.6 / NA 49A175V (GCC→GTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/20);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATG  >  NC_000913/1658495‑1658680
                                                                                                 |                                                                                        
cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCt                                                                                        >  1:243679/1‑100 (MQ=255)
  gaagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcg                                                                                      >  2:274418/1‑100 (MQ=255)
     gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcg                                                                                   >  1:36198/1‑100 (MQ=255)
       aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGt                                                                                 >  2:193756/1‑100 (MQ=255)
       aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGt                                                                                 <  1:292696/100‑1 (MQ=255)
        aaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTc                                                                                <  1:148532/100‑1 (MQ=255)
             tCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATAGTTGGCCGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACg                                                                           <  2:680563/100‑1 (MQ=255)
              cGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGc                                                                          >  1:722712/1‑100 (MQ=255)
              cGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGc                                                                          >  2:751829/1‑100 (MQ=255)
               gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGcc                                                                         >  1:1145399/1‑100 (MQ=255)
                aaCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCt                                                                        >  2:324589/1‑100 (MQ=255)
                    atatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgatg                                                                    >  1:1435741/1‑100 (MQ=255)
                      atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGaa                                                                  >  2:1171637/1‑100 (MQ=255)
                        ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACt                                                                >  2:435172/1‑100 (MQ=255)
                          cgaTGAGGGTGTAAATAATCCCTTCTCATCGGGGATCAATGAGGGCAATATGACCCGCTCTTCGCCGTCAGTCTCGGCGGTCAAACGCCTGATGAACTg                                                               <  2:1300093/97‑1 (MQ=255)
                             tgagtttGTACATATTCACTGCTCATCGGGGATCGTTGGCGGCAGTATATCCCCCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGc                                                              <  2:1275514/91‑1 (MQ=255)
                             tGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAc                                                           <  1:981114/100‑1 (MQ=255)
                               aGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcgg                                                         >  1:1230439/1‑100 (MQ=255)
                                ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                        <  2:1188035/100‑1 (MQ=255)
                                 gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcg                                                       >  1:584454/1‑100 (MQ=255)
                                   tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGt                                                     >  1:506822/1‑100 (MQ=255)
                                      atatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCa                                                  >  1:464396/1‑100 (MQ=255)
                                      atatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                   <  2:1120100/99‑1 (MQ=255)
                                       tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                 <  2:1230434/100‑1 (MQ=255)
                                          atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCa                                              >  1:784813/1‑100 (MQ=255)
                                          atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCa                                              >  2:343074/1‑100 (MQ=255)
                                             cAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAcc                                           >  2:1542547/1‑100 (MQ=255)
                                             cAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAcc                                           >  1:1568608/1‑100 (MQ=255)
                                              aGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCa                                          <  1:1488533/100‑1 (MQ=255)
                                               gTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAg                                         <  2:1164281/100‑1 (MQ=255)
                                               gTACCCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAg                                         <  2:1314398/100‑1 (MQ=255)
                                                 cctgtaccgggATGGTTGGCGGCAATCAGTCCCGCTATTCGCCAGCAGTCTCGGCGGTCAAACGCCGGATGGACTGCTACGGCGGTTCACTCAACCa                                          <  2:329400/89‑1 (MQ=255)
                                                 aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                       <  1:1174042/100‑1 (MQ=255)
                                                  ctctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTAt                                      >  2:1027534/1‑100 (MQ=255)
                                                   tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                     <  1:980076/100‑1 (MQ=255)
                                                     ttCGGGGATCGTTGGCGGCAATATGCCCCGCTCTTCGCCATGAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                   <  2:10123/100‑1 (MQ=255)
                                                             tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                           <  1:324589/100‑1 (MQ=255)
                                                                  ggcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTg                      >  1:537389/1‑100 (MQ=255)
                                                                    cggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcg                    >  1:454651/1‑100 (MQ=255)
                                                                     ggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcgc                   <  2:506821/100‑1 (MQ=255)
                                                                          tatGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt              >  2:1336001/1‑100 (MQ=255)
                                                                              aCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc           >  2:997411/1‑99 (MQ=255)
                                                                              aCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt          >  1:1434721/1‑100 (MQ=255)
                                                                                ccGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt        >  1:1386531/1‑100 (MQ=255)
                                                                                 cGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg       >  1:671735/1‑100 (MQ=255)
                                                                                 cGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg       <  2:94301/100‑1 (MQ=255)
                                                                                    tctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa    <  2:795966/100‑1 (MQ=255)
                                                                                     ctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt   >  2:721307/1‑100 (MQ=255)
                                                                                      ttCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATg  <  1:385339/100‑1 (MQ=255)
                                                                                                 |                                                                                        
CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATG  >  NC_000913/1658495‑1658680

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