Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_pgi 3,815,809 Δ1 bp 19.4% intergenic (‑41/+25) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi3,815,8090G.19.4% 144.3 / 46.4 67intergenic (‑41/+25)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base . (6/7);  total (31/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 = 2.41e-01

GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCC  >  NC_000913_3_hsa_pgi/3815664‑3815946
                                                                                                                                                 |                                                                                                                                          
gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCg                                                                                                                                         >  2:21873/1‑149 (MQ=255)
 aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTACTCTTACTTTTCTCCAGCCAAAAAAAAGGCGa                                                                                                                                        >  2:487562/1‑149 (MQ=255)
 aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTAAGGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAACAGGGGa                                                                                                                                        >  2:914135/1‑149 (MQ=255)
   gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTCAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCTTCTTACTTTTCTCCATCACAAAAAAAGTCGACt                                                                                                                                      >  1:167991/1‑149 (MQ=255)
     aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtca                                                                                                                                    <  1:609190/149‑1 (MQ=255)
        aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatca                                                                                                                                 >  1:836577/1‑149 (MQ=255)
              atttGCGCCCGTGTTTCAGCGAAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTGTAACGTCACAAAAAAGGCGAATCATCAGTcccc                                                                                                                           >  2:373612/2‑149 (MQ=255)
                     cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa                                                                                                                    >  2:83698/1‑149 (MQ=255)
                      ccGGATTTCGGCGGATGCTCGCGAGACTTTAAAACCTGCTTGCTAAGCGCAAATGCAATAAACTGGCGCTGATGTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGCCTAATCCGTCCCACTAAAAAt                                                                                                                   >  2:241607/1‑149 (MQ=255)
                      ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                   >  2:474026/1‑149 (MQ=255)
                       cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCAACTCATCAGTCGCCTTAAAAATc                                                                                                                  >  1:283264/1‑149 (MQ=255)
                            ttCTGCGTAAACTCTCCAAATTTTAACATCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                                                             <  1:980504/149‑1 (MQ=255)
                               aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑cc                                                                                                         <  1:21873/149‑1 (MQ=255)
                               aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑cc                                                                                                         >  2:788396/1‑149 (MQ=255)
                                            cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCT                                                                                            >  1:600336‑M1/1‑149 (MQ=255)
                                             aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCgc                                                                                           >  2:224867‑M1/1‑147 (MQ=255)
                                             aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTG                                                                                           >  2:978352‑M1/1‑149 (MQ=255)
                                             aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGC                                                                                          >  1:971026‑M1/1‑149 (MQ=255)
                                                    aaaacCTGCTTGCTAAGCCCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTa                                                                                   <  2:288565/146‑1 (MQ=255)
                                                    aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCGGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACTAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTa                                                                                   <  1:777245/149‑1 (MQ=255)
                                                       acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCCCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctg                                                                                 >  1:362882‑M1/1‑140 (MQ=255)
                                                        ccTGTTTGCTAAGCGCAAAGTCTATAGACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTATAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CTAGCGTCGCCTTCTGCGTCGCTACaa                                                                                <  2:820160/149‑1 (MQ=255)
                                                        acTGCTCGCAAAGCGCAAATTCAACAAACTGGCGCTGATACTGTTTCATTCCTTCGCCCCTCATCTTACTTTTCTACATACCAAAAAAAGGCGACTCATCAGTCGCATTAAAAATCAGTTCC‑CCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                <  2:491646/148‑1 (MQ=255)
                                                         cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                               >  2:580595/1‑149 (MQ=255)
                                                           gcttgctAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATgg                                                                             <  1:644511/149‑1 (MQ=255)
                                                             tgtgtaagcCCAACTTCAGTAAACTAGGGCTTATATTGTGTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATTAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttc                                                                           <  2:1001009‑M1/145‑16 (MQ=255)
                                                              tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                          <  1:580595/149‑1 (MQ=255)
                                                               gctAAGCGCAAATTAAATAAACTGGGGCTGATATGGTTTCATGGCTTCGCGCCTCATCTTACTTTTCTAAAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                         <  2:18125/149‑1 (MQ=255)
                                                                ctAAGCTCAAATTCAATTGACTGGCGCTGATATGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                        <  1:621924‑M1/149‑19 (MQ=255)
                                                                 tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGATAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                       <  1:755976/149‑1 (MQ=255)
                                                                 tAAGCGCAAATTCAATAAACTGGCGCTGACATGGTTTCATGCCTTCGCACCTCATCGTACTTTTCTACAGACAAAAAAAA‑GCGACCCATCAGTCGCCTTAAAAATCAGCTTG‑CCAGCGCCGCCCTCTGCGTCGCTACCATGgatccgat                                                                      >  1:202154/1‑149 (MQ=255)
                                                                     cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAAt                                                                                                                   >  2:513104/1‑101 (MQ=255)
                                                                     cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAAt                                                                                                                   <  1:513104/101‑1 (MQ=255)
                                                                       cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                        >  1:84205‑M1/1‑124 (MQ=255)
                                                                       cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                        <  2:84205‑M1/142‑19 (MQ=255)
                                                                        aaaTTCAATAAACTGGCGCTGATATGGTTCCAAACCCTGGGTCCTCATCTTACTCTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACACTGGCTTCGATTCCCCTc                                                               >  2:468143/1‑149 (MQ=255)
                                                                           ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTCTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                             <  1:984358/149‑1 (MQ=255)
                                                                            tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                                                            >  1:148371‑M1/1‑119 (MQ=255)
                                                                             caactaccTGTCGCTTACATGGTTGCATGCCTTGGCTCCTCATCTTTCTTTTCTACAGACAAAAAAAAGGCGACTCATCCGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTTCGTCGCTACATTGGATTCGATTCCCCTCCGGc                                                           <  2:407317/142‑1 (MQ=255)
                                                                                aaacatggcgcTGATCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                        <  1:226491‑M1/144‑35 (MQ=255)
                                                                                  aaCTGGCGCTGTTATGGTTTCATGCCTTCGCTCCTCATCTTATTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCagag                                                      <  1:347032/149‑1 (MQ=255)
                                                                                      ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaa                                                  <  2:635208/149‑1 (MQ=255)
                                                                                         tcTGATATCGTGTCATGGCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                               <  2:316072/148‑1 (MQ=255)
                                                                                         gcTGATATGGTTACATAACTTAGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctggattctaattttccgcccggcttcggtcatcact                                               >  2:460204‑M1/1‑106 (MQ=255)
                                                                                          gtgtgatGGTTTTACTCCTCCGCTGCTCATCTTACATTTCTATAGACAAAAGAAA‑GCGACTCATCAGTCGCCTCAAAAATCCGTTTG‑CCAGCGCCGCCTTCTGCGTGGCTACAATGTATTCGATTCCCCTCGGGCCAGAGCAAACAAGa                                             <  2:952585/144‑1 (MQ=255)
                                                                                           tGATATGGTTTTATGCCTTCGCTCCTCATCTTACTTTTCTACAGAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                             <  2:354708/149‑1 (MQ=255)
                                                                                            gATATTGTGTTATGACTTCGCTCCTCAACTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt                                            <  2:271547/149‑1 (MQ=255)
                                                                                             ttatTGTTTCCAGCCTTCGCTGCTTATCTTAAGTTTCTACAGACATAAAAAATGCGATTTATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCTATTTCCCTCGGGCCAGAGCCAACAAGATg                                           <  2:630379/148‑1 (MQ=255)
                                                                                                        tGCCTTCGCTCTTCATCTTTCTTTTCTACTGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑ACAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                <  1:260017/149‑1 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                        <  2:297227/137‑1 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                        >  1:297227/1‑137 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGCGGCAGAGCCAACAAGATGAGTAGCTCTTCATgg                            <  2:971026/149‑1 (MQ=255)
                                                                                                                ctcctcATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                       <  2:202154/149‑1 (MQ=255)
                                                                                                                  cgtcATCTTACTTTTCTACAGAAAAAAAAAAGGCGACTCATCAGCTGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTTTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                      <  1:136763/147‑1 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                     <  2:148371‑M1/149‑70 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATTGATTCGATTCCCCTCGGGCCAGAGCCACCAAGCTGATTAGCTCTTCATGGGTGCAcc                     >  2:236526/1‑148 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGAAAAAATAAAGGCGACTCATCAGTCGCCTT‑AAAATCAGGCCGTGAAGGGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                     <  1:484157‑M1/149‑70 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                        >  2:539692‑M1/1‑70 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                        <  1:539692‑M1/104‑35 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCACCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCAGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc           >  1:282988/1‑149 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CAAGCGCCGCCTTATGAcgaccactaaacttcaaggaatcccccgtcttcggtcatcactacgttcatgtcggtctctgcggcaaagtcttcaacg           >  2:362806‑M1/1‑70 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCGCAAAAATCCGTTTG‑CCAGCGCCACCTACTGCGTAGCTACAATGGATTCGATTACCCTCGGGCCAGAGCCAACAAGATGAGTACCTCTTCACCGGTCACCGCCCCGCCtac           >  2:223097/1‑147 (MQ=255)
                                                                                                                               atCTACAGAGAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCCCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg         <  2:927041/148‑1 (MQ=255)
                                                                                                                                  tACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAGTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGGTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt     >  2:518848/1‑149 (MQ=255)
                                                                                                                                   aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt     <  2:525043/149‑1 (MQ=255)
                                                                                                                                     aGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAACTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTccc  >  1:197271/1‑149 (MQ=255)
                                                                                                                                      gACAAAAAAAAGGCGACTCATCAGCCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc  <  2:1000636/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTG‑CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCC  >  NC_000913_3_hsa_pgi/3815664‑3815946

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