Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,716,271 A→G I129V (ATA→GTA)  unknown → Ferrochelatase, protoheme ferro‑lyase (EC 4.99.1.1)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,716,2710AG100.0% 182.7 / NA 60I129V (ATA→GTA) unknownFerrochelatase, protoheme ferro‑lyase (EC 4.99.1.1)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/27);  total (33/27)

CGCTGCTCAAAGCCATTATTTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTATATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCTTTCCATGGCATTCCGTTACGTTACGAAAACACCGGCGATTATTA  >  1/1716029‑1716512
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
cGCTGCTCAAAGCCATTATTTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCg                                                                                                                                                                                                                                           <  1:386699/251‑1 (MQ=255)
  ctgctCAAAGCCATTATTTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCg                                                                                                                                                                                                                                           >  1:431111/1‑249 (MQ=255)
             cATTATTTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAAt                                                                                                                                                                                                                              <  1:326886/251‑1 (MQ=255)
                 attTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAAt                                                                                                                                                                                                                              >  1:642829/1‑247 (MQ=255)
                 attTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAg                                                                                                                                                                                                                          >  1:726300/1‑251 (MQ=255)
                 attTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTCCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAg                                                                                                                                                                                                                          >  1:891014/1‑251 (MQ=255)
                  ttttGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGt                                                                                                                                                                                                                         <  1:315837/251‑1 (MQ=255)
                    ttGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAg                                                                                                                                                                                                                       <  1:508309/251‑1 (MQ=255)
                     tGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGc                                                                                                                                                                                                                      <  1:725163/251‑1 (MQ=255)
                       cccTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGcac                                                                                                                                                                                                                    <  1:364095/251‑1 (MQ=255)
                       cccTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGcac                                                                                                                                                                                                                    <  1:207207/251‑1 (MQ=255)
                         cTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGcacca                                                                                                                                                                                                                  >  1:316697/1‑251 (MQ=255)
                                 ccAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATAc                                                                                                                                                                                                                            >  1:82295/1‑233 (MQ=255)
                                  cAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGgcg                                                                                                                                                                                                         <  1:322608/251‑1 (MQ=255)
                                      cGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACGAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATAc                                                                                                                                                                                                                            >  1:454242/1‑228 (MQ=255)
                                            gCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGa                                                                                                                                                                                               >  1:530162/1‑251 (MQ=255)
                                             cAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGAc                                                                                                                                                                                              <  1:824046/251‑1 (MQ=255)
                                              aaaaaaCTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACt                                                                                                                                                                                             <  1:719630/251‑1 (MQ=255)
                                                     tACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGc                                                                                                                                                                                      <  1:131364/251‑1 (MQ=255)
                                                                   gACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGaaa                                                                                                                                                                         <  1:267336/250‑1 (MQ=255)
                                                                     ccAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGc                                                                                                                                                                                      >  1:424147/1‑235 (MQ=255)
                                                                      cAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAaacaac                                                                                                                                                                     <  1:551422/251‑1 (MQ=255)
                                                                                tCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCAt                                                                                                                                                           <  1:308112/251‑1 (MQ=255)
                                                                                tCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCAt                                                                                                                                                           <  1:700602/251‑1 (MQ=255)
                                                                                tCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCAt                                                                                                                                                           <  1:25936/251‑1 (MQ=255)
                                                                                tCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCAt                                                                                                                                                           <  1:884413/251‑1 (MQ=255)
                                                                                          ttGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTc                                                                                                                                                 <  1:679525/251‑1 (MQ=255)
                                                                                           tGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCg                                                                                                                                                >  1:716134/1‑251 (MQ=255)
                                                                                                tAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGC‑CTGAAACAACAGCGTGGCATCGTGCCGTTCGATttca                                                                                                                                          >  1:158758/1‑251 (MQ=255)
                                                                                                           caaAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTa                                                                                                                                >  1:265165/1‑251 (MQ=255)
                                                                                                                                tttCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGAAGCACCGCCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTatat                                                                                                           <  1:345574/251‑1 (MQ=255)
                                                                                                                                  tCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTa                                                                                                         <  1:783826/251‑1 (MQ=255)
                                                                                                                                         caacaaAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATAtt                                                                                                          >  1:4678/1‑243 (MQ=255)
                                                                                                                                         caacaaAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGcc                                                                                                    >  1:737217/1‑249 (MQ=255)
                                                                                                                                                aCCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCtt                                                                                                                                                                                         >  1:408468/1‑157 (MQ=255)
                                                                                                                                                aCCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATttcattca                                                                                                                                      >  1:890750/1‑208 (MQ=255)
                                                                                                                                                                            gACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAag                                                               >  1:606313/1‑251 (MQ=255)
                                                                                                                                                                             aCTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGTGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAaga                                                              >  1:655961/1‑251 (MQ=255)
                                                                                                                                                                             aCTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCAt                                                                                                                                                           <  1:276740/158‑1 (MQ=255)
                                                                                                                                                                             aCTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGAttt                                                                                                                                            >  1:82448/1‑173 (MQ=255)
                                                                                                                                                                             aCTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGccc                                                                                                   >  1:615855/1‑214 (MQ=255)
                                                                                                                                                                             aCTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAaga                                                              >  1:372641/1‑251 (MQ=255)
                                                                                                                                                                                         ccGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATc                                                                                         >  1:791036/1‑212 (MQ=255)
                                                                                                                                                                                           gTCAATACCCGGCGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGaaa                                                                                                                                                                         <  1:858269/130‑1 (MQ=255)
                                                                                                                                                                                           gTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGgc                                                                                                                                                                                                          <  1:517831/97‑1 (MQ=255)
                                                                                                                                                                                           gTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAAcc                                                                                                                     >  1:284281/1‑182 (MQ=255)
                                                                                                                                                                                               aTACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTgtat                                                                                                                                                                                                                                                <  1:371763/55‑1 (MQ=255)
                                                                                                                                                                                               aTACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCttt                                            >  1:481780/1‑251 (MQ=255)
                                                                                                                                                                                                           gTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCg                                                                                                                                                                                                                                           >  1:405335/1‑48 (MQ=255)
                                                                                                                                                                                                           gTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATcc                                                                                        <  1:71313/195‑1 (MQ=255)
                                                                                                                                                                                                           gTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTAc                                                      >  1:47495/1‑229 (MQ=255)
                                                                                                                                                                                                           gTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCTTTCCATGGCATTc                                 >  1:544808/1‑250 (MQ=255)
                                                                                                                                                                                                           gTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGATTATATTACAGCCCTGATTGAATCCTGCaa                                                                                   >  1:450242/1‑200 (MQ=255)
                                                                                                                                                                                                                aaaaTTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATc                                                                    >  1:491909/1‑210 (MQ=255)
                                                                                                                                                                                                                  aaTTAATAGATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCTTTCCATGGCATTCcgttacgt                         <  1:537555/251‑1 (MQ=255)
                                                                                                                                                                                                                          gATAAAAAAGTGGATCACATTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGAc                                                                                                                                                                                              >  1:583554/1‑78 (MQ=255)
                                                                                                                                                                                                                                          acaTTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCTTTCCATGGCATTCCGTTACGTTACGAAAACACCGGCGAttatta  >  1:399695/1‑250 (MQ=255)
                                                                                                                                                                                                                                          acaTTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCTTTCCATGGCATTCCGTTACGTTACGAAAACACCGGCGAttatta  >  1:745096/1‑250 (MQ=255)
                                                                                                                                                                                                                                           caTTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATAtt                                                                                                          <  1:401465/145‑1 (MQ=255)
                                                                                                                                                                                                                                           caTTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCaa                                                                                   <  1:848633/168‑1 (MQ=255)
                                                                                                                                                                                                                                           caTTATTGTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCttt                                            >  1:475250/1‑207 (MQ=255)
                                                                                                                                                                                                                                             ttattgTATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGAtt                                                                                              <  1:247968/155‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
CGCTGCTCAAAGCCATTATTTTGCCCTTGCGCTCCAAACGTATTGCAAAAAACTACCGTGCCATTTGGACCAAACAAGGTTCGCCGCTACTTGCCATAACTCAACAACAAAAAGAAGCGTTAGCGATTTTTCTGCAACAACAAAACCTTAATATGGAAATCGAAATCGGCATGACTTACGGCAATCCGTCAATACCCGACGCGGTACAAAAATTAATAGATAAAAAAGTGGATCACATTATTATATTGCCGCTTTATCCGCAATACAGTAGCACCACCACAGGCGCTGTGTTTGACTCCTTCGCTGCGGCGCTGAAACAACAGCGTGGCATCGTGCCGTTCGATTTCATTCATTCTTATCAATTAAACCAAGGTTATATTACAGCCCTGATTGAATCCTGCAAAGTGCGGTTGAAATCCGAAGAGTTTTTACTGTTCTCTTTCCATGGCATTCCGTTACGTTACGAAAACACCGGCGATTATTA  >  1/1716029‑1716512

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