Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,679,525 G→T A1074D (GCT→GAT)  yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,679,5250GT97.4% 106.4 / ‑4.2 39A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (20/17);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGA  >  NC_000913_3_hsa_tpiA/2679384‑2679655
                                                                                                                                             |                                                                                                                                  
tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGAtt                                                                                                                             >  1:708646/1‑149 (MQ=255)
     attaATAATCGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAATACAACAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCg                                                                                                                        <  1:163688/149‑1 (MQ=255)
        aacaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGCGAATAATGCCGCTTTGCGCCAATCCGGATTAATCGcca                                                                                                                     >  2:336236/1‑148 (MQ=255)
              ggCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAtt                                                                                                               >  2:681272/1‑149 (MQ=255)
                 ttATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAttaca                                                                                                            >  2:417748/1‑146 (MQ=255)
                    ttCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgcgc                                                                                                         <  1:518145/149‑1 (MQ=255)
                       tAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAg                                                                                                      >  1:220548/1‑149 (MQ=255)
                                 gtgATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCata                                                                                            >  2:427371/1‑149 (MQ=255)
                                     tgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc                                                                                        <  1:427371/149‑1 (MQ=255)
                                      gaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAGAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAttacc                                                                                       >  1:28162/1‑146 (MQ=255)
                                        gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa                                                                                     <  1:344628/149‑1 (MQ=255)
                                              gcgcGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcc                                                                               <  2:220548/149‑1 (MQ=255)
                                                gcgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                             <  1:102577/149‑1 (MQ=255)
                                                         tCAAATTTTGCGGTGTATAGGGTTTGGTGGACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGTCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGGCTCATATTGCGTAACCTCCCCGAGTCCcgc                                                                    <  2:820869/149‑1 (MQ=255)
                                                           aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCTCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCGCGGGTCCCGCgg                                                                  >  2:269588/1‑149 (MQ=255)
                                                             atatTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGACCCACAGAc                                                                >  2:73373/1‑149 (MQ=255)
                                                                 tACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCa                                                            <  1:417748/149‑1 (MQ=255)
                                                                    aGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCa                                                         <  2:221256/149‑1 (MQ=255)
                                                                            gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCCGCATCGccacc                                                 >  1:11648/1‑149 (MQ=255)
                                                                                      aaCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                                                                                                <  1:423911/60‑1 (MQ=255)
                                                                                      aaCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                                                                                                >  2:423911/1‑60 (MQ=255)
                                                                                         gctcttAAAAATAAAAGGCATCACGGTGGTGAATCATGCAGCTTTCGACCAATACGGATTAATCTCCGTCAATTGCGCGCGAGCAGACTCATATTGAGCAACCTCCCCAAATAACGCGGACACAATCATTCCCCCCATCAGGCAttttt                                    >  2:358018/5‑149 (MQ=255)
                                                                                         gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCCTTCAGATGTGCGCGATTCATCTCCGGCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAttttt                                    >  1:197644/1‑149 (MQ=255)
                                                                                                 aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgc                            >  2:608868/1‑149 (MQ=255)
                                                                                                 aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgc                            >  1:830100/1‑149 (MQ=255)
                                                                                                  aaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCCTCGCCCCCATCAGCCATTTTTCTTTAtgct                           >  1:212430/1‑149 (MQ=255)
                                                                                                        aGGCTTCACGTTGGTGAATAATGCCGTTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCATACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGTATTTTTCTTTATGCTGCtgtt                     <  2:205140/149‑1 (MQ=255)
                                                                                                        aGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgtt                     <  1:681272/149‑1 (MQ=255)
                                                                                                        aGGCATCACGGTGGTGAATAATGCCGCCTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCTTATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgtt                     <  2:256712/149‑1 (MQ=255)
                                                                                                            aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCGAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTaa                 >  1:155778/1‑149 (MQ=255)
                                                                                                             tCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATc                                                       <  1:502307/110‑1 (MQ=255)
                                                                                                             tCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATc                                                       >  2:502307/1‑110 (MQ=255)
                                                                                                              cACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCACCTGCTTTTCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAAt               >  1:193822/1‑149 (MQ=255)
                                                                                                               aCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATc              <  1:30676/149‑1 (MQ=255)
                                                                                                                 tgggttgAATGATGCCGCTTTCGGCTAACCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCCACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGc            <  2:479477/144‑1 (MQ=255)
                                                                                                                   tggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCGGTCAATTCAGCGCGGGCAGAATCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCGTCTGTTGTGAATCGcga          >  2:117536/1‑147 (MQ=255)
                                                                                                                   tggtgAATAATGCCCCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCa          <  1:414262/149‑1 (MQ=255)
                                                                                                                       gAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCCGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCCTACa      >  1:158001/1‑149 (MQ=255)
                                                                                                                           aatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGa  <  1:73373/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                  
TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGA  >  NC_000913_3_hsa_tpiA/2679384‑2679655

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