Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 3,909,054 A→G 5.6% L129P (CTC→CCC)  pstA ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA3,909,0540AG 5.6% 218.3 / 3.1 73L129P (CTC→CCC) pstAphosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base G (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

CCAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGC  >  NC_000913_3_hsa_tpiA/3908905‑3909202
                                                                                                                                                     |                                                                                                                                                    
ccAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagaga                                                                                                                                                      <  3:9229/150‑1 (MQ=255)
 cAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAGCAGACCAACCACAATCGACGGCGCagagag                                                                                                                                                     >  1:375962/1‑150 (MQ=255)
 cAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagagag                                                                                                                                                     >  1:15587/1‑150 (MQ=255)
 cAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCagagag                                                                                                                                                     >  3:238054/1‑150 (MQ=255)
   gTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCa                                                                                                                                                   <  1:129170/150‑1 (MQ=255)
       cAGCATGTTCTCGGTGGTGCGGGTAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAAt                                                                                                                                               >  2:293782/1‑150 (MQ=255)
          cATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTc                                                                                                                                            <  1:280192/150‑1 (MQ=255)
           aTGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCg                                                                                                                                           <  3:250899/150‑1 (MQ=255)
           aTGTTCGCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCCCAATCGACGGCGCAGAGAGCAGAATGTCg                                                                                                                                           <  1:212478/150‑1 (MQ=255)
            tGTTCTCGGTGGTGCGGATAACAATCAGCACCTGCAACAACGCCAGGGCAATCACGCTCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATTGACGGTGCAGAGTGCAGATTGTCGt                                                                                                                                          >  2:259872/1‑150 (MQ=255)
               tgccGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATAACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTaa                                                                                                                                       <  2:471523/147‑1 (MQ=255)
                ctcGGTGGTCGGCGTAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGACCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAGACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATCACGTTaag                                                                                                                                      >  4:68955/1‑149 (MQ=255)
                  cGGTGGTGCGGATATCAATCGGCACCTGCAACAACGCCAGGGCGATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCAACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGa                                                                                                                                    <  2:308473/150‑1 (MQ=255)
                        tgCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGaa                                                                                                                              >  1:45364/1‑150 (MQ=255)
                         gCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAAt                                                                                                                             >  1:399600/1‑150 (MQ=255)
                           ggATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCa                                                                                                                           >  1:548922/1‑150 (MQ=255)
                                aCAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCt                                                                                                                      <  4:326126/150‑1 (MQ=255)
                                 cAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGACTGTCGTTAATGAAGCGAATCACTTCTg                                                                                                                     >  1:430479/1‑150 (MQ=255)
                                 cAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTg                                                                                                                     >  2:313750/1‑150 (MQ=255)
                                      ggCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgccagc                                                                                                                <  3:135695/150‑1 (MQ=255)
                                        cagcTGCAACAACGACAGGGCAATCACGCCCGCCCAGCCGGAGACGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACATCTgccagcca                                                                                                              <  2:351332/147‑1 (MQ=255)
                                        cACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTgccagcca                                                                                                              >  1:325158/1‑150 (MQ=255)
                                         aCCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCCCCACAACGGCGTCAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGccagccag                                                                                                             <  4:91996/150‑1 (MQ=255)
                                          ccTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAgg                                                                                                            >  1:439670/1‑150 (MQ=255)
                                          ccTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTCTAAACAAACAGACCAACCACAATCGACCGCGCAGAGAGCACAATGTCGTCAATGAAGCGAATCACCTCTGCCAGCCAgg                                                                                                            >  4:349972/1‑150 (MQ=255)
                                              caacaaCGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGAttt                                                                                                        >  3:118181/1‑150 (MQ=255)
                                                acaacGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTAc                                                                                                      <  4:276802/150‑1 (MQ=255)
                                                         gggCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATATGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATAtt                                                                                             >  4:147871/1‑150 (MQ=255)
                                                              aTCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGcc                                                                                        >  1:221299/1‑150 (MQ=255)
                                                              aTCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGcc                                                                                        >  1:519413/1‑150 (MQ=255)
                                                                 aCGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGa                                                                                     >  4:348204/1‑150 (MQ=255)
                                                                         ccAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATccc                                                                             >  2:500504/1‑150 (MQ=255)
                                                                           aGCCGGAGAAGTGCTCCATCTGCGCCACCACAGTGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCg                                                                           <  1:231860/150‑1 (MQ=255)
                                                                             ccGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGcc                                                                         >  1:385554/1‑150 (MQ=255)
                                                                               ggAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCAt                                                                       <  4:3881/150‑1 (MQ=255)
                                                                                   aaGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGGGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATg                                                                   >  3:148331/1‑150 (MQ=255)
                                                                                   aaGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGGGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATg                                                                   >  3:148330/1‑150 (MQ=255)
                                                                                   aaGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATg                                                                   >  4:183384/1‑150 (MQ=255)
                                                                                     gTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCCGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGcc                                                                 >  4:279189/1‑150 (MQ=255)
                                                                                      tGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGccc                                                                <  1:1192/150‑1 (MQ=255)
                                                                                      tGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGccc                                                                >  2:352568/1‑150 (MQ=255)
                                                                                      tGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAACCGACGGCGCAGAGAGCAGAATGGCGTTAACGAAGCGAATCACTTCTGCCACCCAGGATTTACGACCATATTCCCCCAGATAAATCCCCGCCATAATGccc                                                                <  4:55835/150‑1 (MQ=255)
                                                                                        cTCCATCTGCGCCACCACAATGGTGTACACAAACAGCCCAACCACAATCGACGGCGCAGAGAGCAGAATGCCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAg                                                              >  4:296624/1‑150 (MQ=255)
                                                                                        cTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAg                                                              <  4:301846/150‑1 (MQ=255)
                                                                                        cTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAg                                                              <  3:315702/150‑1 (MQ=255)
                                                                                             tCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCCGACTTTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAgcggcg                                                         >  1:276454/1‑150 (MQ=255)
                                                                                                 cgcCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTAcc                                                     <  2:10024/150‑1 (MQ=255)
                                                                                                    caccacAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGaa                                                  <  3:284470/150‑1 (MQ=255)
                                                                                                      ccacAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATa                                                >  1:283142/1‑150 (MQ=255)
                                                                                                        acaATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATAcc                                              >  2:162261/1‑150 (MQ=255)
                                                                                                        acaATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATAcc                                              >  4:384618/1‑150 (MQ=255)
                                                                                                         caATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGCCCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCg                                             <  2:89143/150‑1 (MQ=255)
                                                                                                            tGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCCTATTCCGCCAGATAAAGCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTgg                                          <  2:198510/150‑1 (MQ=255)
                                                                                                            tGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTgg                                          <  3:349972/150‑1 (MQ=255)
                                                                                                             ggTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAAGGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATCTTCCGCCAGATAAAGCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGc                                         <  2:198106/150‑1 (MQ=255)
                                                                                                               tgtAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGTAGAATGTCGTTAATGAAGCGAATCACTTGTGCCAGCCAGGACTTACGACCTTATTCCTCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGGATACCGTGGccc                                       <  2:15587/150‑1 (MQ=255)
                                                                                                               tgtAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGccc                                       <  4:238054/150‑1 (MQ=255)
                                                                                                                  aaacaaacAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCcaca                                    <  2:439670/150‑1 (MQ=255)
                                                                                                                   aacaaacaGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACaa                                   >  2:51578/1‑150 (MQ=255)
                                                                                                                     caaacaGACCAACCACAATCGACGGCGCAGAGGGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACaaaa                                 <  4:148330/150‑1 (MQ=255)
                                                                                                                     caaacaGACCAACCACAATCGACGGCGCAGAGGGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACaaaa                                 <  4:148331/150‑1 (MQ=255)
                                                                                                                      aaacaGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAAt                                >  4:331841/1‑150 (MQ=255)
                                                                                                                         caGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTaa                             <  3:257558/150‑1 (MQ=255)
                                                                                                                         caGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTaa                             <  1:193211/150‑1 (MQ=255)
                                                                                                                         caGACCAACCACAATCGACGGAACAGAGAGCAGAATGTCGTAAATGAAGAGAATCACTTCCGCAAGCTAGGATTTACGACCAGATTCCGCCAGACAATTCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAACATTaa                             <  4:23034/150‑1 (MQ=255)
                                                                                                                          aGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAAc                            <  3:348204/150‑1 (MQ=255)
                                                                                                                            accaaccaCAATCGACGGCGCAGAGAGCAGACTGTCGTTAATGCAACGAATCACTTCTGCCAGCCGGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACac                          >  4:67007/1‑149 (MQ=255)
                                                                                                                            accaaccaCAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAg                          >  2:465269/1‑150 (MQ=255)
                                                                                                                             ccaaccaCAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGc                         <  4:118181/150‑1 (MQ=255)
                                                                                                                                  cacaATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCt                    >  1:505015/1‑150 (MQ=255)
                                                                                                                                  cacaATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCt                    >  3:410/1‑150 (MQ=255)
                                                                                                                                       tCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCg               >  3:209399/1‑150 (MQ=255)
                                                                                                                                          aCGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCACGCTTCACGACCCTATTCAGCCAGATAAATCCTCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCa            <  1:259872/150‑1 (MQ=255)
                                                                                                                                           cGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCag           <  1:959/150‑1 (MQ=255)
                                                                                                                                             gcgcAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCagag         >  2:320686/1‑150 (MQ=255)
                                                                                                                                              cgcAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGc        >  2:489370/1‑150 (MQ=255)
                                                                                                                                              cgcAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGc        <  3:147871/150‑1 (MQ=255)
                                                                                                                                                   agagCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCg   >  1:548137/1‑150 (MQ=255)
                                                                                                                                                   agagCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCg   >  3:209727/1‑150 (MQ=255)
                                                                                                                                                    gagCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGc  <  1:207834/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCACCTGCAACAACGCCAGGGCAATCACGCCCGCCCAGCCGGAGAAGTGCTCCATCTGCGCCACCACAATGGTGTAAACAAACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGC  >  NC_000913_3_hsa_tpiA/3908905‑3909202

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: