Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,586 2 bp→GC coding (814‑815/963 nt) hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,5860CG100.0% 572.8 / NA 185R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (106/79);  total (106/79)
*NC_0009133,987,5870GC100.0% 572.7 / NA 185R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (106/79);  total (106/79)

Only 100 of 193 total aligned reads displayed.

AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGG  >  NC_000913/3987499‑3987672
                                                                                       ||                                                                                     
aGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcgcgccg                                                                                   <  1:2420188/93‑1 (MQ=255)
aGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcgcgccg                                                                                   >  1:2694100/1‑93 (MQ=255)
 gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                  <  2:4634424/93‑1 (MQ=255)
 gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                  >  1:2329387/1‑93 (MQ=255)
 gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                  >  1:1707503/1‑93 (MQ=255)
 gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                  >  2:712566/1‑93 (MQ=255)
  ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                 <  1:2871231/93‑1 (MQ=255)
   tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                <  1:517542/93‑1 (MQ=255)
    cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCa                                                                               <  1:3901355/93‑1 (MQ=255)
     ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAg                                                                              >  1:2011112/1‑93 (MQ=255)
     ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAg                                                                              >  1:3578077/1‑93 (MQ=255)
        tGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc                                                                           >  2:2204583/1‑93 (MQ=255)
         gCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgca                                                                          <  2:3466929/93‑1 (MQ=255)
          ccAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                         >  1:2540741/1‑93 (MQ=255)
           cAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacg                                                                        >  1:2617663/1‑93 (MQ=255)
            aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGGTCACCGCGAATAATCTGCGCACAGTCCGGCGCGCCCACCACggcacgc                                                                       >  2:4214912/1‑93 (MQ=255)
             gCGAACTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgca                                                                      <  2:53599/93‑1 (MQ=255)
               gAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                    >  1:3391181/1‑93 (MQ=255)
                 aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  >  2:1941813/1‑93 (MQ=255)
                 aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  <  2:4582820/93‑1 (MQ=255)
                  aTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAg                                                                 >  2:2833727/1‑93 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                <  2:2747577/93‑1 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                >  1:1501317/1‑93 (MQ=255)
                     cccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAtt                                                              >  2:33507/1‑93 (MQ=255)
                      ccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAttt                                                             >  2:3668840/1‑93 (MQ=255)
                      ccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAttt                                                             >  2:3285861/1‑93 (MQ=255)
                       cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            >  1:1150527/1‑93 (MQ=255)
                       cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            >  2:2084661/1‑93 (MQ=255)
                        aTTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCg                                                           >  1:812432/1‑93 (MQ=255)
                         tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGc                                                          >  1:1684578/1‑93 (MQ=255)
                          ttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGcc                                                         >  2:139594/1‑93 (MQ=255)
                           tGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCa                                                        >  1:2297912/1‑93 (MQ=255)
                             ttCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATc                                                      >  2:96853/1‑93 (MQ=255)
                              tCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCa                                                     >  1:3142853/1‑93 (MQ=255)
                               cGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                    <  1:1028107/93‑1 (MQ=255)
                                ggCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAAt                                                   <  2:2128570/93‑1 (MQ=255)
                                 gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                  >  2:2898763/1‑93 (MQ=255)
                                 gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                  >  1:1064010/1‑93 (MQ=255)
                                   aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAg                                                <  1:4028833/93‑1 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                               >  2:3035669/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                               >  1:4305010/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                               >  1:1231995/1‑93 (MQ=255)
                                     cTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCt                                              >  1:4164495/1‑93 (MQ=255)
                                      ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             >  1:2964710/1‑93 (MQ=255)
                                      ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             >  2:3105999/1‑93 (MQ=255)
                                      ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             >  1:351239/1‑93 (MQ=255)
                                          ggcgCACCGCGGAGTTCACCGCGAATAATCTGCGAACCATCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCcgcg                                         >  1:2792068/1‑93 (MQ=255)
                                           gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                        >  1:4366991/1‑93 (MQ=255)
                                             gcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg                                      >  1:2172306/1‑93 (MQ=255)
                                              cACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                     <  2:400963/93‑1 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    >  1:1038959/1‑93 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    >  2:1671246/1‑93 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    <  2:3408228/93‑1 (MQ=255)
                                               aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    >  1:3923052/1‑93 (MQ=255)
                                                 cgcgGCGTTCACCGCGAATAATCTGCGAGCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAc                                  >  1:313061/1‑93 (MQ=255)
                                                  gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAcc                                 <  2:1348764/93‑1 (MQ=255)
                                                  gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAcc                                 >  1:882750/1‑93 (MQ=255)
                                                   cggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCa                                <  1:2311260/93‑1 (MQ=255)
                                                   cggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCa                                >  2:3300232/1‑93 (MQ=255)
                                                    ggcgGTCGCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa                               <  1:1399115/93‑1 (MQ=255)
                                                      cgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAtt                             <  2:1814776/93‑1 (MQ=255)
                                                       gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                            <  1:4384062/93‑1 (MQ=255)
                                                       gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                            <  2:2075427/93‑1 (MQ=255)
                                                       gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                            >  2:2330040/1‑93 (MQ=255)
                                                        ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGACTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                           >  1:1976944/1‑93 (MQ=255)
                                                         tCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                          >  2:2143466/1‑93 (MQ=255)
                                                         tCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                          <  2:1000141/93‑1 (MQ=255)
                                                          cACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCa                         <  1:335979/93‑1 (MQ=255)
                                                            ccGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCTATAATCTCGGCGTAGCTACCAATTGGCAcc                       >  2:874484/1‑93 (MQ=255)
                                                            ccGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAcc                       <  1:1118544/93‑1 (MQ=255)
                                                              gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                     <  1:53352/93‑1 (MQ=255)
                                                              gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                     >  2:1727532/1‑93 (MQ=255)
                                                              gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                     <  2:3705489/93‑1 (MQ=255)
                                                              gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                     >  2:103976/1‑93 (MQ=255)
                                                                gAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGAc                   >  2:3458928/1‑93 (MQ=255)
                                                                gAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGAc                   >  2:445506/1‑93 (MQ=255)
                                                                 aataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACa                  >  1:1580619/1‑93 (MQ=255)
                                                                   taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                >  1:3546541/1‑93 (MQ=255)
                                                                   taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                >  2:3788317/1‑93 (MQ=255)
                                                                   taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                <  1:1523024/93‑1 (MQ=255)
                                                                    aatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATcc               <  2:1601901/93‑1 (MQ=255)
                                                                       cTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGcc            <  1:397160/93‑1 (MQ=255)
                                                                        tGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt           >  2:4310318/1‑93 (MQ=255)
                                                                        tGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt           >  2:1250031/1‑93 (MQ=255)
                                                                        tGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt           >  2:2566015/1‑93 (MQ=255)
                                                                         gCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCCCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt          <  1:632740/93‑1 (MQ=255)
                                                                           gAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCg        >  1:1467047/1‑93 (MQ=255)
                                                                             aCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag      >  1:3093271/1‑93 (MQ=255)
                                                                             aCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag      >  2:3737993/1‑93 (MQ=255)
                                                                             aCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag      >  2:1837233/1‑93 (MQ=255)
                                                                              ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     <  2:3191728/93‑1 (MQ=255)
                                                                              ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     >  2:3165238/1‑93 (MQ=255)
                                                                              ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     <  2:157776/93‑1 (MQ=255)
                                                                               cGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAc    <  2:2951061/93‑1 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   >  2:3843795/1‑93 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   <  1:4071411/93‑1 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   >  2:2365761/1‑93 (MQ=255)
                                                                                 tCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAGCTGACATCCGCCTTCGAGACgg  >  1:1639409/1‑93 (MQ=255)
                                                                                 tCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg  <  2:572618/93‑1 (MQ=255)
                                                                                 tCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCCGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg  <  1:2516792/93‑1 (MQ=255)
                                                                                       ||                                                                                     
AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGG  >  NC_000913/3987499‑3987672

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