Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,430 T→C I322I (ATT→ATC cysG → uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,4300TC99.4% 1132.8 / ‑3.9 355I322I (ATT→ATCcysGuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (178/173);  total (180/175)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

Only 100 of 371 total aligned reads displayed.

AGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAG  >  NC_000913/3496340‑3496519
                                                                                          |                                                                                         
aGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATccc                                                                                         >  1:10023554/1‑93 (MQ=255)
aGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATccc                                                                                         >  1:3277629/1‑93 (MQ=255)
aGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATccc                                                                                         <  1:7586766/93‑1 (MQ=255)
 ggCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCg                                                                                        >  1:1107594/1‑93 (MQ=255)
 ggCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCg                                                                                        >  2:6770019/1‑93 (MQ=255)
  gCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGt                                                                                       <  2:662755/93‑1 (MQ=255)
  gCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGt                                                                                       <  1:7412717/93‑1 (MQ=255)
   cAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGtt                                                                                      >  1:3748710/1‑93 (MQ=255)
   cAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGtt                                                                                      <  2:5403947/93‑1 (MQ=255)
     aaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtct                                                                                    >  1:415096/1‑93 (MQ=255)
      aCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtctc                                                                                   <  1:3684299/93‑1 (MQ=255)
       cgcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCg                                                                                  <  1:4630651/93‑1 (MQ=255)
       cgcgTGGTGCCGCTGAAAGGTGGCGATCCTTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCg                                                                                  <  1:2174172/93‑1 (MQ=255)
        gcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCgg                                                                                 >  1:9582134/1‑93 (MQ=255)
        gcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCgg                                                                                 <  2:2098559/93‑1 (MQ=255)
        gcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCgg                                                                                 >  2:6126260/1‑93 (MQ=255)
           tggtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCggtgg                                                                              >  1:3515409/1‑93 (MQ=255)
             gtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGtt                                                                            >  2:186958/1‑93 (MQ=255)
             gtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGtt                                                                            >  1:6119266/1‑93 (MQ=255)
               gCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTcc                                                                          <  2:6305019/93‑1 (MQ=255)
                cGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCg                                                                         >  1:7810674/1‑93 (MQ=255)
                 gcatGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCgg                                                                        <  2:8072506/90‑1 (MQ=255)
                    tGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTa                                                                     >  2:1135135/1‑93 (MQ=255)
                     gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTAt                                                                    <  1:4283465/93‑1 (MQ=255)
                     gAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTAt                                                                    <  1:685261/93‑1 (MQ=255)
                      aaaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTAtt                                                                   >  1:5564824/1‑93 (MQ=255)
                      aaaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTAtt                                                                   <  1:636059/93‑1 (MQ=255)
                       aaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTa                                                                  >  1:4406784/1‑93 (MQ=255)
                       aaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTa                                                                  <  1:5492093/93‑1 (MQ=255)
                         ggTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTAcc                                                                >  2:4667211/1‑93 (MQ=255)
                           tGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGc                                                              >  1:3062908/1‑93 (MQ=255)
                           tGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGc                                                              >  2:4254304/1‑93 (MQ=255)
                           tGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGGGGTTCCGGGGATTACCGc                                                              >  2:3453938/1‑93 (MQ=255)
                            ggCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCa                                                             <  1:7768091/93‑1 (MQ=255)
                             gCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAg                                                            <  2:4394611/93‑1 (MQ=255)
                              cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACTCTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGc                                                           <  1:2350657/93‑1 (MQ=255)
                              cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGc                                                           >  1:8860813/1‑93 (MQ=255)
                              cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGc                                                           >  2:3160326/1‑93 (MQ=255)
                               gATCCGTTTATTTTTGGCCGTGGTAGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCt                                                          >  2:2245867/1‑93 (MQ=255)
                                aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCtt                                                         >  1:698954/1‑93 (MQ=255)
                                aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCtt                                                         >  1:1500009/1‑93 (MQ=255)
                                  ccGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCt                                                       <  1:1936092/93‑1 (MQ=255)
                                  ccGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCt                                                       >  2:9517381/1‑93 (MQ=255)
                                   cGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTg                                                      <  2:7415389/93‑1 (MQ=255)
                                   cGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTg                                                      <  2:6405445/93‑1 (MQ=255)
                                    gTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTgg                                                     >  1:577267/1‑93 (MQ=255)
                                    gTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTgg                                                     >  1:7118012/1‑93 (MQ=255)
                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGt                                                    >  1:813813/1‑93 (MQ=255)
                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGt                                                    >  1:2634526/1‑93 (MQ=255)
                                      ttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGtt                                                   <  2:7665097/93‑1 (MQ=255)
                                      ttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGtt                                                   <  2:1308449/93‑1 (MQ=255)
                                       tATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTg                                                  >  2:5575615/1‑93 (MQ=255)
                                         tttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGct                                                >  1:5051685/1‑93 (MQ=255)
                                          ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGctc                                               >  1:6742467/1‑93 (MQ=255)
                                          ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGctc                                               <  2:7487378/93‑1 (MQ=255)
                                             tGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATACCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGc                                            <  2:3720105/93‑1 (MQ=255)
                                               gCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCt                                          <  2:4592149/93‑1 (MQ=255)
                                                 cGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTAt                                        <  1:8506953/93‑1 (MQ=255)
                                                 cGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTAt                                        <  1:5387640/93‑1 (MQ=255)
                                                  gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                       >  1:3226629/1‑93 (MQ=255)
                                                      tggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                   >  1:82195/1‑93 (MQ=255)
                                                        gCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTa                                 <  2:2448471/93‑1 (MQ=255)
                                                        gCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTa                                 <  2:8722685/93‑1 (MQ=255)
                                                          gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                               >  1:4459196/1‑93 (MQ=255)
                                                          gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                               >  1:8676557/1‑93 (MQ=255)
                                                           aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                              <  2:8150358/93‑1 (MQ=255)
                                                           aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                              >  1:6622473/1‑93 (MQ=255)
                                                           aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                              >  2:5040698/1‑93 (MQ=255)
                                                            agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTTGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             >  2:8426210/1‑93 (MQ=255)
                                                            agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             >  2:6823165/1‑93 (MQ=255)
                                                            agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             >  2:5979298/1‑93 (MQ=255)
                                                             gagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCa                            >  2:9311425/1‑93 (MQ=255)
                                                              agCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAc                           >  1:1685047/1‑93 (MQ=255)
                                                              agCTGGAAACACTGTGCAACGCGCGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAc                           <  2:5502572/93‑1 (MQ=255)
                                                                 tGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCa                        >  1:5266619/1‑93 (MQ=255)
                                                                  ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCAc                       >  1:9514036/1‑93 (MQ=255)
                                                                  ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCAc                       >  1:6224368/1‑93 (MQ=255)
                                                                  ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCAc                       <  1:1738528/93‑1 (MQ=255)
                                                                  ggAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGAATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAATCag                       >  1:9038593/1‑92 (MQ=255)
                                                                   gAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACg                      >  2:7071754/1‑93 (MQ=255)
                                                                    aaaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGc                     >  1:1124313/1‑93 (MQ=255)
                                                                     aaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCa                    <  2:7977470/93‑1 (MQ=255)
                                                                      acacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCAt                   <  1:342244/93‑1 (MQ=255)
                                                                      acacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCAt                   >  2:864578/1‑93 (MQ=255)
                                                                       cacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                  <  2:9751988/93‑1 (MQ=255)
                                                                        acTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcg                 >  2:8283316/1‑93 (MQ=255)
                                                                         cTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgc                >  2:3014538/1‑93 (MQ=255)
                                                                          tgtgCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgcg               >  1:7932835/1‑93 (MQ=255)
                                                                          tgtgCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgcg               >  1:1823059/1‑93 (MQ=255)
                                                                             gcgaCGCGGGTCACCCGATATCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAtt            <  2:8971237/90‑1 (MQ=255)
                                                                              cAACGCGGGTATCCCGTTCTGTGTGGTTCCGGGTATTACCGCGGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa           <  2:272793/93‑1 (MQ=255)
                                                                              cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa           <  2:5328195/93‑1 (MQ=255)
                                                                              cAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa           <  2:4746529/93‑1 (MQ=255)
                                                                               aaCGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTAt          <  2:1686279/93‑1 (MQ=255)
                                                                                aCGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATg         >  2:5198166/1‑93 (MQ=255)
                                                                                 cgcgGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc        <  1:2832663/93‑1 (MQ=255)
                                                                                     ggTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCag    >  2:1471922/1‑93 (MQ=255)
                                                                                     ggTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCag    >  1:7230860/1‑93 (MQ=255)
                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga   <  2:8886950/93‑1 (MQ=255)
                                                                                       tATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCagag  >  1:9248598/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAG  >  NC_000913/3496340‑3496519

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