Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 1,207,234 G→T 7.8% A108S (GCC→TCC)  ymfQ → prohage e14 tail protein homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,207,2340GT 7.8% 171.2 / 5.5 52A108S (GCC→TCC) ymfQprohage e14 tail protein homolog
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base T (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

AAACCGCTGGGAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATGCCAGCCACCACCAACTCCACC  >  NC_000913_3_hsa_tpiA/1207089‑1207374
                                                                                                                                                 |                                                                                                                                            
aaaCCGCTGGGAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTTCCCt                                                                                                                                          >  1:100112/1‑150 (MQ=255)
aaaCCGCTGGGAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTTCCCt                                                                                                                                          >  1:100020/1‑150 (MQ=255)
aaaCCGCTGGGAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCt                                                                                                                                          >  4:383139/1‑150 (MQ=255)
  aCCGCTGGGAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTgg                                                                                                                                        >  2:559379/1‑150 (MQ=255)
  aCCGCTGGGAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTgg                                                                                                                                        >  2:83147/1‑150 (MQ=255)
  aCCGCTGGGAGCGTCTGTGCGGTCTGCCGGAGGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGAATTTTATCTTGCACAGCTTGCTGCCCTgg                                                                                                                                        >  4:226365/1‑150 (MQ=255)
          gAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCa                                                                                                                                >  4:286317/1‑150 (MQ=255)
          gAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTACCTTGCACAGCTTGCTGCCCTGGGCAGACCa                                                                                                                                >  1:375284/1‑150 (MQ=255)
            gCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGa                                                                                                                              <  3:100993/150‑1 (MQ=255)
                  gtgCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGccac                                                                                                                        >  1:566929/1‑150 (MQ=255)
                   tgCGGTCTGCCGGATGAATGTATTCCGGCGTGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGACGGTTAACCTCGCGGGCGGCATCAACGAGGATTTGTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGccacc                                                                                                                       <  4:355147/150‑1 (MQ=255)
                          tGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCAcgc                                                                                                                <  1:320073/150‑1 (MQ=255)
                                  gAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGAc                                                                                                        <  1:563791/150‑1 (MQ=255)
                                     tGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACaaa                                                                                                     >  2:524923/1‑150 (MQ=255)
                                       tATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAg                                                                                                   >  1:424920/1‑150 (MQ=255)
                                              gcggGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTc                                                                                            >  1:262680/1‑150 (MQ=255)
                                              gcggGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTc                                                                                            <  3:375226/150‑1 (MQ=255)
                                                  gAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAGCCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCt                                                                                        >  2:389364/1‑150 (MQ=255)
                                                        aGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCAt                                                                                  <  1:137641/150‑1 (MQ=255)
                                                            ccTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGc                                                                              <  1:431390/150‑1 (MQ=255)
                                                             cTTCGCCAGCGTCAGCAACTGCTGGATGCGAATGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAACAGCACTTTCACCTGCTCATCGGcc                                                                             <  2:305170/150‑1 (MQ=255)
                                                                cGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGt                                                                          >  1:387088/1‑150 (MQ=255)
                                                                  ccAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCGCCCTCACGCGCTACGACAACAGCACTTTCACCTGCTCATCGGCCTGttc                                                                        >  4:264133/1‑148 (MQ=255)
                                                                  ccAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTAc                                                                        >  2:208595/1‑150 (MQ=255)
                                                                  ccAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTAc                                                                        <  4:322336/150‑1 (MQ=255)
                                                                     gCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTACCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGa                                                                     <  3:226820/150‑1 (MQ=255)
                                                                       gTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGAcg                                                                   <  1:591261/150‑1 (MQ=255)
                                                                           gCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGt                                                               <  1:284397/150‑1 (MQ=255)
                                                                              aCGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTTCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGaa                                                            <  2:100112/150‑1 (MQ=255)
                                                                              aCGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTTCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGaa                                                            <  2:100020/150‑1 (MQ=255)
                                                                               cGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAAt                                                           <  2:509190/150‑1 (MQ=255)
                                                                                   tGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATgcgc                                                       >  4:9481/1‑150 (MQ=255)
                                                                                        gCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGaa                                                  <  2:378712/150‑1 (MQ=255)
                                                                                           aaGGTTAACCTGGCGGGCGGCATCAACGAGCATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATgg                                               <  1:375808/150‑1 (MQ=255)
                                                                                              gTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCgg                                            >  3:233440/1‑150 (MQ=255)
                                                                                              gTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACCCCAAAAGCACTTTCACCTGCCCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGCCgg                                            >  2:428575/1‑150 (MQ=255)
                                                                                                tAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTa                                          <  4:323676/150‑1 (MQ=255)
                                                                                                tAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTa                                          <  2:240757/150‑1 (MQ=255)
                                                                                                tAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTa                                          <  2:240716/150‑1 (MQ=255)
                                                                                                  aCCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTAtt                                        >  1:432946/1‑150 (MQ=255)
                                                                                                               cATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCaa                           >  4:101257/1‑150 (MQ=255)
                                                                                                               cATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCaa                           <  2:159730/150‑1 (MQ=255)
                                                                                                                      gAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATgcca                    <  1:473230/150‑1 (MQ=255)
                                                                                                                         gATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATgccagcc                 >  1:227172/1‑150 (MQ=255)
                                                                                                                          aTTTTTATCTTGCAGAGCTTGCTGCCCGGGGCAGACCAGATGCCACCATCACGCTATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATgccagcca                <  2:590399/150‑1 (MQ=255)
                                                                                                                          aTTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATgccagcca                >  4:234684/1‑150 (MQ=255)
                                                                                                                          aTTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATgccagcca                >  2:466461/1‑150 (MQ=255)
                                                                                                                           tttttATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATGCCAGccac               <  1:247770/150‑1 (MQ=255)
                                                                                                                               tATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATGCCAGccaccacc           >  1:440850/1‑150 (MQ=255)
                                                                                                                                 tCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATGCCAGCCACCACCaa         <  2:67797/150‑1 (MQ=255)
                                                                                                                                   ttGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCCCATCGGCCTGTACTCACGCGGTGAATGCGCCGGAATGGCGGTATTACTGCCAGGCCAACACGCCAGCCACCACCAAcc       >  4:51429/1‑149 (MQ=255)
                                                                                                                                        caGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATGCCAGCCACCACCAACTCCAcc  <  3:383139/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
AAACCGCTGGGAGCGTCTGTGCGGTCTGCCGGATGAATGTATTCCGGCGGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGGTATTACTGGCAGGTCAACATGCCAGCCACCACCAACTCCACC  >  NC_000913_3_hsa_tpiA/1207089‑1207374

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