Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,224,188 C→T R928C (CGT→TGT)  metH → B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,224,1880CT85.6% 743.7 / 78.8 278R928C (CGT→TGT) metHB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand):  ref base C (25/15);  new base T (114/124);  total (139/139)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 278 total aligned reads displayed.

GAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCG  >  NC_000913/4224097‑4224278
                                                                                           |                                                                                          
gAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCg                                                                                           >  2:2087340/1‑93 (MQ=255)
 aaaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGt                                                                                          <  2:2649721/93‑1 (MQ=255)
  aaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGtc                                                                                         <  2:5975578/93‑1 (MQ=255)
  aaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGtc                                                                                         >  1:5457208/1‑93 (MQ=255)
  aaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGtc                                                                                         >  2:1356510/1‑93 (MQ=255)
   aCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGtct                                                                                        <  2:2929186/93‑1 (MQ=255)
   aCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGtct                                                                                        >  2:2567670/1‑93 (MQ=255)
    ccGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGtctc                                                                                       <  2:7454526/93‑1 (MQ=255)
     cgcgCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCg                                                                                      >  2:7334554/1‑93 (MQ=255)
     cgcgCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCg                                                                                      <  1:4592133/93‑1 (MQ=255)
       cgcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGc                                                                                    <  1:587452/93‑1 (MQ=255)
       cgcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGc                                                                                    >  2:5173406/1‑93 (MQ=255)
        gcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCg                                                                                   >  2:5235804/1‑93 (MQ=255)
           caccacCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGc                                                                                >  2:4678623/1‑93 (MQ=255)
           caccacCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGc                                                                                >  2:3479571/1‑93 (MQ=255)
            accaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCa                                                                               <  1:7428808/93‑1 (MQ=255)
             ccaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAg                                                                              <  1:6167622/93‑1 (MQ=255)
              caccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAgg                                                                             >  2:6647492/1‑93 (MQ=255)
              caccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAgg                                                                             <  1:4888548/93‑1 (MQ=255)
               accGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGa                                                                            >  2:589991/1‑93 (MQ=255)
                 cGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAg                                                                          >  1:745238/1‑93 (MQ=255)
                 cGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAg                                                                          <  2:2406167/93‑1 (MQ=255)
                  ggTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGt                                                                         >  2:4211827/1‑93 (MQ=255)
                   gTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTc                                                                        <  1:3300289/93‑1 (MQ=255)
                    tCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCg                                                                       >  2:1796303/1‑93 (MQ=255)
                     cACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGa                                                                      <  1:5050244/93‑1 (MQ=255)
                        gCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGc                                                                   <  1:6790933/93‑1 (MQ=255)
                        gCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGc                                                                   >  2:3989372/1‑93 (MQ=255)
                        gCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGc                                                                   >  1:6678493/1‑93 (MQ=255)
                        gCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGc                                                                   >  1:3009801/1‑93 (MQ=255)
                        gCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCCGTGGCGCACCGTCTCCGCGCGCAGGAAGTCGAAGc                                                                   >  1:5275223/1‑93 (MQ=255)
                         cTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGGCGAAGcc                                                                  >  1:4563147/1‑93 (MQ=255)
                          tGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCa                                                                 <  1:5642949/93‑1 (MQ=255)
                          tGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCa                                                                 >  1:6319302/1‑93 (MQ=255)
                             aaGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCa                                                              <  2:7191068/93‑1 (MQ=255)
                               gcggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCGGGAAGTCGAAGCCAGCATc                                                            >  1:7523869/1‑93 (MQ=255)
                               gcggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATc                                                            <  1:431195/93‑1 (MQ=255)
                               gcggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATc                                                            >  2:6519795/1‑93 (MQ=255)
                                cggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCg                                                           <  2:4732720/93‑1 (MQ=255)
                                 ggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGa                                                          >  1:3504329/1‑93 (MQ=255)
                                 ggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGa                                                          <  1:5346556/93‑1 (MQ=255)
                                   cgcgcgATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGaaa                                                        <  2:4999215/93‑1 (MQ=255)
                                   cgcgcgATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGaaa                                                        <  1:357370/93‑1 (MQ=255)
                                   cgcgcgATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGaaa                                                        <  2:4412825/93‑1 (MQ=255)
                                    gcgcgATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAAc                                                       <  2:1037579/93‑1 (MQ=255)
                                    gcgcgATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAAc                                                       >  2:1509797/1‑93 (MQ=255)
                                        gATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTg                                                   >  1:7486965/1‑93 (MQ=255)
                                        gATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTg                                                   >  1:2627310/1‑93 (MQ=255)
                                         aTAACGATTTCGCTTTTGACTGGTAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGc                                                  <  2:6263632/93‑1 (MQ=255)
                                          tAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCg                                                 >  2:6960480/1‑93 (MQ=255)
                                          tAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCg                                                 >  1:7159734/1‑93 (MQ=255)
                                          tAACGATTTCGCTTTTGACTGGCAGGCTTACACCCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGGCGAAGCCAGCATCGAAACGCTGCg                                                 >  1:2028314/1‑93 (MQ=255)
                                           aaCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGt                                                >  2:4041969/1‑93 (MQ=255)
                                            aCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTa                                               >  1:1318963/1‑93 (MQ=255)
                                            aCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTa                                               <  2:5550211/93‑1 (MQ=255)
                                            aCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTa                                               <  1:949661/93‑1 (MQ=255)
                                             cGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTaa                                              <  1:7068849/93‑1 (MQ=255)
                                              gATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAAt                                             >  2:3916026/1‑93 (MQ=255)
                                              gATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAAt                                             >  1:3880455/1‑93 (MQ=255)
                                                tttCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTa                                           >  2:7034456/1‑93 (MQ=255)
                                                  tCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACa                                         >  1:1603510/1‑93 (MQ=255)
                                                    gCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATc                                       >  1:288359/1‑93 (MQ=255)
                                                    gCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATc                                       >  2:1643433/1‑93 (MQ=255)
                                                    gCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATc                                       <  2:6158143/93‑1 (MQ=255)
                                                    gCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATc                                       <  2:3177292/93‑1 (MQ=255)
                                                      ttttGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGa                                     >  1:6379943/1‑93 (MQ=255)
                                                      ttttGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGa                                     <  1:4843846/93‑1 (MQ=255)
                                                      ttttGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGa                                     >  2:191187/1‑93 (MQ=255)
                                                      ttttGACTGGCAGGCTTACACGCCGCCGGGGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGGAATTACATCGa                                     >  1:2750008/1‑93 (MQ=255)
                                                        ttGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACt                                   <  2:452259/93‑1 (MQ=255)
                                                         tGACTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTg                                  >  1:7306206/1‑93 (MQ=255)
                                                            cTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGac                               >  1:5690607/1‑93 (MQ=255)
                                                            cTGGCAGGCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGac                               >  1:3738437/1‑93 (MQ=255)
                                                                  ggCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGtt                         >  1:2111375/1‑93 (MQ=255)
                                                                  ggCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGtt                         >  1:6878037/1‑93 (MQ=255)
                                                                  ggCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGtt                         >  2:1200252/1‑93 (MQ=255)
                                                                   gCTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTc                        <  1:6430405/93‑1 (MQ=255)
                                                                    cTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCt                       >  2:4272221/1‑93 (MQ=255)
                                                                    cTTACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCt                       <  2:3782430/93‑1 (MQ=255)
                                                                     ttACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCtt                      >  1:2324787/1‑93 (MQ=255)
                                                                     ttACACGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCtt                      <  1:1109512/93‑1 (MQ=255)
                                                                       acacGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTa                    <  1:5232703/93‑1 (MQ=255)
                                                                       acacGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTa                    <  2:7111427/93‑1 (MQ=255)
                                                                       acacGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTa                    <  1:2712460/93‑1 (MQ=255)
                                                                        cacGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTAt                   >  1:5998744/1‑93 (MQ=255)
                                                                         acGCCGCCGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATg                  >  2:5811093/1‑93 (MQ=255)
                                                                           gccgccGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGAc                >  1:6985809/1‑93 (MQ=255)
                                                                           gccgccGGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGAc                >  1:4120767/1‑93 (MQ=255)
                                                                            ccgccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGAcc               >  1:3486039/1‑93 (MQ=255)
                                                                             cgccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCt              <  1:4705438/93‑1 (MQ=255)
                                                                             cgccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCt              <  1:5114621/93‑1 (MQ=255)
                                                                              gccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTg             <  1:1917978/93‑1 (MQ=255)
                                                                              gccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTg             <  2:3580549/93‑1 (MQ=255)
                                                                              gccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTg             >  1:3098937/1‑93 (MQ=255)
                                                                               ccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTgg            >  2:6020325/1‑93 (MQ=255)
                                                                               ccgGTGGCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTgg            >  1:6932972/1‑93 (MQ=255)
                                                                                    ggCGCACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCt       <  1:5580245/93‑1 (MQ=255)
                                                                                       gcACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGc    <  1:2475191/93‑1 (MQ=255)
                                                                                        cACTGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGcc   >  2:4945972/1‑93 (MQ=255)
                                                                                         aCCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCg  >  2:4561822/1‑93 (MQ=255)
                                                                                           |                                                                                          
GAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCG  >  NC_000913/4224097‑4224278

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