Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 3,229,001 G→A P38P (CCG→CCA ygjK → alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA3,229,0010GA95.7% 215.8 / 0.1 69P38P (CCG→CCAygjKalpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (28/38);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

ATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGAAGAGTACA  >  NC_000913_3_hsa_tpiA/3228856‑3229146
                                                                                                                                                 |                                                                                                                                                 
atatTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGt                                                                                                                                               <  2:296009/150‑1 (MQ=255)
  atTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTAc                                                                                                                                             <  1:80315/150‑1 (MQ=255)
    ttCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCTAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACAt                                                                                                                                           <  2:498239/150‑1 (MQ=255)
              aataaGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGgattac                                                                                                                                 <  4:56534/150‑1 (MQ=255)
              aataaGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGgattac                                                                                                                                 <  3:369528/150‑1 (MQ=255)
                 aaGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGgattacgat                                                                                                                              >  2:159449/1‑150 (MQ=255)
                       tGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTacgac                                                                                                                        >  4:41321/1‑150 (MQ=255)
                       tGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTGTTTCCGCCCACGCCGCTAACGCCGGCAATTATAAAAACGTGATTAACCGCACTCCCGCGCCACAGTCCAGGAAGGATCACCACTaccac                                                                                                                        >  2:762534/1‑150 (MQ=255)
                          aaTTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTTATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGAt                                                                                                                     >  3:198247/1‑150 (MQ=255)
                          aaTTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGAt                                                                                                                     <  4:268315/150‑1 (MQ=255)
                            tttCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCa                                                                                                                   >  2:145180/1‑150 (MQ=255)
                            tttCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCCGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCa                                                                                                                   >  2:255967/1‑150 (MQ=255)
                                aTGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAg                                                                                                               >  4:135303/1‑150 (MQ=255)
                                  gAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCg                                                                                                             <  1:758456/150‑1 (MQ=255)
                                     aaTAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGttt                                                                                                          <  3:266041/150‑1 (MQ=255)
                                         aaaaCTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                                                                                                      <  4:198247/150‑1 (MQ=255)
                                           aaCTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATcc                                                                                                    >  4:143420/1‑150 (MQ=255)
                                              tATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGtt                                                                                                 <  4:457428/150‑1 (MQ=255)
                                               aTTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGttt                                                                                                >  2:263368/1‑150 (MQ=255)
                                                   tAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCg                                                                                            <  2:270133/150‑1 (MQ=255)
                                                      cGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtc                                                                                         >  2:278916/1‑150 (MQ=255)
                                                        ccAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc                                                                                       <  1:278916/150‑1 (MQ=255)
                                                         cAGTAACCTGCGCTCCGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCg                                                                                      <  1:255967/150‑1 (MQ=255)
                                                           gTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGa                                                                                    <  1:145180/150‑1 (MQ=255)
                                                            tAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAg                                                                                   <  3:176463/150‑1 (MQ=255)
                                                              aCCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGcc                                                                                 >  2:6295/1‑150 (MQ=255)
                                                                cTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg                                                                               >  2:50606/1‑150 (MQ=255)
                                                                 tGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTgg                                                                              <  2:191068/150‑1 (MQ=255)
                                                                 tGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTgg                                                                              >  2:288978/1‑150 (MQ=255)
                                                                  gcgcTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGc                                                                             >  1:542373/1‑150 (MQ=255)
                                                                     ctctGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATg                                                                          <  2:86705/150‑1 (MQ=255)
                                                                        tgctgATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTc                                                                       >  2:174897/1‑150 (MQ=255)
                                                                                   ttttCCGCCCATGCCGCTAACGCCGACAAGTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCCCGGACCCTGGCATGGTCATCTGTTGccc                                                            >  4:192796/1‑149 (MQ=255)
                                                                                          cccATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGcc                                                     <  3:143420/150‑1 (MQ=255)
                                                                                          cccATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGcc                                                     <  2:51220/150‑1 (MQ=255)
                                                                                              tGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTaa                                                 <  4:207311/150‑1 (MQ=255)
                                                                                              tGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTaa                                                 <  3:41321/150‑1 (MQ=255)
                                                                                               gCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAac                                                <  2:698359/150‑1 (MQ=255)
                                                                                                ccGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAaca                                               <  4:208296/150‑1 (MQ=255)
                                                                                                ccGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAaca                                               <  3:135303/150‑1 (MQ=255)
                                                                                                ccGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTAAATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAaca                                               <  1:159449/150‑1 (MQ=255)
                                                                                                     aaCGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATg                                          <  1:709946/150‑1 (MQ=255)
                                                                                                     aaCGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATg                                          <  1:709981/150‑1 (MQ=255)
                                                                                                        gCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTc                                                                       <  2:485705/118‑1 (MQ=255)
                                                                                                        gCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTc                                                                       >  1:485705/1‑118 (MQ=255)
                                                                                                         ccGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGggcg                                      >  4:121886/1‑150 (MQ=255)
                                                                                                          cGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGc                                                                             >  1:380253/1‑110 (MQ=255)
                                                                                                          cGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGc                                                                             <  2:380253/110‑1 (MQ=255)
                                                                                                          cGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGggcgg                                     >  1:725827/1‑150 (MQ=255)
                                                                                                              aaTTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGgatta                                                                                                                                  <  3:92173/53‑1 (MQ=255)
                                                                                                              aaTTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGgatta                                                                                                                                  >  4:92173/1‑53 (MQ=255)
                                                                                                              aaTTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc                                                                                       >  1:670641/1‑96 (MQ=255)
                                                                                                              aaTTATAAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc                                                                                       <  2:670641/96‑1 (MQ=255)
                                                                                                                   taAAAACGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCggg                            >  1:126307/1‑150 (MQ=255)
                                                                                                                     aaaaCGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCg                          <  1:342438/150‑1 (MQ=255)
                                                                                                                       aaCGTGATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGtt                        <  2:47615/150‑1 (MQ=255)
                                                                                                                            gATTAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCgct                   <  2:542373/150‑1 (MQ=255)
                                                                                                                              ttAACCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCgctgc                 >  3:162616/1‑150 (MQ=255)
                                                                                                                                aaCCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGctgctg               >  1:660851/1‑150 (MQ=255)
                                                                                                                                 aCCGTACTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGa              >  3:360568/1‑150 (MQ=255)
                                                                                                                                      aCTGGCGCGCCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGaa         <  1:164050/150‑1 (MQ=255)
                                                                                                                                         ggCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGAagag      >  1:515261/1‑150 (MQ=255)
                                                                                                                                          gcgcgcCACAGTACATGAAGGATTACGATTACGACGATCACCa                                                                                                                >  1:599201/1‑43 (MQ=255)
                                                                                                                                          gcgcgcCACAGTACATGAAGGATTACGATTACGACGATCACCa                                                                                                                <  2:599201/43‑1 (MQ=255)
                                                                                                                                          gcgcgcCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTc                                                                                             <  2:760951/62‑1 (MQ=255)
                                                                                                                                          gcgcgcCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTc                                                                                             >  1:760951/1‑62 (MQ=255)
                                                                                                                                          gcgcgcCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATACGCTGTTTCATCTAGGCTGCTGGCATCGTGATCTTGGGCCAGACGGTCCTTACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGAAGAGt     >  1:360752/1‑150 (MQ=255)
                                                                                                                                           cgcgcCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGAAGAGTa    <  1:174897/150‑1 (MQ=255)
                                                                                                                                            gcgcCACAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGAAGAGTAc   <  2:725827/150‑1 (MQ=255)
                                                                                                                                             cgccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGAAGAGTACa  <  1:618706/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATATTTCTACTGGTAATAAGGACTGTAATTTCATGAAAATAAAAACTATTTTAACGCCAGTAACCTGCGCTCTGCTGATAAGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGCATGGTCATCTGTTGCCAGACGGCCCTAACACCATGGGCGGCTTTCCGGGCGTTGCGCTGCTGACGGAAGAGTACA  >  NC_000913_3_hsa_tpiA/3228856‑3229146

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: