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% 1059.7 / NA 329V313V (GTA→GTGhemXuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (158/171);  total (158/171)

Only 100 of 343 total aligned reads displayed.

AAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985651‑3985830
                                                                                          |                                                                                         
 aGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAcc                                                                                        <  1:5247486/93‑1 (MQ=255)
 aGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAcc                                                                                        >  1:4545404/1‑93 (MQ=255)
 aGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAcc                                                                                        >  2:6552221/1‑93 (MQ=255)
  gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       <  2:6328936/93‑1 (MQ=255)
  gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       >  1:1570973/1‑93 (MQ=255)
  gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       >  1:4763419/1‑93 (MQ=255)
  gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       <  1:3694616/93‑1 (MQ=255)
   aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCa                                                                                      <  2:3850475/93‑1 (MQ=255)
   aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCa                                                                                      <  2:2717359/93‑1 (MQ=255)
   aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCa                                                                                      >  2:5511445/1‑93 (MQ=255)
    tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     <  1:2026531/93‑1 (MQ=255)
    tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     >  1:3832604/1‑93 (MQ=255)
     ccATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAgg                                                                                    >  1:5827366/1‑93 (MQ=255)
     ccATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAgg                                                                                    >  1:1893821/1‑93 (MQ=255)
      cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTGGTAAGCACGCACCCAGGt                                                                                   >  2:1038453/1‑93 (MQ=255)
      cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                   >  2:5403567/1‑93 (MQ=255)
      cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                   <  2:3284534/93‑1 (MQ=255)
       ctCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  <  1:3657086/92‑1 (MQ=255)
       aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  >  2:5811891/1‑93 (MQ=255)
        tCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTgg                                                                                 <  2:2150674/93‑1 (MQ=255)
        tCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTgg                                                                                 >  2:4091946/1‑93 (MQ=255)
        tCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTgg                                                                                 >  1:2368784/1‑93 (MQ=255)
          gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACGCAGGTGgag                                                                               >  1:5120212/1‑93 (MQ=255)
          gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgag                                                                               >  1:4941124/1‑93 (MQ=255)
          gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgag                                                                               >  1:4247013/1‑93 (MQ=255)
           ggaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              <  2:5342955/92‑1 (MQ=255)
           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              >  1:2727502/1‑93 (MQ=255)
           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              >  1:1685134/1‑93 (MQ=255)
           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              <  1:5737858/93‑1 (MQ=255)
            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                             >  1:546619/1‑93 (MQ=255)
            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                             >  1:4154495/1‑93 (MQ=255)
             atatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACg                                                                            >  1:1063701/1‑93 (MQ=255)
              tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                           >  1:5943757/1‑93 (MQ=255)
              tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                           <  1:6642999/93‑1 (MQ=255)
               atTTTGCTGGCTTAACTGGTCCACCTCGGCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGtt                                                                          <  2:2963453/93‑1 (MQ=255)
                ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                         <  1:5985388/93‑1 (MQ=255)
                   tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                      <  2:2468618/93‑1 (MQ=255)
                    gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     >  2:131783/1‑93 (MQ=255)
                    gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     <  2:822004/93‑1 (MQ=255)
                      tGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgc                                                                   >  2:2348954/1‑93 (MQ=255)
                      tGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgc                                                                   >  2:2066771/1‑93 (MQ=255)
                         cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                                >  1:2611722/1‑93 (MQ=255)
                         cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                                <  2:5674804/93‑1 (MQ=255)
                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTTGCATCATCAGGATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               >  1:3379633/1‑93 (MQ=255)
                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               <  1:5604577/93‑1 (MQ=255)
                          ttAACTGATCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               <  2:86623/93‑1 (MQ=255)
                            aaCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTgg                                                             >  1:1153440/1‑93 (MQ=255)
                            aaCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTgg                                                             >  2:708274/1‑93 (MQ=255)
                              cTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCg                                                           >  1:1316904/1‑93 (MQ=255)
                                ggTCCACCTCGTCGAGGAACTCTTTGGTGGTGTCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                         <  2:6204659/93‑1 (MQ=255)
                                ggTCCACCTCGTCGAGGAACGCTTTGGGGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                         >  1:294434/1‑93 (MQ=255)
                                 gTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGGGGAGCCGTTCTCCAGCGCCTGGCGATa                                                        >  1:1838713/1‑93 (MQ=255)
                                  tCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATaa                                                       >  2:4718586/1‑93 (MQ=255)
                                     aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    >  2:3540481/1‑93 (MQ=255)
                                     aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    >  2:3665811/1‑93 (MQ=255)
                                       cTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtctc                                                  <  1:3970938/93‑1 (MQ=255)
                                        tcgtctAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 <  2:4643233/93‑1 (MQ=255)
                                        tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 <  1:2274687/93‑1 (MQ=255)
                                        tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 <  2:5977130/93‑1 (MQ=255)
                                         cgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCtt                                                >  2:4391421/1‑93 (MQ=255)
                                           tcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTcc                                              <  1:6401175/93‑1 (MQ=255)
                                            cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCt                                             >  1:244722/1‑93 (MQ=255)
                                            cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCt                                             <  2:2872451/93‑1 (MQ=255)
                                              aGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                           <  2:6121758/93‑1 (MQ=255)
                                               ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          <  1:5388791/93‑1 (MQ=255)
                                               ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          >  2:4806304/1‑93 (MQ=255)
                                                gAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTg                                         <  1:6324762/93‑1 (MQ=255)
                                                  aCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGCGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                       >  2:31409/1‑93 (MQ=255)
                                                     cTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcgcg                                    <  2:149927/93‑1 (MQ=255)
                                                       ttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGt                                  >  1:3305925/1‑93 (MQ=255)
                                                        tggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTa                                 <  2:3521575/93‑1 (MQ=255)
                                                         ggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTAc                                >  2:5034449/1‑93 (MQ=255)
                                                         ggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTAc                                <  1:4713879/93‑1 (MQ=255)
                                                          gtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACa                               <  1:949199/93‑1 (MQ=255)
                                                            ggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGc                             <  2:1582081/93‑1 (MQ=255)
                                                             gtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCt                            <  1:728088/93‑1 (MQ=255)
                                                             gtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCt                            >  1:2530467/1‑93 (MQ=255)
                                                               ggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtg                          >  2:431717/1‑93 (MQ=255)
                                                                gCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgt                         <  2:5230268/93‑1 (MQ=255)
                                                                 catcatCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgtg                        >  2:1173819/1‑93 (MQ=255)
                                                                  atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       >  2:1855360/1‑93 (MQ=255)
                                                                     atcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgctgc                    >  1:2680214/1‑93 (MQ=255)
                                                                      tcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCg                   >  1:6227241/1‑93 (MQ=255)
                                                                        aGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAc                 <  2:39198/93‑1 (MQ=255)
                                                                         gTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAcc                >  2:4115973/1‑93 (MQ=255)
                                                                              gtagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag           <  1:3588784/93‑1 (MQ=255)
                                                                               tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg          >  1:5171495/1‑93 (MQ=255)
                                                                               tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg          >  2:3588686/1‑93 (MQ=255)
                                                                                agtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGc         <  1:436565/93‑1 (MQ=255)
                                                                                 gtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCg        <  2:3398124/93‑1 (MQ=255)
                                                                                  taAGCACGCCCCCAGGTGGAGGCGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       <  1:4065791/93‑1 (MQ=255)
                                                                                  taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       >  2:4602235/1‑93 (MQ=255)
                                                                                  taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       <  1:821767/93‑1 (MQ=255)
                                                                                   aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      <  1:5544331/93‑1 (MQ=255)
                                                                                   aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      <  1:6124496/93‑1 (MQ=255)
                                                                                   aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      >  1:1013326/1‑93 (MQ=255)
                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGTTACAGCTTGTGCTGCTACCAGCAGGCgaga     <  1:3895861/93‑1 (MQ=255)
                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga     >  2:5101321/1‑93 (MQ=255)
                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga     <  2:2615356/93‑1 (MQ=255)
                                                                                     gcacgcacCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGagag    >  1:1400139/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985651‑3985830

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