Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 2,637,439 C→A 9.8% intergenic (‑83/+35) der ← / ← bamB GTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,637,4390CA 9.8% 86.4 / 2.7 39intergenic (‑83/+35)der/bamBGTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (11/24);  new base A (2/2);  total (13/26)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

ACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAA  >  NC_000913_3_hsa_tpiA/2637294‑2637583
                                                                                                                                                 |                                                                                                                                                
gagggtAAATGACGTGGATTTTGCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCGTGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAg                                                                                                                                               <  2:529242/144‑1 (MQ=255)
aCGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAg                                                                                                                                               >  2:472910/1‑149 (MQ=255)
  ggTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGa                                                                                                                                             <  2:352763/149‑1 (MQ=255)
     tAAATAACGTGGGTGTGCCGGCGTGAGGGAGCCAGACAAGTGCGAGCATAGGTACCATGTTTAAATCCTCTTGTTTATAAATCCTCAGACAACGCACGCTATATTGGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGcc                                                                                                                                          <  2:219373/149‑1 (MQ=255)
         tAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGttt                                                                                                                                      <  1:88147/149‑1 (MQ=255)
         caaCGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAACCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGttt                                                                                                                                      <  1:573234/148‑1 (MQ=255)
             gTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACACGAAAACGCCCCCTTTCCAGGAGCCGCTTTCaa                                                                                                                                  >  1:418712/1‑149 (MQ=255)
                        tACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAATCCACCCTTCTATAAATCACCATCCAACGCACGCTCTATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACa                                                                                                                       >  2:628828/1‑149 (MQ=255)
                           gTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGCACAGGAGCCGTTTTCAAAGAGAACGGCagag                                                                                                                    >  1:488080/1‑149 (MQ=255)
                            ttAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCACAGTGAACGAAAaaga                                                                                                                   >  2:536319/1‑149 (MQ=255)
                               gggCGCTCCACAATCGCGACCACAGGCACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAAAAGGAAAACGGCCCCTGTCCAGGAGACGTTTTCAAAGTGAACGACAGAGACGa                                                                                                                <  1:113770/149‑1 (MQ=255)
                                              gcgACCACAGGTACCCTGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATaga                                                                                                 <  1:291684/149‑1 (MQ=255)
                                                         tACCATGGTTAAAGCCTCATTTTTATAAATCAGCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTc                                                                                      <  2:127688/149‑1 (MQ=255)
                                                         caCCAGGTTTAAAGCGTCATGTTTATAAATCATCAGACAACGCACGTGATATTCGCGTCGTTGGTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTGCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTc                                                                                      <  2:127878/148‑1 (MQ=255)
                                                                     aGCCTCACTTGTAAAAAACATGAGACTACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCt                                                                          <  1:641026/149‑1 (MQ=255)
                                                                      gCCTCATTTTTATAAATCATCAGACAACGCCCGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTg                                                                         <  1:420052/149‑1 (MQ=255)
                                                                        cTCATTTTTATAACTGATCAGACATCGGCCGTTTGATTCGCGTGGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCTTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGa                                                                       <  2:71517/149‑1 (MQ=255)
                                                                        cTCATTTCTATAAATCATCAGACAACGCACGCAACATTCGCGTCGCTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGa                                                                       <  1:381737/149‑1 (MQ=255)
                                                                         tCATTTTTATAAGTCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGAt                                                                      <  2:172368/149‑1 (MQ=255)
                                                                            tttttATAAACCATTCGCCAACCGGCGCTATATGCGCGGTGTTGTTAAAACCAGGAACACGGCCCGTGTCCAGGGTCCGTTTACAAAGTGACCGACGGAGACGATTAACGTGTAATAGAGTACACGGTTCAGTCTTTTGCCTGGATcag                                                                   <  2:475688/149‑1 (MQ=255)
                                                                                     atcatcAGACAACTCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGCCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATATAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                          <  2:271231/149‑1 (MQ=255)
                                                                                          caGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAg                                                     >  1:231817/1‑149 (MQ=255)
                                                                                                    acgcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCgg                                           <  2:656462/149‑1 (MQ=255)
                                                                                                    acgcTATATTCGCGTCGTTGGTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTAAGCGAGAGAGACGATTAAGGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGCTCAGACGTTTGCCGTCATCGGCAACCgg                                           >  1:283736/1‑149 (MQ=255)
                                                                                                      gcTATATGCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCTGtt                                         <  2:555204/149‑1 (MQ=255)
                                                                                                             tCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTg                                  <  1:472380/149‑1 (MQ=255)
                                                                                                             cgctgtcgtTGTGAAAGAAAGGAAACGGGCCTCGGTCCAGGAGCGGTTTTCAAATTGACCGACAGGTACGATTAACGTGGAATTGAGGTCATGTTTCCGTCTTTTGCCTGGATCAGCATTTTGCCGTCAGCGGCAACCGGTTCAGTCTg                                  <  2:384712/145‑1 (MQ=255)
                                                                                                                  tcgtTGTTAAAAACACGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAAc                             <  1:229693/149‑1 (MQ=255)
                                                                                                                     ttgttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGCTTCAGTCTGCAAACCgg                          >  2:661096/1‑149 (MQ=255)
                                                                                                                     ttgttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTGTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCgg                          <  1:651974/149‑1 (MQ=255)
                                                                                                                       gttAAAAACCGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGCGTACACGGTTCCGTCTTTTGCCCGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGaa                        <  1:122987/149‑1 (MQ=255)
                                                                                                                       gttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGTCGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGaa                        <  1:689808/149‑1 (MQ=255)
                                                                                                                        ttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGCTTAACGTGTAATAGCGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGCAAc                       >  1:638235/1‑149 (MQ=255)
                                                                                                                        ttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACAGGTTCAGTCTGAAAGCCGGCAc                       >  1:638594/1‑149 (MQ=255)
                                                                                                                               tttgAAAACGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGACACCGGGACTCTCaat                >  2:11985/4‑148 (MQ=255)
                                                                                                                                  gAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACttt             >  1:251927/1‑149 (MQ=255)
                                                                                                                                        gggcccTGGCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACGCGGTTCCGTCTTCTGCCGGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTg       <  1:647545/146‑1 (MQ=255)
                                                                                                                                         gCCCCTGGCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGACTACTGTATATTGTTGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTgg      <  1:568290/149‑1 (MQ=255)
                                                                                                                                             cTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCATCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCaa  >  1:122838/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAA  >  NC_000913_3_hsa_tpiA/2637294‑2637583

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