Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,822 2 bp→CG coding (578‑579/963 nt) hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,8220GC100.0% 519.5 / NA 160G193A (GGC→GCGhemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (85/75);  total (85/75)
*NC_0009133,987,8230CG100.0% 514.1 / NA 160G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (85/75);  total (85/75)

Only 100 of 167 total aligned reads displayed.

AGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTT  >  NC_000913/3987736‑3987908
                                                                                      ||                                                                                     
agcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGGGGAAACgcggcg                                                                                  >  2:2365155/1‑93 (MQ=255)
   tCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGa                                                                               >  1:4745663/1‑93 (MQ=255)
   tCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGa                                                                               >  2:2248375/1‑93 (MQ=255)
    cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                              >  1:1056794/1‑93 (MQ=255)
    cgcgAGGGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                              >  1:1305934/1‑93 (MQ=255)
      cgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATa                                                                            >  2:3378300/1‑93 (MQ=255)
         gTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGt                                                                         <  2:3061361/93‑1 (MQ=255)
         gTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGt                                                                         <  2:2914249/93‑1 (MQ=255)
           gCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAATAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGa                                                                       <  2:4348034/93‑1 (MQ=255)
           gCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGa                                                                       <  1:3015575/93‑1 (MQ=255)
            cGTGAATCATCAAGGCGGCAGTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                                      <  2:2845725/93‑1 (MQ=255)
               gAATCATCAAGGTGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCACCGCGGCGCGAATACGTGACTcc                                                                   >  2:224302/1‑93 (MQ=255)
                aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                  <  2:3576273/93‑1 (MQ=255)
                 atcatcAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                                 >  1:238907/1‑93 (MQ=255)
                  tcatcaAGGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                                <  1:1042199/93‑1 (MQ=255)
                  tcatcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                                <  1:4227501/93‑1 (MQ=255)
                  tcatcaAGGCGGCATTAAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                                <  1:1442309/93‑1 (MQ=255)
                       aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                           >  1:2874470/1‑93 (MQ=255)
                       aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                           >  1:611308/1‑93 (MQ=255)
                        aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                          >  1:157771/1‑93 (MQ=255)
                        aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                          >  1:3829495/1‑93 (MQ=255)
                        aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAACTaac                                                          >  1:434438/1‑92 (MQ=255)
                         agcggaATTCAATACCCACCGCAACTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAc                                                         <  2:349341/92‑1 (MQ=255)
                           cggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGt                                                       <  2:639696/93‑1 (MQ=255)
                            ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                                      >  1:1926082/1‑93 (MQ=255)
                            ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                                      >  1:1608491/1‑93 (MQ=255)
                            ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                                      <  2:3685774/93‑1 (MQ=255)
                              cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                                    >  2:1322403/1‑93 (MQ=255)
                              cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                                    <  2:3165581/93‑1 (MQ=255)
                              cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                                    >  2:2235312/1‑93 (MQ=255)
                               aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACCTGACTCCAGACCTAAACGTTTTa                                                   >  2:2606173/1‑93 (MQ=255)
                                ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGGCCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAGACGTTTTAg                                                  >  2:3701376/1‑93 (MQ=255)
                                ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGTCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                  >  2:170774/1‑93 (MQ=255)
                                 tCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                 <  1:1844346/93‑1 (MQ=255)
                                  taaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                <  1:4295297/92‑1 (MQ=255)
                                  cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                <  2:2041261/93‑1 (MQ=255)
                                  cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                <  2:3301684/93‑1 (MQ=255)
                                   aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                               >  1:1027903/1‑93 (MQ=255)
                                   aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                               <  1:1989106/93‑1 (MQ=255)
                                    aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCg                                              <  2:119460/93‑1 (MQ=255)
                                      aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                            >  2:3483624/1‑93 (MQ=255)
                                      aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                            <  1:2968672/93‑1 (MQ=255)
                                      aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                            <  1:2968671/93‑1 (MQ=255)
                                        ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                          >  2:482760/1‑93 (MQ=255)
                                        ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTACGGCTa                                          >  2:1438939/1‑93 (MQ=255)
                                         cACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                         >  2:2873719/1‑93 (MQ=255)
                                         cACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                         <  2:1993108/93‑1 (MQ=255)
                                         cACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                         <  2:4061571/93‑1 (MQ=255)
                                          aCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                                        >  1:496538/1‑93 (MQ=255)
                                          aCCGCACCTTGTCCTACCGCCAGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGGTTTAGTACGGCTACg                                        >  2:3403036/1‑93 (MQ=255)
                                           ccGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                                       >  2:1393886/1‑93 (MQ=255)
                                            cGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                                      <  2:3871299/93‑1 (MQ=255)
                                            cGCACCGTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                                      >  2:4606930/1‑93 (MQ=255)
                                             gCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                                     >  1:951918/1‑93 (MQ=255)
                                              cACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCaa                                    >  2:3204694/1‑93 (MQ=255)
                                                ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                  <  1:3747302/93‑1 (MQ=255)
                                                 cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                                 <  1:4002470/93‑1 (MQ=255)
                                                 cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                                 >  1:3472491/1‑93 (MQ=255)
                                                  ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                >  2:799127/1‑93 (MQ=255)
                                                  ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                <  2:284702/93‑1 (MQ=255)
                                                  ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                <  2:4338341/93‑1 (MQ=255)
                                                   tGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                               <  1:3583204/93‑1 (MQ=255)
                                                   tGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                               >  1:1226790/1‑93 (MQ=255)
                                                     tcatACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                             <  2:2531915/90‑1 (MQ=255)
                                                     tCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                             <  1:2143085/93‑1 (MQ=255)
                                                     tCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGAGCCTCCAGACCTAAACGCTTTAGTCCGGCTAAGGCAAGAatgac                             >  2:1770598/1‑92 (MQ=255)
                                                      ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                            <  2:1732562/93‑1 (MQ=255)
                                                      ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                            <  2:3589567/93‑1 (MQ=255)
                                                      ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                            <  2:40067/93‑1 (MQ=255)
                                                      ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                            >  2:1651172/1‑93 (MQ=255)
                                                      ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                            <  2:1742532/93‑1 (MQ=255)
                                                       cTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATgg                           >  2:1321798/1‑93 (MQ=255)
                                                       cTACCGCCGGAAGAGAGATCTCGGGTGGCAACGAGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATgg                           >  1:4314415/1‑93 (MQ=255)
                                                        tACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                          >  2:4084913/1‑93 (MQ=255)
                                                         aCCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCa                         <  2:1877890/93‑1 (MQ=255)
                                                         aCCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCa                         >  1:82913/1‑93 (MQ=255)
                                                           cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                       <  2:2550930/93‑1 (MQ=255)
                                                           cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                       <  1:4651312/93‑1 (MQ=255)
                                                            gccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg                      <  1:2895429/93‑1 (MQ=255)
                                                             ccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt                     >  2:4111970/1‑93 (MQ=255)
                                                              cgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa                    >  2:1664807/1‑93 (MQ=255)
                                                               ggAATAGAGATCTCGTGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTTTGGCTACGTCAAGAATGATGGCATCGTAt                   >  2:4333209/1‑93 (MQ=255)
                                                               ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                   <  2:2659492/93‑1 (MQ=255)
                                                               ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                   >  1:3180852/1‑93 (MQ=255)
                                                               ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGCCTACGGCAGGAATGATGGCATCGTAt                   >  1:4027849/1‑93 (MQ=255)
                                                                  agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAGGTTTTAGTCCGGCTACGTCACGAATAAAGGCATCGTATCCg                >  2:357394/1‑93 (MQ=255)
                                                                   gagagaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc               >  1:2783762/1‑93 (MQ=255)
                                                                   gagagaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc               >  2:4390999/1‑93 (MQ=255)
                                                                    agagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGcc              <  1:3373805/93‑1 (MQ=255)
                                                                    agagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGcc              >  1:222025/1‑93 (MQ=255)
                                                                       gaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGtt           <  2:3818405/92‑1 (MQ=255)
                                                                       gaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGtt           <  1:2584751/92‑1 (MQ=255)
                                                                           tcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATcc       <  2:1189264/93‑1 (MQ=255)
                                                                           tcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATcc       <  1:1707997/93‑1 (MQ=255)
                                                                            cGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCa      >  2:2484192/1‑93 (MQ=255)
                                                                             gggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAg     >  2:155336/1‑93 (MQ=255)
                                                                              ggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGt    <  1:446552/93‑1 (MQ=255)
                                                                              ggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGt    <  2:1033044/93‑1 (MQ=255)
                                                                                tGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGCTGGCATCGTATTCGCCGTTATCCAGttt  >  1:2279129/1‑93 (MQ=255)
                                                                                tGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGttt  >  1:2564801/1‑93 (MQ=255)
                                                                                      ||                                                                                     
AGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTT  >  NC_000913/3987736‑3987908

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