Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 1,734,464 G→C 20.3% T10T (ACG→ACC mepH → murein DD‑endopeptidase, space‑maker hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,734,4640GC20.3% 42.5 / 5.2 24T10T (ACG→ACCmepHmurein DD‑endopeptidase, space‑maker hydrolase
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base C (3/2);  total (13/11)
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.16e-01

ATAAAAGCTATTAACTTTCTCTTCTTCTATGCATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGCAGATAAAAAGAAAAGCACGGCAACCACCAAAAAGACCCAGAA  >  NC_000913_3_hsa_tpiA/1734335‑1734607
                                                                                                                                 |                                                                                                                                               
aTAAAATCTATTAACTTTCTCTTCTTTGATGCGTTCGAATGCTCAAGTTTTGTAAATCAGACGCAGGCACGATAGACCTGTCTTTACAGAGGGACGCTTAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTa                                                                                                                              <  2:315931/149‑1 (MQ=255)
    aaGCTATTAACTTTCTCTTCTTCTATGAATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTaccac                                                                                                                          <  2:361343/149‑1 (MQ=255)
    aaGCTATCAACTTACTCTTCTTCTATGTATGAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTaccac                                                                                                                          <  2:361177/149‑1 (MQ=255)
        tATTAACTTTCTCTTCTTCTATGCATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTg                                                                                                                      <  1:3688/149‑1 (MQ=255)
        tATTAACTTTCTCTTCTTCTATGCATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGAACCCTCAGTAGCGCGGATAAACCGAATTTCGATCACCCTCTGTGCTTTACTTTTTACCCCCCTg                                                                                                                      >  2:527175/1‑149 (MQ=255)
                ttttcttcttctATGCATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAAc                                                                                                              <  2:648850/146‑1 (MQ=255)
                   tcttcttctATGCATTAGAATCATCAAGTTTTGTAAACCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGcc                                                                                                           >  1:363235/1‑149 (MQ=255)
                   tattcttctAAGCATTAGAATCATCAAGTTTTGTAAAACAGAAGCAGGCATGATAGACCTGACTTTACAGATGGACGCTCAGTGGCGCGGATAAAACGTATTTCAATCACGCCCTGTACTTTACTTTTTCCCACCCCTCCTTTAAcccc                                                                                                           >  2:230080/3‑149 (MQ=255)
                        ttctATGCATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTCTTTCGATCACGCTCTGTGCTTTACATTTTACCACCCTGCCTTTAACGCCTATgg                                                                                                      >  1:567319/1‑149 (MQ=255)
                         tctATGCATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGc                                                                                                     >  1:396818/1‑149 (MQ=255)
                                      atcatcAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAaag                                                                                        <  1:178517/149‑1 (MQ=255)
                                                         cAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACCCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTg                                                                     <  2:118861/149‑1 (MQ=255)
                                                                      gATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCatat                                                        <  2:478162/149‑1 (MQ=255)
                                                                           aCCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGAGAGAGTTCTGCTACCACTCATATCACCa                                                   >  2:438882/1‑149 (MQ=255)
                                                                           aCCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAAACGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCAATCATATCCCCa                                                   >  1:387335/1‑149 (MQ=255)
                                                                                  tttACAGAGGGACGCTCAGTGGCGCGGATAAACCGTGTTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGc                                            >  1:388338/1‑149 (MQ=255)
                                                                                     acatacggACGCTCAGTGGGGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTAATACCCTGCCTTTATCGTTTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTATCATTCATCTCACCAAAAAAGCAGa                                         <  2:515962/143‑1 (MQ=255)
                                                                                     aCAGAGTGACGCTCAGTGGCGCGGATCATCCCTATTTCGATCACCCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGCAGa                                         <  1:168674/149‑1 (MQ=255)
                                                                                     aCAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGGTCTGCTACCCCTCATATCACCAACAAAGCAGa                                         >  2:424768/1‑149 (MQ=255)
                                                                                                        cgcgGATAAACCGTATTTCGATCACCCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGCAGATAAAATCAAAAacacacca                      >  2:657885/1‑145 (MQ=255)
                                                                                                         gcgGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGCAGATAAAAAGAAAAGCACGGCaa                     >  1:208488/1‑149 (MQ=255)
                                                                                                                       tttCGATCCCGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGCAGATAAAAAGAAAAGCACGGCAACCACCAAAAAGAcc       <  1:370946/149‑1 (MQ=255)
                                                                                                                       tttCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGCAGAGAAAAAGAAAAGCACGGCAACCACCAAGAAGAcc       >  1:107812/1‑149 (MQ=255)
                                                                                                                            aTCACCCTCTGTGCTTTACTTTTTACAACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGCGAGTTCTGCTACCACTCATATCACCAAAAAAACAGATAAAAAGAAAAGCACGGCAACCACCAAAAAGACCCAGaa  >  1:626469/1‑149 (MQ=255)
                                                                                                                                 |                                                                                                                                               
ATAAAAGCTATTAACTTTCTCTTCTTCTATGCATTAGAATCATCAAGTTTTGTAAATCAGACGCAGGCATGATAGACCTGCCTTTACAGAGGGACGCTCAGTGGCGCGGATAAACCGTATTTCGATCACGCTCTGTGCTTTACTTTTTACCACCCTGCCTTTAACGCCTATGGCCCATGCTTCAAAGCAAGCCAGGGAGAGTTCTGCTACCACTCATATCACCAAAAAAGCAGATAAAAAGAAAAGCACGGCAACCACCAAAAAGACCCAGAA  >  NC_000913_3_hsa_tpiA/1734335‑1734607

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