Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 4,001,994 T→C 6.1% L190P (CTG→CCG)  corA → magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,001,9940TC 6.1% 209.0 / 3.8 66L190P (CTG→CCG) corAmagnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base C (2/2);  total (29/37)
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.54e-01

GAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGAAGGTGAACTTCC  >  NC_000913_3_hsa_tpiA/4001849‑4002143
                                                                                                                                                 |                                                                                                                                                     
gAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGCCCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTgtgt                                                                                                                                                   <  3:326892/150‑1 (MQ=255)
gAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTgtgt                                                                                                                                                   <  4:181988/150‑1 (MQ=255)
 aaCAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTc                                                                                                                                                  >  2:363083/1‑150 (MQ=255)
 aaCAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCCGTGTc                                                                                                                                                  >  4:60946/1‑150 (MQ=255)
 aaCAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCCGTGTc                                                                                                                                                  >  3:118556/1‑150 (MQ=255)
    aGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGa                                                                                                                                               <  1:90029/150‑1 (MQ=255)
    aGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGa                                                                                                                                               <  2:123025/150‑1 (MQ=255)
    aGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGa                                                                                                                                               <  1:253720/150‑1 (MQ=255)
     gTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGAt                                                                                                                                              <  1:484314/150‑1 (MQ=255)
     gTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGAt                                                                                                                                              <  1:98171/150‑1 (MQ=255)
      ttGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATg                                                                                                                                             <  2:137876/150‑1 (MQ=255)
       tGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATgg                                                                                                                                            <  2:117188/150‑1 (MQ=255)
                aaaTTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAgc                                                                                                                                   >  2:501351/1‑150 (MQ=255)
                  aTTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAgcgc                                                                                                                                 >  3:383253/1‑150 (MQ=255)
                  aTTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAgcgc                                                                                                                                 <  3:225386/150‑1 (MQ=255)
                        aaTATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGTGCTc                                                                                                                           <  2:210882/150‑1 (MQ=255)
                             cTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACtt                                                                                                                      >  3:51587/1‑150 (MQ=255)
                                   cGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGt                                                                                                                >  2:319638/1‑150 (MQ=255)
                                   cGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGt                                                                                                                >  1:534614/1‑150 (MQ=255)
                                       cTGGGGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCGCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCCCAGCTTCCTGGTGCGt                                                                                                            <  2:527187/150‑1 (MQ=255)
                                         ggAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCTCGCTCAACTTCCTGATGCGTaa                                                                                                          >  4:168316/1‑150 (MQ=255)
                                         ggAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACACAGCGCGCGCTCAACTTCCTGGTGCGTaa                                                                                                          >  4:276006/1‑150 (MQ=255)
                                          gAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTaaa                                                                                                         >  3:237769/1‑150 (MQ=255)
                                                ttGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGt                                                                                                   >  4:249671/1‑150 (MQ=255)
                                                    gCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTAc                                                                                               <  4:223680/150‑1 (MQ=255)
                                                    gCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTAc                                                                                               <  2:213936/150‑1 (MQ=255)
                                                     ccGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTAcc                                                                                              >  1:484106/1‑150 (MQ=255)
                                                             ttATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCGGGTGCGTAAAGCGCGTTTACCGGGTGGGc                                                                                      >  2:367031/1‑150 (MQ=255)
                                                                          tCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGgc                                                                         >  3:184506/1‑150 (MQ=255)
                                                                            aGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGgcgc                                                                       >  4:103473/1‑150 (MQ=255)
                                                                            aGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGgcgc                                                                       >  2:161994/1‑150 (MQ=255)
                                                                              ggCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGt                                                                     >  2:554354/1‑150 (MQ=255)
                                                                                cGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGa                                                                   <  4:260449/150‑1 (MQ=255)
                                                                                    gAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATc                                                               <  3:374303/150‑1 (MQ=255)
                                                                                    gAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTAGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATc                                                               <  2:125716/150‑1 (MQ=255)
                                                                                       tACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTg                                                            >  2:462536/1‑150 (MQ=255)
                                                                                        acgacgAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGc                                                           <  1:309625/150‑1 (MQ=255)
                                                                                        acgacgAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTCCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGc                                                           <  1:255298/150‑1 (MQ=255)
                                                                                         cgacgaGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCg                                                          >  4:365898/1‑150 (MQ=255)
                                                                                            ggtcgcgcTCTGCACTCAGGCGGAACTGGAAGAGAGCGGCTGGAAAGTTCGCCTGTGCCCGAGGGATACCCGGCGCGCGCGCAACTTCGTGGTGCGTAAAGCGGGTTGACCGGGTGGGCAGCCGGTTCAGGCGCGTGAAATCCTGCgagt                                                       <  4:338197/146‑2 (MQ=255)
                                                                                                 cgcTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCg                                                  <  4:121450/150‑1 (MQ=255)
                                                                                                  gcTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGa                                                 >  3:202569/1‑150 (MQ=255)
                                                                                                  gcTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGa                                                 >  4:25046/1‑150 (MQ=255)
                                                                                                  gcTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGGCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTCTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGa                                                 >  2:270806/1‑150 (MQ=255)
                                                                                                   ctctCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGaa                                                <  2:215235/150‑1 (MQ=255)
                                                                                                     ctcCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATc                                              <  2:487972/150‑1 (MQ=255)
                                                                                                     ctcCAACCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTAAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGTAGGCGCGTGAAATCCTGCGAGATATCGAATc                                              <  2:311533/150‑1 (MQ=255)
                                                                                                      tcCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAGTCCTGCGAGATATCGAATcc                                             <  2:357197/150‑1 (MQ=255)
                                                                                                      tcCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATcc                                             <  1:363083/150‑1 (MQ=255)
                                                                                                        cACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCCGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCct                                           <  3:60946/150‑1 (MQ=255)
                                                                                                        cACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCCGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCct                                           <  4:118556/150‑1 (MQ=255)
                                                                                                                  gAACTNGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCANGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCAt                                 <  1:51917/150‑1 (MQ=255)
                                                                                                                    aCTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATaa                               <  1:365061/150‑1 (MQ=255)
                                                                                                                    aCTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATaa                               <  1:177804/150‑1 (MQ=255)
                                                                                                                    aCTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATaa                               <  1:464911/150‑1 (MQ=255)
                                                                                                                     cTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAAc                              >  4:213374/1‑150 (MQ=255)
                                                                                                                        gAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGaa                           >  1:219907/1‑150 (MQ=255)
                                                                                                                               tCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGt                    <  2:194278/150‑1 (MQ=255)
                                                                                                                                  gCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTcc                 <  1:308571/150‑1 (MQ=255)
                                                                                                                                   cTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGTTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCa                <  1:85599/150‑1 (MQ=255)
                                                                                                                                     ggAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGa              <  3:223794/150‑1 (MQ=255)
                                                                                                                                      gAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGaa             <  3:6026/150‑1 (MQ=255)
                                                                                                                                       aaaGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGAAg            >  1:62996/1‑150 (MQ=255)
                                                                                                                                           ttCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGAAGGTGa        >  1:110062/1‑150 (MQ=255)
                                                                                                                                           ttCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGAAGGTGa        >  1:110077/1‑150 (MQ=255)
                                                                                                                                               ccTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGAAGGTGAACtt    >  4:41949/1‑150 (MQ=255)
                                                                                                                                                 tgtgTCTGATGGATGCCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGAAGGTGAACTTcc  <  4:216347/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCGATGAGTACGACGAGGCGCTCTCCACTCTGGCGGAACTGGAAGATATCGGCTGGAAAGTTCGCCTGTGTCTGATGGATACCCAGCGCGCGCTCAACTTCCTGGTGCGTAAAGCGCGTTTACCGGGTGGGCAACTGGAGCAGGCGCGTGAAATCCTGCGAGATATCGAATCCCTGCTGCCGCATAACGAATCCCTGTTCCAGAAGGTGAACTTCC  >  NC_000913_3_hsa_tpiA/4001849‑4002143

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