Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 4,187,572 G→T M747I (ATG→ATT rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,187,5720GT98.9% 331.8 / NA 91M747I (ATG→ATTrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (51/39);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACC  >  NC_000913_3_hsa_tpiA/4187428‑4187717
                                                                                                                                                 |                                                                                                                                                 
tAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACATATTCGTCAGCTTGCTGGTATGCGTGGTCAGATTGCGaa                                                                                                                                              <  3:639891/150‑1 (MQ=255)
tAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGaa                                                                                                                                              >  3:757244/1‑150 (MQ=255)
      gCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCGACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGa                                                                                                                                        >  3:686619/1‑150 (MQ=255)
       cAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGAt                                                                                                                                       <  2:1209905/150‑1 (MQ=255)
           cTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCt                                                                                                                                   <  2:320074/150‑1 (MQ=255)
            tGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTc                                                                                                                                  >  1:1043052/1‑150 (MQ=255)
             gAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTcc                                                                                                                                 <  1:1334884/150‑1 (MQ=255)
                 ccGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCcatca                                                                                                                             >  2:828888/1‑150 (MQ=255)
                   gTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCatcatc                                                                                                                           >  1:1340485/1‑150 (MQ=255)
                   gTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCatcatc                                                                                                                           >  1:1157582/1‑150 (MQ=255)
                   gTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCatcatc                                                                                                                           <  2:1157582/150‑1 (MQ=255)
                        tAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAAc                                                                                                                      <  3:293085/150‑1 (MQ=255)
                                      aGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGa                                                                                                        <  1:264322/150‑1 (MQ=255)
                                        gAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGCGCGCGTGGTTCTGCGGCACAGACTCGACAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAAc                                                                                                      >  4:124847/1‑150 (MQ=255)
                                             gaAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCg                                                                                                 <  4:55901/150‑1 (MQ=255)
                                              aaGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACcgc                                                                                                          >  2:243094/1‑140 (MQ=255)
                                              aaGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACcgc                                                                                                          <  1:243094/140‑1 (MQ=255)
                                              aaGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGt                                                                                                <  4:367132/150‑1 (MQ=255)
                                                 cAGGTCTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGaa                                                                                             <  1:749545/150‑1 (MQ=255)
                                                  aGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                            >  2:830325/1‑150 (MQ=255)
                                                    gTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGt                                                                                          >  4:452277/1‑150 (MQ=255)
                                                          ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                                                                    >  2:717853/1‑150 (MQ=255)
                                                            cAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGt                                                                                  <  2:159196/150‑1 (MQ=255)
                                                             aaCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGaa                                                                                                                                              <  1:461562/89‑1 (MQ=255)
                                                             aaCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGaa                                                                                                                                              >  2:461562/1‑89 (MQ=255)
                                                              aCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTAc                                                                                >  3:572433/1‑150 (MQ=255)
                                                               cAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACt                                                                               >  1:518562/1‑150 (MQ=255)
                                                               cAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACt                                                                               <  1:891678/150‑1 (MQ=255)
                                                                 gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTc                                                                                         <  3:226148/138‑1 (MQ=255)
                                                                 gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTc                                                                                         >  4:226148/1‑138 (MQ=255)
                                                                    tCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGaa                                                                                                                                              >  3:666398/1‑82 (MQ=255)
                                                                    tCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGaa                                                                                                                                              <  4:666398/82‑1 (MQ=255)
                                                                       aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACt                                                                       >  1:647540/1‑150 (MQ=255)
                                                                        cATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTcc                                                                             >  2:239976/1‑143 (MQ=255)
                                                                        cATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTcc                                                                             <  1:239976/143‑1 (MQ=255)
                                                                        cATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACtt                                                                      <  1:813571/150‑1 (MQ=255)
                                                                         atgatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTc                                                                     >  1:824351/1‑150 (MQ=255)
                                                                           gatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACg                                                                                   >  1:1253321/1‑134 (MQ=255)
                                                                           gatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACg                                                                                   <  2:1253321/134‑1 (MQ=255)
                                                                            atgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtc                                                                  <  2:818558/150‑1 (MQ=255)
                                                                              ggCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtctc                                                                >  1:670513/1‑150 (MQ=255)
                                                                                ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcca                                                              >  1:134887/1‑150 (MQ=255)
                                                                                ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcca                                                              <  4:151778/150‑1 (MQ=255)
                                                                                ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcca                                                              >  3:195250/1‑150 (MQ=255)
                                                                                 cGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccac                                                             >  4:98577/1‑150 (MQ=255)
                                                                                     tCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccacccac                                                         >  2:345927/1‑150 (MQ=255)
                                                                                       cGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcc                                                               <  2:522808/142‑1 (MQ=255)
                                                                                       cGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcc                                                               >  1:522808/1‑142 (MQ=255)
                                                                                        tatGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGt                                                      <  3:439153/148‑1 (MQ=255)
                                                                                         gTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCGTTCTTTTATGCTGGTGCTTATTTCGAGGCCGGAGTGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCTTGATGATCTGAACGTACTCCAGCACTTCATCGCCACCCACGGTg                                                     >  2:1356248/1‑150 (MQ=255)
                                                                                            cgcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTc                                                  >  2:109437/1‑150 (MQ=255)
                                                                                             gcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCg                                                                                                                                          <  1:772496/61‑1 (MQ=255)
                                                                                             gcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCg                                                                                                                                          >  2:772496/1‑61 (MQ=255)
                                                                                             gcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTACATCTCCACCCACGGTGCTCg                                                 >  1:694740/1‑150 (MQ=255)
                                                                                               gTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTa                                               >  4:569659/1‑150 (MQ=255)
                                                                                                  gTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAg                                            >  2:209212/1‑150 (MQ=255)
                                                                                                      tGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACcgc                                                                                                          >  1:1176832/1‑84 (MQ=255)
                                                                                                      tGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACcgc                                                                                                          <  2:1176832/84‑1 (MQ=255)
                                                                                                      tGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCt                                        >  4:264540/1‑150 (MQ=255)
                                                                                                          gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGTGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCg                                    <  1:828888/150‑1 (MQ=255)
                                                                                                          gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGaa                                                                                                       <  2:948954/83‑1 (MQ=255)
                                                                                                          gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGaa                                                                                                       >  1:948954/1‑83 (MQ=255)
                                                                                                            acaGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGa                                  <  2:1043052/150‑1 (MQ=255)
                                                                                                              aGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGAt                                                                                                                                       >  1:692841/1‑47 (MQ=255)
                                                                                                              aGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGAt                                                                                                                                       <  2:692841/47‑1 (MQ=255)
                                                                                                               gATTCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATAc                               <  3:283314/150‑1 (MQ=255)
                                                                                                                 ttCGTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCg                             <  3:98577/150‑1 (MQ=255)
                                                                                                                    gTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCAc                          <  3:452277/150‑1 (MQ=255)
                                                                                                                    gTCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCAc                          >  2:229058/1‑150 (MQ=255)
                                                                                                                     tCAGCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACt                         >  3:540066/1‑150 (MQ=255)
                                                                                                                        gCTTGCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGaa                      >  2:968183/1‑150 (MQ=255)
                                                                                                                            gCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACt                  <  3:124847/150‑1 (MQ=255)
                                                                                                                            gCTGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACt                  >  2:1280198/1‑150 (MQ=255)
                                                                                                                              tGGTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGc                >  1:114156/1‑150 (MQ=255)
                                                                                                                               tgTATGCGTGGGCTGATTGCGAAGCCGGATTGCTCCATCAGCGAAACGCCATTCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCg               <  2:336928/149‑1 (MQ=255)
                                                                                                                                gTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCcat                                                                                                                               <  2:736323/37‑1 (MQ=255)
                                                                                                                                gTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCcat                                                                                                                               >  1:736323/1‑37 (MQ=255)
                                                                                                                                gTATGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGa              >  2:900053/1‑150 (MQ=255)
                                                                                                                                  aagCGTGGTCTGATTGCGAAGCCGGATGGCTCAATCATCGAAACACCAATCACCGCGAACTTACGTGACCGTATGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAAc            >  2:961121/3‑150 (MQ=255)
                                                                                                                                  aTGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTg                                                     <  2:1241774/109‑1 (MQ=255)
                                                                                                                                  aTGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTg                                                     >  1:1241774/1‑109 (MQ=255)
                                                                                                                                   tGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCa                                                                                                                  <  1:671399/47‑1 (MQ=255)
                                                                                                                                   tGCGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCa                                                                                                                  >  2:671399/1‑47 (MQ=255)
                                                                                                                                     cGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTc                                                                     <  1:487473/90‑1 (MQ=255)
                                                                                                                                     cGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTc                                                                     >  2:487473/1‑90 (MQ=255)
                                                                                                                                     cGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTcc         >  3:67268/1‑150 (MQ=255)
                                                                                                                                      gTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCAc                                                                                                             >  3:501431/1‑49 (MQ=255)
                                                                                                                                      gTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCAc                                                                                                             <  4:501431/49‑1 (MQ=255)
                                                                                                                                      gTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCg        >  1:244886/1‑150 (MQ=255)
                                                                                                                                            cTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAc   <  4:435227/150‑1 (MQ=255)
                                                                                                                                             tGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAcc  >  1:1333380/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACC  >  NC_000913_3_hsa_tpiA/4187428‑4187717

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