Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 4,546,943 C→A S622R (AGC→AGA fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,546,9430CA99.1% 372.2 / NA 106S622R (AGC→AGAfimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (47/57);  total (48/58)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

AACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACG  >  NC_000913_3_hsa_tpiA/4546797‑4547087
                                                                                                                                                  |                                                                                                                                                
aaCTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATg                                                                                                                                               >  3:130210/1‑150 (MQ=255)
 aCTGGACGCTCAGCTATAGCCTGACGAGAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGt                                                                                                                                              >  2:383225/1‑150 (MQ=255)
   tGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTGCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTca                                                                                                                                            <  1:1057527/150‑1 (MQ=255)
   tGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTca                                                                                                                                            >  1:1199548/1‑150 (MQ=255)
   tGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTca                                                                                                                                            >  1:729080/1‑150 (MQ=255)
     gACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTcaca                                                                                                                                          >  2:568124/1‑150 (MQ=255)
      aCGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACATCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCacac                                                                                                                                         <  3:613320/150‑1 (MQ=255)
       cGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACg                                                                                                                                        <  4:729166/150‑1 (MQ=255)
        gCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTGCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTTCCAGCTACAGAATGTCACACGa                                                                                                                                       >  4:72282/1‑150 (MQ=255)
         cTCAGCTATAGCCTGACGAAAACCGCCTGGCAAAAAGGCCGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGAt                                                                                                                                      <  3:219360/150‑1 (MQ=255)
          tctgCTATAGCTTGACGAAAATCTCCTGGCAAAAAGGACGTGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGATAGTGCCTGCTACTGTTTGTCTCACGAtc                                                                                                                                     <  4:130210/147‑1 (MQ=255)
          tCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGAtc                                                                                                                                     >  4:509377/1‑150 (MQ=255)
          tCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGAtc                                                                                                                                     >  1:749970/1‑150 (MQ=255)
                ataGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACg                                                                                                                               >  4:534318/1‑150 (MQ=255)
                ataGCCTGACGAAAAACGACTGGAAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACg                                                                                                                               <  2:749063/150‑1 (MQ=255)
                    ccTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGCGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCg                                                                                                                           >  2:33790/1‑150 (MQ=255)
                     cTGACGAAGAACGCCGGGCAAAAAGCCCGGGATCAGATGTGAGCGCGTAATGTCAATATTCCTGTCAGCCACTGGCTGCGTTCTGACAGTAAATCGCAGTGGCGACATGCCAGTGCCCGCTACAGAATGTCACACGATCTCAACGGTCgg                                                                                                                          <  2:360766/150‑1 (MQ=255)
                        aCGAAAAACGCCAGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATg                                                                                                                       <  2:318410/150‑1 (MQ=255)
                         cGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGa                                                                                                                      <  1:568124/150‑1 (MQ=255)
                          gAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGAc                                                                                                                     <  1:745281/150‑1 (MQ=255)
                             aaaCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCaa                                                                                                                  >  4:249180/1‑150 (MQ=255)
                             aaaCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCGGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCaa                                                                                                                  >  3:72700/1‑150 (MQ=255)
                              aaCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAAt                                                                                                                 <  1:683745/150‑1 (MQ=255)
                               aCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGAt                                                                                                                        <  3:170118/142‑1 (MQ=255)
                               aCGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGAt                                                                                                                        >  4:170118/1‑142 (MQ=255)
                                   cTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATctggc                                                                                                            <  3:167737/150‑1 (MQ=255)
                                       cAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGgt                                                                                                        <  3:413561/150‑1 (MQ=255)
                                         aagaGGACGTGATCAGAGGTTAGCGCTATACGTCAATATTCCATTCAGCAACCGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCACAACGGTCGGATGCCCAATCTGGCTGgtgt                                                                                                      <  2:639802/147‑1 (MQ=255)
                                         aaaaGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGgtgt                                                                                                      >  1:420174/1‑150 (MQ=255)
                                                   gATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACg                                                                                            >  2:683800/1‑150 (MQ=255)
                                                   gATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACg                                                                                            >  2:683805/1‑150 (MQ=255)
                                                   gATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGCGTATACGGTACg                                                                                            >  2:1049017/1‑150 (MQ=255)
                                                     tCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGtt                                                                                          >  1:37882/1‑150 (MQ=255)
                                                        gATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCt                                                                                       <  2:397140/150‑1 (MQ=255)
                                                        gATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCt                                                                                       >  4:668995/1‑150 (MQ=255)
                                                        gATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACACGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCt                                                                                       <  1:1020286/150‑1 (MQ=255)
                                                         aTGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTg                                                                                      <  2:25717/150‑1 (MQ=255)
                                                            ttAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGaa                                                                                   <  4:709642/150‑1 (MQ=255)
                                                            ttAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGaa                                                                                   <  1:767398/150‑1 (MQ=255)
                                                            ttAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGaa                                                                                   <  3:534318/150‑1 (MQ=255)
                                                            ttAGCGCTTAACGCCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGaa                                                                                   <  1:615987/150‑1 (MQ=255)
                                                             tAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACCTGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAg                                                                                  <  4:32381/150‑1 (MQ=255)
                                                              aGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGa                                                                                 <  1:988521/150‑1 (MQ=255)
                                                               gcgcTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGac                                                                                <  2:440855/150‑1 (MQ=255)
                                                               gcgcTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGac                                                                                >  1:220188/1‑150 (MQ=255)
                                                                   ttAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGacaaca                                                                            <  1:203445/150‑1 (MQ=255)
                                                                      aCGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAAt                                                                                                                                                <  3:521881/79‑1 (MQ=255)
                                                                        gTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTc                                                                       >  1:300974/1‑150 (MQ=255)
                                                                        gTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTc                                                                       >  2:826969/1‑150 (MQ=255)
                                                                               ttCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGgtg                                                                                                       <  4:506712/111‑1 (MQ=255)
                                                                               ttCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGgtg                                                                                                       >  3:506712/1‑111 (MQ=255)
                                                                               ttCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCATCtata                                                                <  4:632646/150‑1 (MQ=255)
                                                                                 ccTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGc                                                              <  2:1170412/150‑1 (MQ=255)
                                                                                  cTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCTGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCg                                                             >  1:831188/1‑150 (MQ=255)
                                                                                   tttCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGc                                                              >  4:367779/1‑148 (MQ=255)
                                                                                         ccACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGgt                                                                                                        <  2:558895/100‑1 (MQ=255)
                                                                                         ccACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGgt                                                                                                        >  1:558895/1‑100 (MQ=255)
                                                                                         ccACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAAc                                                      >  1:549111/1‑150 (MQ=255)
                                                                                          cACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGac                                                                                <  1:295482/123‑1 (MQ=255)
                                                                                          cACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGac                                                                                >  2:295482/1‑123 (MQ=255)
                                                                                            cTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCgg                                                   >  4:46291/1‑150 (MQ=255)
                                                                                            cTGGCTGCGTTCTGACAGTAAATCTCAGCGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCgg                                                   <  1:33790/150‑1 (MQ=255)
                                                                                             tGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGc                                                  >  3:753186/1‑150 (MQ=255)
                                                                                               gCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTc                                                                       >  1:1135081/1‑127 (MQ=255)
                                                                                               gCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTc                                                                       <  2:1135081/127‑1 (MQ=255)
                                                                                                     ttCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCgg                                          >  3:200226/1‑150 (MQ=255)
                                                                                                      tCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCggg                                         <  1:965369/150‑1 (MQ=255)
                                                                                                      tCTGACAGTAAATCTCAGTGGCGACATGCCAGGGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCggg                                         <  2:1199548/150‑1 (MQ=255)
                                                                                                        tGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCggggg                                       >  4:660011/1‑150 (MQ=255)
                                                                                                         gACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAGCCGGCTATGCCGGGGGa                                      <  4:72700/150‑1 (MQ=255)
                                                                                                          aCAGTAACTCACAGTGGCGACATGCCAGTGTCAGCAACAGACTGTCATACGATCTCAACGGTCGGATGGCCAATCTGGCTGGTGGATACGGTACGTTGCTGGAAGACAACAACCTCGGCTATAGCGTGCAAACCGGCTATGCCGGGGGAg                                     <  2:1068892/150‑1 (MQ=255)
                                                                                                          aCAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATctgg                                                                                                             <  4:658862/78‑1 (MQ=255)
                                                                                                          aCAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATctgg                                                                                                             >  3:658862/1‑78 (MQ=255)
                                                                                                          aCAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAg                                     >  3:320176/1‑150 (MQ=255)
                                                                                                                aaTCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAGCCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATg                               <  2:884798/150‑1 (MQ=255)
                                                                                                                aaTCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATg                               <  2:471724/150‑1 (MQ=255)
                                                                                                                 aTCTCAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATgg                              >  3:43757/1‑150 (MQ=255)
                                                                                                                    tcAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGaaa                           <  1:1029609/150‑1 (MQ=255)
                                                                                                                    tcAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGaaa                           >  3:600627/1‑150 (MQ=255)
                                                                                                                    tcAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGaaa                           >  2:268015/1‑150 (MQ=255)
                                                                                                                     cAGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGaaaa                          >  1:19596/1‑149 (MQ=255)
                                                                                                                      aGTGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATa                         <  3:249180/150‑1 (MQ=255)
                                                                                                                        tGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCggg                                         >  4:120245/1‑132 (MQ=255)
                                                                                                                        tGGCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCggg                                         <  3:120245/132‑1 (MQ=255)
                                                                                                                          gCGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCgg                     <  4:41874/150‑1 (MQ=255)
                                                                                                                           cGACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGGGGCGATGGAAATAGCGGa                    >  3:197097/1‑150 (MQ=255)
                                                                                                                            gACATGCCAGTGCCAGCTACAGAATGTCTCACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGaa                   <  3:685269/150‑1 (MQ=255)
                                                                                                                            gACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGaa                   >  3:764834/1‑150 (MQ=255)
                                                                                                                            gACATGCCAGTGCCAGCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGaa                   <  2:812380/150‑1 (MQ=255)
                                                                                                                                 gccagtgccagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACa              >  3:177505/1‑150 (MQ=255)
                                                                                                                                 gccagtgccagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACa              <  2:220188/150‑1 (MQ=255)
                                                                                                                                 gccagtgccagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACa              >  2:1063568/1‑150 (MQ=255)
                                                                                                                                  ccagtgccagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAg             >  3:306213/1‑150 (MQ=255)
                                                                                                                                   cagtgccagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAgg            <  1:1063568/150‑1 (MQ=255)
                                                                                                                                        ccagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCgg                                                   >  3:28688/1‑106 (MQ=255)
                                                                                                                                         cagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGc      >  1:836441/1‑150 (MQ=255)
                                                                                                                                         cagCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGc      <  2:300974/150‑1 (MQ=255)
                                                                                                                                           gCTACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCa    <  3:46291/150‑1 (MQ=255)
                                                                                                                                             tACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACg  >  4:28806/1‑150 (MQ=255)
                                                                                                                                             tACAGAATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACg  <  3:348210/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AACTGGACGCTCAGCTATAGCCTGACGAAAAACGCCTGGCAAAAAGGACGGGATCAGATGTTAGCGCTTAACGTCAATATTCCTTTCAGCCACTGGCTGCGTTCTGACAGTAAATCTCAGTGGCGACATGCCAGTGCCAGCTACAGCATGTCACACGATCTCAACGGTCGGATGACCAATCTGGCTGGTGTATACGGTACGTTGCTGGAAGACAACAACCTCAGCTATAGCGTGCAAACCGGCTATGCCGGGGGAGGCGATGGAAATAGCGGAAGTACAGGCTACGCCACG  >  NC_000913_3_hsa_tpiA/4546797‑4547087

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: