Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,538,717 A→G I246I (ATT→ATC)  unknown ← Porphobilinogen deaminase (EC 2.5.1.61)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,538,7170AG100.0% 104.4 / NA 39I246I (ATT→ATC) unknownPorphobilinogen deaminase (EC 2.5.1.61)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/1);  total (38/1)

GCTGACCGCGTTCATCGGGGCGGGTGACCAAAACTGCCATACGTTTTCTTCCTTACGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCAATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAGCCAAAATAATGGCAT  >  1/1538476‑1538961
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gctgACCGCGTTCATCGGGGCGGGTGACCAAAACTGCCATACGTTTTCTTCCTTACGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCACCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCa                                                                                                                                                                                                                                               >  1:357159/1‑249 (MQ=255)
       gcgTTCATCGGGGCGGGTGACCAAAACTGCCATACGTTTTCTTCCTTACGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAg                                                                                                                                                                                                                                       >  1:316613/1‑250 (MQ=255)
            cATCGGGGCGGGTGACCAAAACTGCCATACGTTTTCTTCCTTACGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAgccgc                                                                                                                                                                                                                                   >  1:526114/1‑249 (MQ=255)
                         gACCAAAACTGCCATACGTTTTCTTCCTTACGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGt                                                                                                                                                                                                                    >  1:656226/1‑251 (MQ=255)
                                   gCCATACGTTTTCTTCCTTACGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCAc                                                                                                                                                                                                                                              >  1:14638/1‑215 (MQ=255)
                                                      aCGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGAt                                                                                                                                                                                                                       >  1:878366/1‑219 (MQ=255)
                                                       cgcgTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAgg                                                                                                                                                                                      >  1:662889/1‑251 (MQ=255)
                                                         cgTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTa                                                                                                                                                                                    >  1:450939/1‑251 (MQ=255)
                                                           tAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGt                                                                                                                                                                                  >  1:503880/1‑251 (MQ=255)
                                                           tAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGt                                                                                                                                                                                  >  1:142376/1‑251 (MQ=255)
                                                                          aTTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAg                                                                                                                                                                   >  1:868728/1‑251 (MQ=255)
                                                                          aTTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAg                                                                                                                                                                   >  1:331717/1‑251 (MQ=255)
                                                                              tATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTgg                                                                                                                                                                >  1:159871/1‑250 (MQ=255)
                                                                                    cccccTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGt                                                                                                                                                         >  1:848679/1‑251 (MQ=255)
                                                                                        cTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGcc                                                                                                                                                                     >  1:606805/1‑235 (MQ=255)
                                                                                            gcCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAgc                                                                                                                                                                                                                                      >  1:219865/1‑166 (MQ=255)
                                                                                            gcCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCg                                                                                                                                                                                                                         >  1:437068/1‑179 (MQ=255)
                                                                                                             aaTTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGt                                                                                                                                   >  1:731696/1‑248 (MQ=255)
                                                                                                                       ccTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACAtt                                                                                                                       >  1:255939/1‑250 (MQ=255)
                                                                                                                       ccTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACAtt                                                                                                                       >  1:881914/1‑250 (MQ=255)
                                                                                                                       ccTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACAtt                                                                                                                       >  1:602825/1‑250 (MQ=255)
                                                                                                                                                                       atCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACAt                                                                       >  1:635468/1‑250 (MQ=255)
                                                                                                                                                                         cgtcgAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCaa                                                                    >  1:591488/1‑251 (MQ=255)
                                                                                                                                                                         cgtcgAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCaa                                                                    >  1:193290/1‑251 (MQ=255)
                                                                                                                                                                                           gCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACg                                                    >  1:587077/1‑249 (MQ=255)
                                                                                                                                                                                             ccACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGAt                                                                                 >  1:251800/1‑218 (MQ=255)
                                                                                                                                                                                             ccACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGct                                                  >  1:211947/1‑249 (MQ=255)
                                                                                                                                                                                             ccACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGct                                                  >  1:482973/1‑249 (MQ=255)
                                                                                                                                                                                                        cgTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACAt                                                                                                                                      >  1:32183/1‑154 (MQ=255)
                                                                                                                                                                                                        cgTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTaa                                     >  1:648624/1‑251 (MQ=255)
                                                                                                                                                                                                        cgTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGATCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTaa                                     >  1:579286/1‑251 (MQ=255)
                                                                                                                                                                                                                cAATTCGCCATCGCGCAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCt                                       >  1:183870/1‑241 (MQ=255)
                                                                                                                                                                                                                            gcgcAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACAt                                                                       >  1:890377/1‑197 (MQ=255)
                                                                                                                                                                                                                            gcgcAACACCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAg                 >  1:95910/1‑251 (MQ=255)
                                                                                                                                                                                                                              cccaccccaaacacccgccgATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAGcc               <  1:816849/237‑1 (MQ=255)
                                                                                                                                                                                                                                 acacCGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAGCCaaa            >  1:740031/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ccgaaaacaCCGCCGATAGGCACTTGGAAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAGGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAGCCAAaat          <  1:608984/243‑1 (MQ=255)
                                                                                                                                                                                                                                    ccGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATa                                                      >  1:42926/1‑206 (MQ=255)
                                                                                                                                                                                                                                    ccGCATAACCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAAcc                                        >  1:306777/1‑220 (MQ=255)
                                                                                                                                                                                                                                      caacacccgccgATGGGCACTTGGCAGCGGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAGCCAAaataat       <  1:388007/245‑1 (MQ=255)
                                                                                                                                                                                                                                           aCCGCCGATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAGCCAAAATAATGGcat  >  1:598505/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GCTGACCGCGTTCATCGGGGCGGGTGACCAAAACTGCCATACGTTTTCTTCCTTACGCGTAAATATCCTGCAAAATTTTATCCGCCCCCTGCGCCAATAGCTGTTCGGCAATTTGCACGCCTAAAACGTCGGCGTTTTCCACCGCGCTTTTGCCTTCCGCGCGAATAATCGTCGAACCGTCCACGCTGCCCACTAAGGCGCGTAAATACAATTCGCCATCGCGCAACACCGCATAACCGCCAATAGGCACTTGGCAGCCGCCTTGTAAGCGATTGTTCATGGCGCGTTCCGCCAGCACGCAAGAGGTAGTTTCCGCATCCGCCAGTGGCGCCAGTAATTGTTGAACCGCAACATCGTCTGTCCGACATTCGATGCCCACCGCGCCTTGTCCGGCAGCAGGCAAGGATTGCTCCACATCAATAAAAGAAGCGATACGCTCTGCTAAACCTAACCGAATCAAACCGGCGGAAGCCAAAATAATGGCAT  >  1/1538476‑1538961

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