Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 654,773 T→C 6.6% Q399R (CAG→CGG)  dcuC ← anaerobic C4‑dicarboxylate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA654,7730TC 6.6% 186.4 / 3.6 62Q399R (CAG→CGG) dcuCanaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base C (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGG  >  NC_000913_3_hsa_tpiA/654625‑654922
                                                                                                                                                    |                                                                                                                                                     
ctctGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTg                                                                                                                                                      <  1:115130/150‑1 (MQ=255)
  ctGTAGCAACGATAACAATCACCAACCCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAGCGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTgca                                                                                                                                                    <  2:59423/150‑1 (MQ=255)
  ctGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTgca                                                                                                                                                    <  4:238539/150‑1 (MQ=255)
    gTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCGTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTgcagc                                                                                                                                                  >  3:149847/1‑150 (MQ=255)
     tAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCCgcagca                                                                                                                                                 >  4:56398/1‑150 (MQ=255)
     tAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCCgcagca                                                                                                                                                 >  4:170219/1‑150 (MQ=255)
       gCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATc                                                                                                                                               <  4:38449/150‑1 (MQ=255)
         aaCGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCgg                                                                                                                                             >  1:574207/1‑150 (MQ=255)
          aCGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCggg                                                                                                                                            <  3:245433/150‑1 (MQ=255)
           cGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGa                                                                                                                                           <  4:272047/150‑1 (MQ=255)
            gATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGAt                                                                                                                                          >  4:369014/1‑150 (MQ=255)
                  aaTCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCaa                                                                                                                                    >  1:150814/1‑150 (MQ=255)
                    tCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTCCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAAAGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAAt                                                                                                                                  <  4:194379/150‑1 (MQ=255)
                       ccaaaccaaCAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACg                                                                                                                               >  2:237345/1‑150 (MQ=255)
                       ccaaaccaaCAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACg                                                                                                                               >  3:92092/1‑150 (MQ=255)
                        caaaccaaCAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGc                                                                                                                              >  1:523825/1‑150 (MQ=255)
                        caaaccaaCAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGc                                                                                                                              >  3:330901/1‑150 (MQ=255)
                         aaaccaaCAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGcc                                                                                                                             <  2:180661/150‑1 (MQ=255)
                          aaccaaCAAGCACCGGTACCGACGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCg                                                                                                                            <  3:320836/150‑1 (MQ=255)
                                 aaGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGGCAAATACGCCGGGCTaag                                                                                                                     >  2:408509/1‑149 (MQ=255)
                                  aGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATg                                                                                                                    <  1:495892/150‑1 (MQ=255)
                                   gCACCGGTACCGACGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGc                                                                                                                   <  2:266273/150‑1 (MQ=255)
                                                  tGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCa                                                                                                    <  2:515602/150‑1 (MQ=255)
                                                  tGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTAACGCCTGCAGCATCGGGATAGTCAAATACGCCGGCTTAATGCCGGAAGAGTGCGCCa                                                                                                    >  2:380160/1‑150 (MQ=255)
                                                     gTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGtt                                                                                                 >  1:322572/1‑150 (MQ=255)
                                                     gTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGtt                                                                                                 >  1:45247/1‑150 (MQ=255)
                                                         tACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCgg                                                                                             >  2:438673/1‑150 (MQ=255)
                                                          aCGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCggg                                                                                            <  1:14324/150‑1 (MQ=255)
                                                            gACTTCAAACGGCGAGATGTTAGCAATCCCGGCAAGCGAAACGAATCCGCCAGAACCGGGCGAACGGGTACGACCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAGATACGCCGGGTTACTGCCGTAAGAGTGCGCCAGTTTCGGGat                                                                                          <  4:188903/150‑1 (MQ=255)
                                                            gACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGCTTCGGGat                                                                                          >  3:71778/1‑150 (MQ=255)
                                                             aCTTCAAACGGCGAGATCTCCGCCATCCCGTCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGatc                                                                                         <  1:54941/150‑1 (MQ=255)
                                                              cTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGatca                                                                                        <  4:193340/150‑1 (MQ=255)
                                                               ttCAAACGGCGAGATATTCGCCATTCCGGCAACCGCAACAACTACACCAGAAAAGAGAGAAAGGGTACGGCCAAGGTTTGACGCAGGCAACATCGGGATATTCAAATACGCTAGGTTCATGCCGGAAAACTTCGAAAGTTTAGGGatcat                                                                                       >  2:446702/1‑150 (MQ=255)
                                                                 cAAACGGCGGGAACTTAGACATCGTGGCACCCGCAATCACTACGCCAGAAACGGGCGAAAGAGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCACATACGCCTGGTATATCCCGGAATTGTGCGCCAGTTTCGCGATCAtct                                                                                     <  4:189119/150‑1 (MQ=255)
                                                                  aaaCGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCAtctc                                                                                    >  3:48594/1‑150 (MQ=255)
                                                                     cGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTcaac                                                                                 >  1:325253/1‑150 (MQ=255)
                                                                     cGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTcaac                                                                                 <  3:30910/150‑1 (MQ=255)
                                                                        cGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACaaa                                                                              <  2:569298/150‑1 (MQ=255)
                                                                         gagaTCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAAc                                                                             <  2:402525/150‑1 (MQ=255)
                                                                                 cGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATaaaa                                                                     <  3:305341/150‑1 (MQ=255)
                                                                                   ccATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACg                                                                   >  2:222061/1‑150 (MQ=255)
                                                                                        ccGGCAACCGCAACGACTACGCCAGAAACGGGCGCAAGGGTACGACCAAGGTTTGACGCCGGCAGCATCAGGATAGTGAAGTGCGCCGGGGTAATGCCGGAAGAGGGCCCCAGTTTCGGGATCATCACAACAAACGCATAAAACGGCGCa                                                              >  2:75353/1‑150 (MQ=255)
                                                                                          ggCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCAtt                                                            <  4:158444/150‑1 (MQ=255)
                                                                                                         tACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGCGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGAATTGCCTGACCCGGTCGt                                             >  2:344693/1‑150 (MQ=255)
                                                                                                           cGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGa                                           <  1:326463/150‑1 (MQ=255)
                                                                                                           cGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGa                                           <  2:400514/150‑1 (MQ=255)
                                                                                                           cGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGa                                           <  3:207626/150‑1 (MQ=255)
                                                                                                           cGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCCGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGa                                           <  3:56398/150‑1 (MQ=255)
                                                                                                           cGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCCGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGa                                           <  3:170219/150‑1 (MQ=255)
                                                                                                                gAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTgcc                                      <  4:148553/150‑1 (MQ=255)
                                                                                                                     gggCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAACTACGCCGGGTTACTGCCGGAAGAGTGCGCCAGTTTCGGGCTCATCTCAACAAACGCATACAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAg                                 >  2:138496/1‑150 (MQ=255)
                                                                                                                      ggCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGc                                >  2:67882/1‑150 (MQ=255)
                                                                                                                            aGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTc                          >  4:241343/1‑150 (MQ=255)
                                                                                                                                 aCGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAAt                     >  3:379433/1‑150 (MQ=255)
                                                                                                                                   ggCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCa                   <  2:374927/150‑1 (MQ=255)
                                                                                                                                     ccAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATCCCGGAAGAGTGCGCCAGCTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCAcc                 <  4:330901/150‑1 (MQ=255)
                                                                                                                                      cAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCa                <  2:386817/150‑1 (MQ=255)
                                                                                                                                          gTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATAc            <  2:322572/150‑1 (MQ=255)
                                                                                                                                          gTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATAc            <  1:438673/150‑1 (MQ=255)
                                                                                                                                            ttGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAATCGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCa          <  2:88526/150‑1 (MQ=255)
                                                                                                                                              gACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGc        >  1:394515/1‑150 (MQ=255)
                                                                                                                                               aCGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCa       >  1:313357/1‑150 (MQ=255)
                                                                                                                                                    tGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAgg  <  1:580088/150‑1 (MQ=255)
                                                                                                                                                    tGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAgg  <  1:268294/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCATCCCGGCAACCGCAACGACTACGCCAGAAACGGGCGAAAGGGTACGGCCAAGGTTTGACGCCTGCAGCATCGGGATAGTCAAATACGCCGGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGG  >  NC_000913_3_hsa_tpiA/654625‑654922

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