Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 1,158,559 2 bp→TC coding (691‑692/1434 nt) ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,158,5590CT92.6% 623.6 / 26.3 203Q231S (CAG→TCG) ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base C (7/8);  new base T (100/88);  total (107/96)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01
*NC_000913_3_hsa_tpiA1,158,5600AC92.6% 686.5 / 31.3 203Q231S (CAG→TCG) ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base A (7/8);  new base C (100/88);  total (107/96)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 205 total aligned reads displayed.

CAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGT  >  NC_000913_3_hsa_tpiA/1158418‑1158707
                                                                                                                                             ||                                                                                                                                                   
cAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTATGCACCACATCTCGAACGTACCGTTCCGGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTc                                                                                                                                              <  2:1399595/150‑1 (MQ=255)
cAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTAATCGCCCGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTc                                                                                                                                              <  3:1432388/150‑1 (MQ=255)
 aGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTcc                                                                                                                                             >  1:576856/1‑150 (MQ=255)
     ggCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACgg                                                                                                                                         >  1:1314771/1‑150 (MQ=255)
     ggCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACgg                                                                                                                                         <  4:1103602/150‑1 (MQ=255)
              ccAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTcc                                                                                                                                <  1:1295968/150‑1 (MQ=255)
               cAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCcg                                                                                                                               >  3:258980/1‑150 (MQ=255)
               cAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCcg                                                                                                                               >  4:1092268/1‑150 (MQ=255)
                  aaCCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGt                                                                                                                            <  1:76255/150‑1 (MQ=255)
                  aaCCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGt                                                                                                                            <  3:1125226/150‑1 (MQ=255)
                    cccGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTta                                                                                                                          <  3:1114939/150‑1 (MQ=255)
                     ccGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTtat                                                                                                                         <  2:452547/150‑1 (MQ=255)
                       ggTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTatat                                                                                                                       <  2:2012005/150‑1 (MQ=255)
                       ggTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTatat                                                                                                                       <  4:700901/150‑1 (MQ=255)
                        gTAGTTGCGTTTGGCATTTACGATTTCATCGAACGTTGCCTGGTACCGTTTGGTCAGCACCACATCTGAAACGTACCTTTCCAGATGCAGATAGGTGAATACACCAATGCATCAGGTTCGGTATTCAACCGCAACATTCCGCGTTTTATg                                                                                                                      >  3:68188/1‑150 (MQ=255)
                         tAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATgg                                                                                                                     <  1:549994/150‑1 (MQ=255)
                            ttGTGTTTGGCATTTACGGTTTGATCGAACGTTGCATGGTACCGTTTGGTCTTCACCTGATCTGGAACTTACCTTTCCTGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCTGTTTTCCACGGCGACATTCCGCGTTTTATGGCgg                                                                                                                  <  4:786610/150‑1 (MQ=255)
                            ttGCGTTTGGCATTTTCGGTTTCATCGAACGTTGCTTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCgg                                                                                                                  <  4:515176/150‑1 (MQ=255)
                            ttGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGttt                                                                                                                                                >  3:1474791/1‑120 (MQ=255)
                            ttGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCgg                                                                                                                  >  1:1536664/1‑150 (MQ=255)
                                       aTTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACt                                                                                                       >  1:1866849/1‑150 (MQ=255)
                                         ttACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAGCGTACCTTTCCGGTTGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGATATTCCGCGTTATATGGCGGGTGACCCGACTGc                                                                                                     <  2:112518/150‑1 (MQ=255)
                                          tACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCg                                                                                                    <  3:915357/150‑1 (MQ=255)
                                           atggTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTc                                                                                                                                              <  1:21902/105‑1 (MQ=255)
                                           aCGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCgg                                                                                                   <  2:98194/150‑1 (MQ=255)
                                           aCGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCgg                                                                                                   <  4:650565/150‑1 (MQ=255)
                                             ggTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGt                                                                                                 >  1:1788390/1‑150 (MQ=255)
                                               tttCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaa                                                                                               >  1:1597411/1‑150 (MQ=255)
                                                    tCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGt                                                                                          >  3:1238724/1‑150 (MQ=255)
                                                    tCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGt                                                                                          >  3:397386/1‑150 (MQ=255)
                                                    tCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGt                                                                                          >  3:528679/1‑150 (MQ=255)
                                                     cGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTc                                                                                         <  1:1441453/150‑1 (MQ=255)
                                                      gAACGCTGCCTGGTACCGCTTGGTCTGCACCACATCTGGAACGTACCTTTCCATATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTCATATGGCGGGTGACACGACTGCGGGTACACTGTCt                                                                                        >  2:1823477/1‑150 (MQ=255)
                                                       aaCGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGa                                                                                                              >  4:1278878/1‑127 (MQ=255)
                                                          gTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtggtg                                                                                    <  1:1221205/150‑1 (MQ=255)
                                                          gTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtggtg                                                                                    <  4:397386/150‑1 (MQ=255)
                                                           ttGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtggtgg                                                                                   >  1:694597/1‑150 (MQ=255)
                                                             gCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGt                                                                                                                <  4:1394123/119‑1 (MQ=255)
                                                             gCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGt                                                                                                                >  3:1394123/1‑119 (MQ=255)
                                                               cTGGTTCCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTc                                                                               <  4:301356/150‑1 (MQ=255)
                                                                           ggTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTtat                                                                                                                         >  2:1710241/1‑96 (MQ=255)
                                                                           ggTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTtat                                                                                                                         <  1:1710241/96‑1 (MQ=255)
                                                                            gTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTGTGGTGGCTTCATGTTCaacgacc                                                                  >  4:1358262/1‑145 (MQ=255)
                                                                            gTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaaaa                                                                     >  4:1173922/1‑147 (MQ=255)
                                                                             tCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCgg                                                                                                                  >  1:1691237/1‑101 (MQ=255)
                                                                             tCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCgg                                                                                                                  <  2:1691237/101‑1 (MQ=255)
                                                                              cTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTAc                                                                >  3:604694/1‑150 (MQ=255)
                                                                               tGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACg                                                               >  2:2422572/1‑150 (MQ=255)
                                                                                gCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCg                                                                                                                   >  4:520149/1‑97 (MQ=255)
                                                                                 caccacATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGt                                                             >  3:614434/1‑150 (MQ=255)
                                                                                  accacATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACgg                                                                                                                                         <  1:1577508/73‑1 (MQ=255)
                                                                                  accacATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACgg                                                                                                                                         >  2:1577508/1‑73 (MQ=255)
                                                                                    cacaTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCgg                                                                                                                  >  3:984965/1‑94 (MQ=255)
                                                                                     acaTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGc                                                         <  4:258980/150‑1 (MQ=255)
                                                                                      caTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCt                                                                             <  3:753183/129‑1 (MQ=255)
                                                                                         cTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGc                                                     >  4:738488/1‑150 (MQ=255)
                                                                                            gAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGc                                                  <  4:618856/150‑1 (MQ=255)
                                                                                            gAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGc                                                  <  4:19697/150‑1 (MQ=255)
                                                                                            gAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGc                                                  <  1:2114316/150‑1 (MQ=255)
                                                                                             aaCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaaaa                                                                     <  3:627767/130‑1 (MQ=255)
                                                                                             aaCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaaaa                                                                     >  4:627767/1‑130 (MQ=255)
                                                                                              aCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTtata                                                                                                                        <  2:988931/78‑1 (MQ=255)
                                                                                              aCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTtata                                                                                                                        >  1:988931/1‑78 (MQ=255)
                                                                                              aCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCg                                                <  2:2247779/150‑1 (MQ=255)
                                                                                              aCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCgg                                                                                                   >  4:1414566/1‑99 (MQ=255)
                                                                                               cGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGc                                               <  2:236871/150‑1 (MQ=255)
                                                                                                gTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtggt                                                                                     <  2:699677/111‑1 (MQ=255)
                                                                                                gTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCa                                              >  3:1024034/1‑150 (MQ=255)
                                                                                                 tACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCaa                                             <  3:729987/150‑1 (MQ=255)
                                                                                                   ccTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAAtt                                           <  2:236003/150‑1 (MQ=255)
                                                                                                     tttCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGa                                                                                                                                      >  2:2063707/1‑57 (MQ=255)
                                                                                                     tttCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTg                                          >  4:102745/1‑149 (MQ=255)
                                                                                                     tttCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGc                                         <  1:2422572/150‑1 (MQ=255)
                                                                                                      ttCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCt                                        >  4:791976/1‑150 (MQ=255)
                                                                                                              gCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTa                                                                 <  1:2298073/117‑1 (MQ=255)
                                                                                                              gCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGc                                               >  3:428702/1‑135 (MQ=255)
                                                                                                              gCAGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGAGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCa                                <  3:1092268/150‑1 (MQ=255)
                                                                                                                aGATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCAct                              <  2:1974109/150‑1 (MQ=255)
                                                                                                                aGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGa                                                                                                                                      <  1:2016148/46‑1 (MQ=255)
                                                                                                                 gATTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCActc                             >  1:2363659/1‑150 (MQ=255)
                                                                                                                  aTTGGTGAATACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCActct                            <  4:604694/150‑1 (MQ=255)
                                                                                                                         aaTACACCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAAc                     >  3:1141745/1‑150 (MQ=255)
                                                                                                                            acacCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTc                                                                               >  2:1858351/1‑89 (MQ=255)
                                                                                                                             cacCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCgg                                                                                                   <  1:1832783/68‑1 (MQ=255)
                                                                                                                              acCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTg                                   >  3:212350/1‑131 (MQ=255)
                                                                                                                              acCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTg                                   <  4:212350/131‑1 (MQ=255)
                                                                                                                              acCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGaa                >  3:618916/1‑150 (MQ=255)
                                                                                                                              acCAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGaa                >  1:876107/1‑150 (MQ=255)
                                                                                                                               ccAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGaaa               >  1:25435/1‑150 (MQ=255)
                                                                                                                                cAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa                                                                                              >  1:2307618/1‑70 (MQ=255)
                                                                                                                                cAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa                                                                                              <  2:2307618/70‑1 (MQ=255)
                                                                                                                                cAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGaaaa              <  2:633436/150‑1 (MQ=255)
                                                                                                                                cAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGaaaa              <  2:408309/150‑1 (MQ=255)
                                                                                                                                cAACGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGaaaa              >  1:633436/1‑150 (MQ=255)
                                                                                                                                 aaCGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAAc             >  2:874063/1‑150 (MQ=255)
                                                                                                                                 aaCGCAGCAGGTTCGGTTTTCCACGGCGACATTCCGCGTTACATGGCGGGTGACCAGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATTTACGGTCTGCCAGCTGCCACAATTGCTATCTGGCACTCTGCTAAACCAGAacac             >  2:123056/1‑150 (MQ=255)
                                                                                                                                     cagcagGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACcgcg         >  3:459321/1‑150 (MQ=255)
                                                                                                                                      agcagGTTCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACcgcgc        >  4:911134/1‑150 (MQ=255)
                                                                                                                                           gtcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAg   >  1:944210/1‑150 (MQ=255)
                                                                                                                                            tcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGt  >  1:1450397/1‑150 (MQ=255)
                                                                                                                                             ||                                                                                                                                                   
CAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGT  >  NC_000913_3_hsa_tpiA/1158418‑1158707

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