Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 1,660,221 T→G I718S (ATC→AGC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,660,2210TG97.5% 114.2 / NA 40I718S (ATC→AGC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (23/16);  minor base A (1/0);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGC  >  NC_000913_3_hsa_tpiA/1660077‑1660364
                                                                                                                                                |                                                                                                                                               
aTGCCTCAACTTTTGAAGGCTGGAATTCTCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAAc                                                                                                                                             <  1:88349/149‑1 (MQ=255)
    cTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACGCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGGGGTGTGGAGCAACCCta                                                                                                                                         >  1:147635/1‑149 (MQ=255)
           tttGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGCGCAACCCTATAGATgc                                                                                                                                  >  2:16363/1‑149 (MQ=255)
             tggaGGGTGGAACGCCCCTGAAAGTAAAACCGTCCCACTGCAACTGTTTGGTTCCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATgcgc                                                                                                                                <  2:610950/146‑1 (MQ=255)
                  gCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGCCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGCTGCGCAGaaa                                                                                                                           >  1:182019/1‑149 (MQ=255)
                              cTGAACGTAGAACCCTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGAAGAAACGTGGGATTggc                                                                                                               >  2:574207/1‑147 (MQ=255)
                                   cGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACgg                                                                                                          <  1:359779/149‑1 (MQ=255)
                                    gTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATCATGATATCCTGAAGGCTGCCTGCTGTCGGAAGGGGTGAAGCAACCCTATAGAGGCGCGGAAACGCGGGTTTGCCAACgcg                                                                                                         >  2:211353/1‑147 (MQ=255)
                                    gTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAACCGAGGGATTGCCCACGGc                                                                                                         >  2:98071/1‑149 (MQ=255)
                                                   tGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCAGTATGATGCGGgtg                                                                                          >  1:301593/1‑149 (MQ=255)
                                                     gaaCTGTTCGGTTGCCATTAAAAATCCCGTTCGCACCCGACCTACGGCAATATGGAGCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                                                        <  2:607973/148‑1 (MQ=255)
                                                      aaCTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGttt                                                                                       <  2:693727/149‑1 (MQ=255)
                                                           tttGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTTTGCAGAAACCCTATAGATGCGAAAAAACGTGGGATTGCCAACGCCAATATGGTGCGAGTGTTTAACCa                                                                                  >  1:131262/1‑149 (MQ=255)
                                                                tttCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTg                                                                             >  2:158372/1‑149 (MQ=255)
                                                                tttCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATGGATATATTCAAGGCTGCATACCGTCGGGGCGCGTGGAGCAAGCCTATAAAAGAGCAGAAACTTGGGATTGCCGAAGGGGATATGGTGAGGCTGTTTAACGATCGTg                                                                             >  2:121841/1‑149 (MQ=255)
                                                                 ttacATTACAAACCCCGTACTAACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTgg                                                                            <  1:425132/146‑1 (MQ=255)
                                                                    cATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGa                                                                         >  1:77664/1‑149 (MQ=255)
                                                                      ttACAAATCCCGTACTCACTCGTCCTACGGCAATATTGATCTACTGAGCGCTGCCTGCCGTCGGGAGGCGTGGAGCCACCCTATACATGCACATAGACGTGGGATTGCCAACGGCGAGAGTGTGCGGGCGTTTATCCATCTTGGCGagc                                                                       >  2:271950/1‑147 (MQ=255)
                                                                      ttACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAg                                                                       >  1:715645/1‑149 (MQ=255)
                                                                      ttACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAg                                                                       >  2:271541/1‑149 (MQ=255)
                                                                      ttACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGATGCCTACAGTTAGGGTGATTTGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAg                                                                       >  2:344003/1‑149 (MQ=255)
                                                                          aaaTCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGGGGTGTGGAACAACCCTATAGATGCGCAGAAACGTGGGACTCCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTAGCGAAGTTCg                                                                   >  1:604812/1‑149 (MQ=255)
                                                                          aaaTCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCg                                                                   >  1:274217/1‑149 (MQ=255)
                                                                           aagcccGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCgg                                                                  <  2:171223/146‑1 (MQ=255)
                                                                               ccGTACTCACTCGACCTCTGGCAATATTGATCTCCTTAAGGCGGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTAc                                                              <  1:191420/149‑1 (MQ=255)
                                                                                      ccgtcGACCGACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCaa                                                       <  2:182019/146‑1 (MQ=255)
                                                                                         tcGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGt                                                    >  1:588283/1‑149 (MQ=255)
                                                                                           gACCTACGGCAATATTGATCTCCTGAACGCCGGGGATCGGCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAACGTaa                                                  >  2:703802/1‑149 (MQ=255)
                                                                                              cTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGAGATGGTGCGGGTGTTTAACCCTCGTGCCGAAGTTCGGCTACCCCCCAAAGTCaccc                                               >  2:32587/1‑147 (MQ=255)
                                                                                                 cGGAAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAAccccac                                            <  2:645158/149‑1 (MQ=255)
                                                                                                             tctcCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGCTGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTCTTCTCCCTg                                >  1:195487/1‑149 (MQ=255)
                                                                                                             tctcCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTAc                                                              <  2:20599/119‑1 (MQ=255)
                                                                                                             tctcCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTAc                                                              >  1:20599/1‑119 (MQ=255)
                                                                                                                           cctgccGTCAGGAGGTGTGGAGCAACCCTATAGATGCGGAGAAACGTGGGATTGCCAACGGCGAGATGGGGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATg                  <  2:100055/149‑1 (MQ=255)
                                                                                                                             tgccGTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCTAAGTTCGGCTACCAGCCAAAGTAACACCACGTATtctc                                    <  2:365703/129‑1 (MQ=255)
                                                                                                                             tgccGTCAGGAGGTGTGGAGCAACACTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATtctc                                    >  1:365703/1‑129 (MQ=255)
                                                                                                                                 gTCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAg            <  1:158372/149‑1 (MQ=255)
                                                                                                                                  tCAGGAGGTGTGGAGCAACCCTATAGTTGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACTACGTATTCTCCCTGGAGTTAGCGCTACgggagcaa           <  2:637873/149‑6 (MQ=255)
                                                                                                                                    aggaggTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGa         <  1:515671/149‑1 (MQ=255)
                                                                                                                                           gtgGAGCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGc  >  1:446488/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGC  >  NC_000913_3_hsa_tpiA/1660077‑1660364

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