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% 920.3 / NA 290R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (177/113);  total (177/113)
*NC_0009133,987,5870GC100.0% 917.1 / NA 290R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (177/113);  total (177/113)

Only 100 of 309 total aligned reads displayed.

AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGG  >  NC_000913/3987499‑3987673
                                                                                       ||                                                                                      
aGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcgcgccg                                                                                    >  2:7578328/1‑93 (MQ=255)
 gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCTAACCGTCCGGCGCGCCGa                                                                                   >  2:8511195/1‑93 (MQ=255)
  ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                  >  1:9646730/1‑93 (MQ=255)
   tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                 <  1:559964/93‑1 (MQ=255)
   tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                 >  2:8834541/1‑93 (MQ=255)
   tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                 >  2:6065951/1‑93 (MQ=255)
    cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCa                                                                                >  1:148825/1‑93 (MQ=255)
    cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCa                                                                                >  2:8449826/1‑93 (MQ=255)
    cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCa                                                                                <  1:8898791/93‑1 (MQ=255)
      tctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgc                                                                              >  2:6174506/1‑93 (MQ=255)
      tctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgc                                                                              <  1:7828202/93‑1 (MQ=255)
      tctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgc                                                                              >  1:4451302/1‑93 (MQ=255)
      tctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgc                                                                              <  1:4593842/93‑1 (MQ=255)
       ctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcg                                                                             >  1:6573163/1‑93 (MQ=255)
         gCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGCCCAGcgca                                                                           >  1:5834214/1‑93 (MQ=255)
         gCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgca                                                                           >  2:2700767/1‑93 (MQ=255)
          ccAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                          >  1:3364755/1‑93 (MQ=255)
           cAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGAGAACCGTCCGGCGCGCCGACCAGcgcacg                                                                         >  1:5276913/1‑93 (MQ=255)
            aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgc                                                                        >  2:2902039/1‑93 (MQ=255)
            aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgc                                                                        <  2:6549428/93‑1 (MQ=255)
             gCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgca                                                                       >  1:306413/1‑93 (MQ=255)
             gCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgca                                                                       >  1:5991458/1‑93 (MQ=255)
                   tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                 >  2:7268000/1‑93 (MQ=255)
                    ccccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCACCGCACGCACCCAgct                                                                >  1:2828408/1‑91 (MQ=255)
                       cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                             <  1:8183398/93‑1 (MQ=255)
                       cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                             >  2:974617/1‑93 (MQ=255)
                        aTTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGAGAACCGTTCGGCGCGCCGACCAGCGCACTTAGCCAGATTTCg                                                            >  2:8281040/1‑93 (MQ=255)
                         tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCCGATTTCGc                                                           >  2:1204186/1‑93 (MQ=255)
                         tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGc                                                           >  1:6877648/1‑93 (MQ=255)
                         tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGc                                                           >  2:7685171/1‑93 (MQ=255)
                           tGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCa                                                         <  1:6285610/93‑1 (MQ=255)
                           tGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCa                                                         >  1:1114260/1‑93 (MQ=255)
                            gTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                        >  2:9232504/1‑93 (MQ=255)
                              tCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCa                                                      <  2:4164481/93‑1 (MQ=255)
                               cGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATcca                                                     >  2:4993974/1‑91 (MQ=255)
                                ggCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAAt                                                    >  1:5500429/1‑93 (MQ=255)
                                 gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                   >  1:3505694/1‑93 (MQ=255)
                                  cATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATaa                                                  <  1:5782177/93‑1 (MQ=255)
                                  cATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATaa                                                  >  1:9399857/1‑93 (MQ=255)
                                  cATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTAGCCATCAATaa                                                  >  2:6606024/1‑93 (MQ=255)
                                   aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAg                                                 >  1:8000268/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCGCCGCGAATAATCTGCGAACCGTCCGGCTCGCCGACCAGCGCCCGCAGCCAGTTTTCGACATCAATAAGc                                                <  2:6910007/93‑1 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                                >  2:5962217/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                                >  2:3641078/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                                >  1:7253220/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGc                                                >  2:196589/1‑93 (MQ=255)
                                    tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGACAGATTTCGCCATCAATAAGc                                                >  1:4921753/1‑93 (MQ=255)
                                     cTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCt                                               <  1:6990840/93‑1 (MQ=255)
                                     cTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCt                                               >  2:153791/1‑93 (MQ=255)
                                        gcggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCgg                                            >  2:3419001/1‑93 (MQ=255)
                                         cggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGc                                           >  2:7427111/1‑93 (MQ=255)
                                         cggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGc                                           <  2:1692377/93‑1 (MQ=255)
                                         cggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGc                                           <  2:2223630/93‑1 (MQ=255)
                                           gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                         >  1:9177260/1‑93 (MQ=255)
                                           gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                         <  1:4874095/93‑1 (MQ=255)
                                            cgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                        <  2:6821106/93‑1 (MQ=255)
                                              cACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                      >  1:7558339/1‑93 (MQ=255)
                                              cACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                      >  1:4633201/1‑93 (MQ=255)
                                              cACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                      >  2:394446/1‑93 (MQ=255)
                                                ccGCGGCGTTCACCGCGAATAATCTGCGGACCGTCCGGCGCGCCGACCAGGGCACGCAGCCAGATTTGGCCATCAATAAGCTCGTCGGAGCTa                                    >  2:6512552/1‑93 (MQ=255)
                                                ccGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGCTGTCGCCATCAATAAGCTCGGCGTAGcga                                    >  1:9003753/1‑91 (MQ=255)
                                                  gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAcc                                  >  2:704944/1‑93 (MQ=255)
                                                   cggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCa                                 >  1:2285925/1‑93 (MQ=255)
                                                   cggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCa                                 <  1:3819513/93‑1 (MQ=255)
                                                    ggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa                                <  1:535235/93‑1 (MQ=255)
                                                    ggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCCCCGACCAGCACACGCAGCCAGATTTCGCCATCAATAAGCTCGGCTTAGCTACCaa                                <  2:3960556/93‑1 (MQ=255)
                                                      cgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAtt                              <  1:9949133/93‑1 (MQ=255)
                                                      cgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAtt                              >  1:1972262/1‑93 (MQ=255)
                                                       gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                             <  1:4318392/93‑1 (MQ=255)
                                                       gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                             >  1:1840276/1‑93 (MQ=255)
                                                       gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                             >  1:156836/1‑93 (MQ=255)
                                                        ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                            >  1:359757/1‑93 (MQ=255)
                                                        ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                            >  2:5063249/1‑93 (MQ=255)
                                                        ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                            >  2:3161581/1‑93 (MQ=255)
                                                        ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                            >  1:780932/1‑93 (MQ=255)
                                                         tCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCGTCAATAAGCTCGGCGGAGCTACCAATTGGc                           >  2:306641/1‑93 (MQ=255)
                                                         tCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                           <  2:5213170/93‑1 (MQ=255)
                                                         tCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                           >  2:3872405/1‑93 (MQ=255)
                                                           aCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                         >  2:509067/1‑93 (MQ=255)
                                                           aCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                         <  2:3330012/93‑1 (MQ=255)
                                                           aCCGCGAAGAATCTGCGAACCGTCCGGCGCTCCGACCAGCTCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                         <  2:4552308/93‑1 (MQ=255)
                                                            ccGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAcc                        >  2:1844906/1‑93 (MQ=255)
                                                              gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                      >  2:1341559/1‑93 (MQ=255)
                                                              gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                      >  1:1272604/1‑93 (MQ=255)
                                                               cgAATAATCTGCGAGCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGa                     <  2:5772355/93‑1 (MQ=255)
                                                               cgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGa                     <  2:9507456/93‑1 (MQ=255)
                                                                   taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                 >  2:4852806/1‑93 (MQ=255)
                                                                    aatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATcc                <  1:9774583/93‑1 (MQ=255)
                                                                    aatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATcc                >  1:4784101/1‑93 (MQ=255)
                                                                       cTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGcc             <  1:2107968/93‑1 (MQ=255)
                                                                        tGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt            <  1:9452794/93‑1 (MQ=255)
                                                                         gCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt           >  2:9320716/1‑93 (MQ=255)
                                                                           gAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCg         >  2:8023047/1‑93 (MQ=255)
                                                                            aaCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCga        >  2:8939916/1‑93 (MQ=255)
                                                                            aaCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCga        <  2:2435388/93‑1 (MQ=255)
                                                                             aCCGTCCGGCGCGCCGACCAGCGCATGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag       >  2:8172226/1‑93 (MQ=255)
                                                                              ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga      <  2:7533014/93‑1 (MQ=255)
                                                                                gTCCGGCGCGCCGACCAGCGCACGCAGCCAGCTTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg    >  1:850153/1‑93 (MQ=255)
                                                                                 gccggcgcgCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg   <  1:2391207/92‑1 (MQ=255)
                                                                                  ccggcgcgCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACggg  <  1:2740650/93‑1 (MQ=255)
                                                                                       ||                                                                                      
AGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGG  >  NC_000913/3987499‑3987673

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