Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,795,396 C→A Q382K (CAG→AAG)  gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,795,3960CA96.2% 69.6 / NA 26Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (12/13);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTAG  >  NC_000913_3_hsa_tpiA/2795258‑2795540
                                                                                                                                          |                                                                                                                                                
gTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTAGCGCCCTGCTGGTTTATTTAGTCATGGCCGTTTCAATGTTGCGGAt                                                                                                                                        <  2:365259/149‑1 (MQ=255)
     ggTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAGTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAGGGTGCGGATGCGCa                                                                                                                                   >  2:244297/1‑149 (MQ=255)
            ctctCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCt                                                                                                                            <  2:364454/149‑1 (MQ=255)
                     ggCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga                                                                                                                   >  2:528050/1‑149 (MQ=255)
                         cgGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaagga                                                                                                               >  1:341807/1‑149 (MQ=255)
                          ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTAAAAGTTGCGGATGCGCATGATTCTGCGAGCAgaaggaa                                                                                                              >  1:54929/1‑149 (MQ=255)
                             aTTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCg                                                                                                           >  2:236587/1‑149 (MQ=255)
                                                  cTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                      <  1:161709/149‑1 (MQ=255)
                                                   tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                   >  1:223203/1‑135 (MQ=255)
                                                   tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                   <  2:223203/135‑1 (MQ=255)
                                                     ttACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGAGCTATCGCCCTGCTGGTTTATTTAGTCATCGCAGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAACGTGGCttt                                                                                   >  1:269656/1‑149 (MQ=255)
                                                      tACGCGCCTGCGAAGGTATTTAAATCTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTCTTTAGTCATCGCCGTTATAAAGTTGCGGATGCGCAAAGTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTa                                                                                  <  1:528050/149‑1 (MQ=255)
                                                             cTGCGAAGGTATTTAGATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTACAAAGTTGCGGCTGCGAAAAATTCTGCGTACAGAAGAAAGCGACATTAGCTTGCGCATGTGGCTTTCTCCgggg                                                                           >  2:351884/1‑149 (MQ=255)
                                                                    ggTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                                    <  1:173944/149‑1 (MQ=255)
                                                                                  tGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTCGCGGATGCGCAAGATTCTGCGAGCAGAAGGACGCGAAATTCGCTTGCGAATGTGACTTTATCCGTGGCTAACCTGGCTAGTCATCGGc                                                      >  2:314578/1‑149 (MQ=255)
                                                                                           gCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc                                             >  2:359501/1‑149 (MQ=255)
                                                                                                   gcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                      >  2:146479/1‑100 (MQ=255)
                                                                                                   gcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                      <  1:146479/100‑1 (MQ=255)
                                                                                                    ctctCGCCCTGCTGGTTTATTTAGTCCGCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGCATACCAAGGCGGGTCATCGGCTTTATTACCTTTGTGTTg                                    <  1:290661/146‑1 (MQ=255)
                                                                                                    cTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                                    <  1:389366/149‑1 (MQ=255)
                                                                                                    cTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                                    >  1:36333/1‑149 (MQ=255)
                                                                                                                    gtatttaGTCATCTCCACTTCAAAGTGGCGGATGCGCAAAATTCTCCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtcc                    <  2:456110/148‑1 (MQ=255)
                                                                                                                                cggcGTTTCAAAGGTGCGGGGGCAAAAAATTCTGGGAGCAGAAGGAACCGAAATTCGCTCGCGAATGTGGCTTTATCCGAGGCTCACCTGGCTAGTCATCGTCTTTATTACCTTGGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAAc        <  2:230134/146‑1 (MQ=255)
                                                                                                                                    gTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGtt    <  2:205237/149‑1 (MQ=255)
                                                                                                                                      ttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTAg  >  2:573297/1‑149 (MQ=255)
                                                                                                                                      ttCAAAGTTGCGGATGCGCAAAATTCTGCAAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTAg  <  1:198327/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
GTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTAG  >  NC_000913_3_hsa_tpiA/2795258‑2795540

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: