Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 4,296,363 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,296,3621.C95.5% 61.2 / ‑3.6 22intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (9/12);  total (9/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
*NC_000913_3_hsa_tpiA4,296,3622.G95.5% 60.1 / ‑3.6 22intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (9/12);  total (9/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 = 9.83e-01

GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913_3_hsa_tpiA/4296278‑4296497
                                                                                     ||                                                                                                                                          
gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCGGGCCTATCTTAACCGt                                                                          >  1:215470/1‑150 (MQ=255)
              cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                            >  1:435115/1‑150 (MQ=255)
                tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCAGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                          <  1:171033/150‑1 (MQ=255)
                    aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                      <  1:574424/150‑1 (MQ=255)
                       ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                   <  1:490199/150‑1 (MQ=255)
                                 aCAGAGCAATACGCTGACTGAGCGTGATTTTGTAGGGCGGATAAGGCGTTTA‑‑CGACGCATCCGACATCCATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                         <  2:341245/150‑1 (MQ=18)
                                     cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCGGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                     <  2:215470/150‑1 (MQ=255)
                                     cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                     <  1:626292/150‑1 (MQ=255)
                                     cgcAATAGCCTGATCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CACTTCCGGCAGTAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                     <  1:332590/150‑1 (MQ=18)
                                      gcAATAGCCTAATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTGGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                    <  1:268670/150‑1 (MQ=18)
                                          tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                <  1:465592/150‑1 (MQ=255)
                                            gCCTGATTTGCCGTGATTTTGTCTGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                              <  2:589383/150‑1 (MQ=21)
                                                       cGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   >  1:131754/1‑150 (MQ=255)
                                                        gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTCCCATACAGCTCaag                  >  1:215497/1‑148 (MQ=21)
                                                        gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGCTATCAAACCCTTCCCTACAGCTCAgg                  >  2:260503/1‑150 (MQ=21)
                                                        gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAAAAACCGCTACCACACACTTCAATGCAGCTCatg                  >  2:648935/1‑148 (MQ=17)
                                                            ttttGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGCCTTATCAGGCCTCTCTTAACCGTTGGCTAATGTTTCAAACTCCTCACACCCGTGATCAAACCCTTCCATACAGCCCACGTTCa              >  1:26715/1‑150 (MQ=18)
                                                              ttGTAGGTCGGATACGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc            <  1:168347/150‑1 (MQ=255)
                                                                  aGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt        <  1:51169/150‑1 (MQ=255)
                                                                    gTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc      >  2:345588/1‑150 (MQ=255)
                                                                        gATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCCgcc  >  2:124539/1‑149 (MQ=255)
                                                                        gATAAGGCGTCTACGC‑‑CGCATCCGACATCAATGCCTGCAAGAAAC‑CCCGCCGCGGCTAATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg  <  1:226094/150‑1 (MQ=255)
                                                                                     ||                                                                                                                                          
GACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTG‑ATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913_3_hsa_tpiA/4296278‑4296497

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: ‑