Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,451 T→C G329G (GGT→GGC cysG → uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,4510TC100.0% 648.8 / NA 204G329G (GGT→GGCcysGuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (108/96);  total (108/96)

Only 100 of 210 total aligned reads displayed.

GAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACA  >  NC_000913/3496361‑3496540
                                                                                          |                                                                                         
gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                                         >  1:4609239/1‑93 (MQ=255)
gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                                         >  1:4174809/1‑93 (MQ=255)
gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAt                                                                                         <  1:1347658/93‑1 (MQ=255)
  aaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTa                                                                                       <  1:458990/93‑1 (MQ=255)
  aaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTa                                                                                       <  2:3016488/93‑1 (MQ=255)
   aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAc                                                                                      >  1:1158427/1‑93 (MQ=255)
   aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAc                                                                                      >  1:2216373/1‑93 (MQ=255)
    ggTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAcc                                                                                     <  2:4260234/93‑1 (MQ=255)
    ggTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAcc                                                                                     <  1:1372172/93‑1 (MQ=255)
    ggTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAcc                                                                                     >  1:2147401/1‑93 (MQ=255)
         cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGc                                                                                <  1:951492/93‑1 (MQ=255)
         cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGc                                                                                <  2:1902499/93‑1 (MQ=255)
           aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCtt                                                                              >  1:388201/1‑93 (MQ=255)
           aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCtt                                                                              >  1:2641449/1‑93 (MQ=255)
            tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                                             <  1:3275727/93‑1 (MQ=255)
            tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                                             <  1:3160487/93‑1 (MQ=255)
            tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                                             <  2:3251754/93‑1 (MQ=255)
            tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                                             >  2:195368/1‑93 (MQ=255)
            tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                                             >  1:2474196/1‑93 (MQ=255)
             ccGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCt                                                                            <  2:4357678/93‑1 (MQ=255)
               gTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTgg                                                                          <  1:3594019/93‑1 (MQ=255)
                 ttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGATTTACGGCAGTTTCTgttg                                                                        >  2:3065940/1‑90 (MQ=255)
                 ttATTTCTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGtt                                                                        >  2:1494771/1‑93 (MQ=255)
                   aTTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGc                                                                      <  2:4019923/93‑1 (MQ=255)
                   aTTGTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGc                                                                      <  1:2025769/93‑1 (MQ=255)
                    tttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGct                                                                     >  2:2441824/1‑93 (MQ=255)
                     ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                                                    >  1:3920525/1‑93 (MQ=255)
                       ttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTg                                                                  >  2:2336597/1‑93 (MQ=255)
                        tGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGc                                                                 <  2:674674/93‑1 (MQ=255)
                         ggCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGcc                                                                >  2:3181599/1‑93 (MQ=255)
                           ccGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTa                                                              >  2:2795907/1‑93 (MQ=255)
                            cGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAt                                                             <  2:4496733/93‑1 (MQ=255)
                              tggtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTc                                                           <  2:2323110/93‑1 (MQ=255)
                              tggtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTc                                                           <  2:713991/93‑1 (MQ=255)
                                gtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCgg                                                         <  2:2867467/93‑1 (MQ=255)
                                gtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCgg                                                         <  1:4049922/93‑1 (MQ=255)
                                     gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                                                    >  2:4073957/1‑93 (MQ=255)
                                      aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                                                   <  2:4633897/93‑1 (MQ=255)
                                       agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                                                  <  2:2952085/93‑1 (MQ=255)
                                       agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                                                  >  1:2901423/1‑93 (MQ=255)
                                       agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                                                  <  1:4449796/93‑1 (MQ=255)
                                        gagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCa                                                 >  1:3179295/1‑93 (MQ=255)
                                        gagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCa                                                 >  2:4421276/1‑93 (MQ=255)
                                        gagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCa                                                 >  1:3029580/1‑93 (MQ=255)
                                         agCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAc                                                >  2:3632745/1‑93 (MQ=255)
                                          gCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACt                                               >  2:2490550/1‑93 (MQ=255)
                                            tGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCa                                             >  2:141027/1‑93 (MQ=255)
                                            tGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCa                                             <  2:746143/93‑1 (MQ=255)
                                             ggAAACACTGTGCAACGCGGGCATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGCTGCTCTGCCTATTAGGGTAGTCCACTCAc                                            >  1:1596906/1‑93 (MQ=255)
                                              gAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACg                                           >  1:2018684/1‑93 (MQ=255)
                                              gAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCACCTTCTGGTTGCTCTGCCTATTCGGGGATTGCACTCACg                                           >  1:2037944/1‑93 (MQ=255)
                                                 acacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCAt                                        >  1:257334/1‑93 (MQ=255)
                                                  cacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                                       >  2:1972099/1‑93 (MQ=255)
                                                    cTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgc                                     >  1:3742507/1‑93 (MQ=255)
                                                    cTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgc                                     <  2:1770117/93‑1 (MQ=255)
                                                    cTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgc                                     <  2:63470/93‑1 (MQ=255)
                                                       tgCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAt                                  <  2:3751976/93‑1 (MQ=255)
                                                       tgCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAt                                  <  1:3376002/93‑1 (MQ=255)
                                                       tgCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAt                                  >  1:1255993/1‑93 (MQ=255)
                                                         cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa                                <  2:1029250/93‑1 (MQ=255)
                                                          aaCGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTAt                               <  2:254401/93‑1 (MQ=255)
                                                          aaCGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTAt                               <  2:144046/93‑1 (MQ=255)
                                                            cgcgGGTGTCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc                             <  2:1684115/93‑1 (MQ=255)
                                                            cgcgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc                             >  2:574114/1‑93 (MQ=255)
                                                              cgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGccc                           <  2:2847734/93‑1 (MQ=255)
                                                              cgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGccc                           >  1:1037014/1‑93 (MQ=255)
                                                                 gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga                        <  1:270687/93‑1 (MQ=255)
                                                                 gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga                        <  1:4486168/93‑1 (MQ=255)
                                                                    tcccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCg                     <  2:1329279/91‑1 (MQ=255)
                                                                    tcccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCg                     >  2:4100419/3‑93 (MQ=255)
                                                                    tcccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCg                     >  1:625590/3‑93 (MQ=255)
                                                                     cccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGt                    <  2:3395735/92‑1 (MQ=255)
                                                                     cccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGt                    >  1:4446519/2‑93 (MQ=255)
                                                                      ctgTTCTCAGTGGTTCCGTGCATTACCGCAGCTTCTGGTTGCTCTGCCTATGCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTa                   <  2:3594699/91‑1 (MQ=255)
                                                                      ccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTa                   >  2:477508/1‑93 (MQ=255)
                                                                      ccGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTa                   <  2:417597/93‑1 (MQ=255)
                                                                        gTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACg                 >  1:3319937/1‑93 (MQ=255)
                                                                        gTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACg                 >  1:4302766/1‑93 (MQ=255)
                                                                          tctcGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCt               <  2:2688186/93‑1 (MQ=255)
                                                                          tctcGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCt               >  2:3936648/1‑93 (MQ=255)
                                                                          tctcGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCt               <  1:4129744/93‑1 (MQ=255)
                                                                           ctcGGTCGTTTCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCtt              <  1:1012034/93‑1 (MQ=255)
                                                                            tcGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTa             <  2:1651112/93‑1 (MQ=255)
                                                                             cGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTaa            >  1:4091480/1‑93 (MQ=255)
                                                                             cGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTaa            >  1:1134366/1‑93 (MQ=255)
                                                                             cGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTaa            <  1:3234682/93‑1 (MQ=255)
                                                                             cGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTaa            <  1:4534484/93‑1 (MQ=255)
                                                                              ggtggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAAt           <  2:1140842/93‑1 (MQ=255)
                                                                               gtggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAAtt          >  1:4683148/1‑93 (MQ=255)
                                                                               gtggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAAtt          <  2:3977351/93‑1 (MQ=255)
                                                                               gtggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAAtt          >  1:1004881/1‑93 (MQ=255)
                                                                                tggtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTa         >  1:746911/1‑93 (MQ=255)
                                                                                  gtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTAcc       <  1:2326028/93‑1 (MQ=255)
                                                                                  gtTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTAcc       <  1:1671887/93‑1 (MQ=255)
                                                                                   ttCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCg      <  1:854073/93‑1 (MQ=255)
                                                                                   ttCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCg      <  1:1876657/93‑1 (MQ=255)
                                                                                     ccGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGa    <  1:950719/93‑1 (MQ=255)
                                                                                      cGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGac   <  2:2423478/93‑1 (MQ=255)
                                                                                      cGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGac   >  1:3990715/1‑93 (MQ=255)
                                                                                       gggCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGaca  >  1:3063326/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACA  >  NC_000913/3496361‑3496540

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