Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 1,196,232 C→T 20.2% T370T (ACC→ACT icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,196,2320CT20.2% 37.5 / 2.3 25T370T (ACC→ACTicdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (12/7);  new base T (2/3);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAA  >  NC_000913_3_hsa_tpiA/1196084‑1196373
                                                                                                                                                    |                                                                                                                                             
gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCACTGTTTGAAGCCACACACGGTACTGCGCCGAAATATGCCGGTGAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGCGATGATGCTGCGCCACGTGGGTTGGAcc                                                                                                                                               >  1:142197/1‑149 (MQ=255)
             tGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCAGGCTCTGTTAGTCTCTCCGCTGAGATGATGCTGCGCCATCTGGGTTGTACTGAAGCGGCTGagc                                                                                                                                  >  2:235949/1‑147 (MQ=255)
                                tGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGTTGGCTTAATTGTTAAAGGTATGGa                                                                                                               >  2:357126/1‑149 (MQ=255)
                                     ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgcg                                                                                                          >  2:433711/1‑149 (MQ=255)
                                      cTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGTaatacac                                                                                                         >  1:225050/1‑144 (MQ=255)
                                      cTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgcgc                                                                                                         >  2:145456/1‑149 (MQ=255)
                                         tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGATCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGTGCAAt                                                                                                      >  2:394536/1‑149 (MQ=255)
                                          ttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTGTTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATc                                                                                                     >  1:212080/1‑149 (MQ=25)
                                             aaGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGATAAAGTAAAGCCGGGCTCTATGATGCTCTCCGCGGAGATGATGCGGCGCCACATGGGTGGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                                                                                  <  1:630540/149‑1 (MQ=37)
                                             aaGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATAGGTTGGACCGAAGCGGCTGACTTACTTGTTAAAGGTATGGAAGGCGCAATCAAc                                                                                                  >  1:566998/1‑149 (MQ=255)
                                               gCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                                                                                                <  1:373088/149‑1 (MQ=255)
                                                  acccacGGAACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGAAGTTGCGCCACATGGGCTCGACCGAAGAGGCTGACTGAATTGTTAGAGCTGTGGACGGAGCAATCAACGCGaa                                                                                             >  2:557070/1‑149 (MQ=255)
                                                          tACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                                     <  2:593252/149‑1 (MQ=255)
                                                               cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                                <  1:163134/149‑1 (MQ=255)
                                                                gcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATTATGCTGAGCCACATGAGTTGTACCGAACCGGCTGACTTAATTGTTAAAGTTATGGAAGGCGGAATCAACACGAAAACCGCACCATATg                                                                               >  2:579194/1‑149 (MQ=255)
                                                                   gAAATATGCCGGTCAGAACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGGTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAGTGAGCAATCAACGCGAATACCGTAGCCTATGACt                                                                            >  2:59013/1‑149 (MQ=255)
                                                                   gAAATATGCCGATAAGGAAAAATTAAATTGTGGCTCTATTATTCTATGCGCTGAGATGATGCTGCGCAAAATGGGTTGGAATGAAGTGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACTTATGACt                                                                            <  2:226602/149‑1 (MQ=18)
                                                                         tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTTGCCGCTTAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaac                                                                      <  2:199509/149‑3 (MQ=14)
                                                                                    aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGACCGTCTgatggat                                                           >  1:54991/1‑149 (MQ=255)
                                                                                             agtcTGGCTCTACCATTCTCCACACTGAGATGGTGCTGCGCCACATGGCTCGGACTGAAGCGCCTGACTTAATTGTCAAAGGTACGGAAAGCGCAATGAACGCGAACACCGTAACCTACGACTTCGAGCGTCTGATGGATGGCGCTaaa                                                  <  2:84957/146‑1 (MQ=255)
                                                                                                                       gagaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGATAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGa                        <  1:645928/149‑1 (MQ=255)
                                                                                                                            gatggctcgccgcATGGGGTGGACCGAAGCGGTTGACTTAGTTGTTGAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                   <  2:460445/142‑1 (MQ=18)
                                                                                                                                ctgcGCGACATGGGTTGGACGGAAGCGGCTGACTTAATTGTGAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATTACTTCGAGTGTCTGATGGATAATGTTTATTTGCTGAAATGTTCAGAGTTTGGTGACGCGatcat               <  2:112085/149‑1 (MQ=17)
                                                                                                                                      cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATAGaaaa         >  1:643042/1‑149 (MQ=255)
                                                                                                                                       acaTGGGCTGGACCGAACCGAATGACTTAATTGTTCGAGGCCTGGAAGGCTCAATCAACGCGAAAACCGTGACCTATGACTGCAAGTGTCCGATGGATGGCGCTAAACTGCTGACATGTTCTGAGTTTGGTGACGCCGTCATCGCAAAc        >  2:71567/1‑149 (MQ=255)
                                                                                                                                             gTTCGACCGAAGTGGCTGACTTAACTGTTAAAGGTATGTAAGGCGCTATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGGGCGAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa  <  2:340882/149‑1 (MQ=17)
                                                                                                                                                    |                                                                                                                                             
GTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAA  >  NC_000913_3_hsa_tpiA/1196084‑1196373

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