Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,985,741 T→C V313V (GTA→GTG hemX ← uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,985,7410TC100.0% 737.1 / NA 232V313V (GTA→GTGhemXuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (112/120);  total (112/120)

Only 100 of 248 total aligned reads displayed.

AAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985651‑3985830
                                                                                          |                                                                                         
aaGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAgcacgcac                                                                                         >  2:410691/1‑93 (MQ=255)
  gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       >  1:1222867/1‑93 (MQ=255)
    tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCGCCCAg                                                                                     <  2:3217189/93‑1 (MQ=255)
    tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     <  1:4204225/93‑1 (MQ=255)
    tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     <  1:1875999/93‑1 (MQ=255)
     ccATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAgg                                                                                    <  2:5676856/93‑1 (MQ=255)
      cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                   <  2:266437/93‑1 (MQ=255)
       aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  >  1:1704871/1‑93 (MQ=255)
       aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  >  1:1818146/1‑93 (MQ=255)
       aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  >  2:2539644/1‑93 (MQ=255)
        tCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTga                                                                                 >  1:1747459/1‑92 (MQ=255)
         cGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGga                                                                                <  2:1492403/93‑1 (MQ=255)
         cGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGga                                                                                >  1:6235752/1‑93 (MQ=255)
           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              <  2:2835010/93‑1 (MQ=255)
           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              >  2:6167658/1‑93 (MQ=255)
            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTTGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                             >  2:4166375/1‑93 (MQ=255)
            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                             >  2:5800487/1‑93 (MQ=255)
              tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                           >  2:5473878/1‑93 (MQ=255)
              tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                           >  2:6017739/1‑93 (MQ=255)
               atTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGtt                                                                          <  2:5429053/93‑1 (MQ=255)
                ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                         <  2:155656/93‑1 (MQ=255)
                ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGGGGTGGCATCATCAGTATCGTAGTAAGAACCCACCCAGGGGGAGACGTtc                                                                         >  2:4340617/1‑93 (MQ=255)
                 tttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtct                                                                        <  1:1313084/93‑1 (MQ=255)
                  ttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtctc                                                                       <  1:2433256/93‑1 (MQ=255)
                   tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                      >  1:3873544/1‑93 (MQ=255)
                    tcTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     <  1:618036/92‑1 (MQ=255)
                    gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     <  1:5294362/93‑1 (MQ=255)
                     cTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAg                                                                    <  2:1234310/93‑1 (MQ=255)
                     cTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAg                                                                    <  1:2574510/93‑1 (MQ=255)
                      tGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgc                                                                   >  2:2806049/1‑93 (MQ=255)
                      tGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgc                                                                   >  1:5835127/1‑93 (MQ=255)
                         cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                                >  2:1699948/1‑93 (MQ=255)
                         cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                                <  2:5322410/93‑1 (MQ=255)
                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               >  1:3217191/1‑93 (MQ=255)
                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               <  2:1626543/93‑1 (MQ=255)
                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               >  2:5708261/1‑93 (MQ=255)
                           tAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTg                                                              <  2:76652/93‑1 (MQ=255)
                             aCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGc                                                            >  1:161276/1‑93 (MQ=255)
                               tGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGa                                                          <  1:2553727/93‑1 (MQ=255)
                                ggTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                         <  2:3277370/93‑1 (MQ=255)
                                 gTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATa                                                        >  2:1093313/1‑93 (MQ=255)
                                  tCCACCTCGTCGAGGAACGCTTTGATGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATaa                                                       <  1:4381078/93‑1 (MQ=255)
                                   ccACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAg                                                      >  1:3430409/1‑93 (MQ=255)
                                   ccACCTCGTCGAGGAAAGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAg                                                      >  2:301634/1‑93 (MQ=255)
                                    cACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGt                                                     <  1:2741414/93‑1 (MQ=255)
                                     cccTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    <  2:1997245/92‑1 (MQ=255)
                                       cTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtctc                                                  <  1:1779289/93‑1 (MQ=255)
                                       cTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtctc                                                  >  2:5168731/1‑93 (MQ=255)
                                       cTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtctc                                                  >  2:4832735/1‑93 (MQ=255)
                                        tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 >  1:5097924/1‑93 (MQ=255)
                                          gtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTc                                               <  2:1827457/93‑1 (MQ=255)
                                          gtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTc                                               >  2:339907/1‑93 (MQ=255)
                                           tcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTcc                                              <  1:4609808/93‑1 (MQ=255)
                                             gAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTg                                            <  1:5925262/93‑1 (MQ=255)
                                              aagAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCCCCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                           <  1:5613210/91‑1 (MQ=255)
                                              aGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                           >  1:5511301/1‑93 (MQ=255)
                                              aGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                           >  1:828871/1‑93 (MQ=255)
                                              aGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                           <  2:5534455/93‑1 (MQ=255)
                                               ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          >  2:4959848/1‑93 (MQ=255)
                                               ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          >  1:2405979/1‑93 (MQ=255)
                                                 aaCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGa                                        <  1:4950994/93‑1 (MQ=255)
                                                  aCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                       >  1:3941593/1‑93 (MQ=255)
                                                    gCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcgc                                     <  1:1896133/93‑1 (MQ=255)
                                                      tttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCgg                                   >  2:4532552/1‑93 (MQ=255)
                                                      tttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCgg                                   >  1:609777/1‑93 (MQ=255)
                                                      tttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCgg                                   >  1:4731941/1‑93 (MQ=255)
                                                         ggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTAc                                >  2:132853/1‑93 (MQ=255)
                                                         ggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTAc                                <  2:1265148/93‑1 (MQ=255)
                                                          gtggtggCATCATCAGTATCGTAGTAAGCACGCGCCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACa                               <  1:6001213/93‑1 (MQ=255)
                                                          gtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACa                               <  2:326643/93‑1 (MQ=255)
                                                            ggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGc                             <  1:3684731/93‑1 (MQ=255)
                                                             gtggGATCATCTGTAGCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCt                            <  2:2382363/93‑1 (MQ=255)
                                                              tggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                           >  1:1067550/1‑93 (MQ=255)
                                                              tggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                           >  2:1443351/1‑93 (MQ=255)
                                                               ggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtg                          >  2:4948069/1‑93 (MQ=255)
                                                                 catcatCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgtg                        <  1:1211490/93‑1 (MQ=255)
                                                                 catcatCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgtg                        >  2:2156926/1‑93 (MQ=255)
                                                                 catcatCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgtg                        <  2:1012122/93‑1 (MQ=255)
                                                                  atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       <  1:1402484/93‑1 (MQ=255)
                                                                  atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       >  2:202317/1‑93 (MQ=255)
                                                                   tcatcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgct                      >  1:3502005/1‑93 (MQ=255)
                                                                   tcatcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgct                      <  2:6253658/93‑1 (MQ=255)
                                                                      tcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCg                   >  1:3613195/1‑93 (MQ=255)
                                                                        aGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGCTAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAc                 >  2:2076820/1‑93 (MQ=255)
                                                                         gTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAcc                >  1:1465306/1‑93 (MQ=255)
                                                                          tATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACca               <  1:3264468/93‑1 (MQ=255)
                                                                            tCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagc             >  2:2638975/1‑93 (MQ=255)
                                                                            tCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagc             <  2:3519698/93‑1 (MQ=255)
                                                                            tCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagc             >  1:372156/1‑93 (MQ=255)
                                                                             cGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagca            <  1:711757/93‑1 (MQ=255)
                                                                             cGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagca            <  1:202230/93‑1 (MQ=255)
                                                                               tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg          >  2:4734247/1‑93 (MQ=255)
                                                                                agtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGc         <  2:2625148/93‑1 (MQ=255)
                                                                                agtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGc         >  1:2222523/1‑93 (MQ=255)
                                                                                agtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGc         >  2:932466/1‑93 (MQ=255)
                                                                                  taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       <  1:5408555/93‑1 (MQ=255)
                                                                                  taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       <  2:5184590/93‑1 (MQ=255)
                                                                                   aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      <  1:3130391/93‑1 (MQ=255)
                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga     <  2:2293194/93‑1 (MQ=255)
                                                                                     gcacgcacCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGagag    <  2:3791878/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985651‑3985830

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