Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,491,266 C→G 15.9% E229Q (GAA→CAA) ‡ oxc ← oxalyl CoA decarboxylase, ThDP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,491,2660CG15.9% 362.0 / 53.3 130E229Q (GAA→CAA) oxcoxalyl CoA decarboxylase, ThDP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (67/42);  new base G (8/13);  total (75/55)
Fisher's exact test for biased strand distribution p-value = 5.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 138 total aligned reads displayed.

CTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGG  >  NC_000913/2491167‑2491365
                                                                                                   |                                                                                                    
cTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtc                                                                                                      >  1:2179332/1‑100 (MQ=255)
  gCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtca                                                                                                     >  2:2860142/1‑99 (MQ=255)
   ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                   >  2:4095671/1‑100 (MQ=255)
   ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                   >  2:3732276/1‑100 (MQ=255)
    cGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatca                                                                                                  >  1:4082185/1‑100 (MQ=255)
    cGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatca                                                                                                  <  2:722840/100‑1 (MQ=255)
     gCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAg                                                                                                 >  2:1685510/1‑100 (MQ=255)
     gCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAg                                                                                                 >  2:4019794/1‑100 (MQ=255)
      cAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGc                                                                                                >  2:2746745/1‑100 (MQ=255)
       aGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCAATTTCAATAAATTCACGAAGCTGTTCATCAGCt                                                                                               >  1:1381017/1‑100 (MQ=255)
        gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCtt                                                                                              >  1:2448263/1‑100 (MQ=255)
         aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTtg                                                                                            <  1:2532930/100‑1 (MQ=255)
           aGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTtgtg                                                                                          <  1:2746761/100‑1 (MQ=255)
           aGTGGATGCGTATCTTCAAGGATCCATTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTtgtg                                                                                          <  2:53590/100‑1 (MQ=255)
            gTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGa                                                                                         <  2:1031445/100‑1 (MQ=255)
              ggATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAAt                                                                                       >  1:3998259/1‑100 (MQ=255)
               gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtt                                                                                      <  2:1577730/100‑2 (MQ=37)
                 tGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAcg                                                                                    >  1:2587765/1‑100 (MQ=255)
                 tGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAc                                                                                     >  2:445345/1‑99 (MQ=255)
                  gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAcgc                                                                                   <  2:3536086/100‑1 (MQ=255)
                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAcgcc                                                                                  <  1:996243/100‑1 (MQ=255)
                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAcgcc                                                                                  >  2:2153990/1‑100 (MQ=255)
                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAcgcc                                                                                  >  1:2951029/1‑100 (MQ=255)
                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAcgc                                                                                   <  2:42613/99‑1 (MQ=255)
                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGATGAACAGCTTCGTGAAtttat                                                                                   <  1:2444899/100‑5 (MQ=37)
                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttatt                                                                                  <  2:3990671/100‑6 (MQ=25)
                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttatt                                                                                  <  2:801221/100‑6 (MQ=25)
                    gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttattg                                                                                 <  1:2714967/100‑7 (MQ=21)
                    gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTCCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttattg                                                                                 <  2:2874243/100‑7 (MQ=21)
                        cTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGcc                                                                             >  2:2453859/1‑100 (MQ=255)
                        cTTCAAGGATCCCTTACGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATTAATTCACAAGCTGATGAACAGCTTCGTGAAtttattgaaa                                                                              <  2:1713885/99‑10 (MQ=18)
                          tCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCtt                                                                           <  1:2124614/100‑1 (MQ=255)
                           cAAGGAGCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCtt                                                                           >  2:3973574/1‑99 (MQ=255)
                            aaGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTg                                                                         <  2:726794/100‑1 (MQ=255)
                            aaGGATCCCTTTCGCCAGAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTg                                                                         >  1:1373640/1‑100 (MQ=255)
                              ggATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGcc                                                                       >  1:323855/1‑100 (MQ=255)
                                atgccgttCGCCATATACATTGGCAGGAATGGAATCTGAGCACTTTCAAGAAATTCACGAGGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCa                                                                      <  2:3799449/93‑1 (MQ=255)
                                 tCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAg                                                                    >  1:3674207/1‑100 (MQ=255)
                                  cccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAgg                                                                   <  1:2193866/100‑1 (MQ=255)
                                  cccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAgg                                                                   <  2:3765515/100‑1 (MQ=255)
                                  cccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAgg                                                                   >  2:2534491/1‑100 (MQ=255)
                                   ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGa                                                                  <  1:1669331/100‑1 (MQ=255)
                                   ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGa                                                                  >  2:3397341/1‑100 (MQ=255)
                                        cGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaatt                                                             >  2:1815300/1‑100 (MQ=255)
                                        cGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaatt                                                             <  1:1391578/100‑1 (MQ=255)
                                         gCCATAGACATTGGCAGGAATGGAATCTTAGCACGTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaatta                                                            <  2:1219657/100‑1 (MQ=255)
                                         gCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaatta                                                            >  2:1727045/1‑100 (MQ=255)
                                          ccATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaa                                                           >  1:504771/1‑100 (MQ=255)
                                          ccATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaa                                                           >  1:1021715/1‑100 (MQ=255)
                                             tAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGCTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATg                                                         >  2:2126212/1‑99 (MQ=255)
                                              aGACATTGGCAGGAATGTAATCTGAGCACTTACAATAAATTGACGAAGCTGTTCATCAGCGTGTGAATGCGCCGCGCCTTTGCCAAGGAtaat                                                              <  2:2556993/93‑1 (MQ=255)
                                              aGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGc                                                       <  2:406058/100‑1 (MQ=255)
                                               gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                      >  1:2322293/1‑100 (MQ=255)
                                               gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                      >  1:1399239/1‑100 (MQ=255)
                                                aCATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCACCTTTGCCAAGGATAATTAATGGCCg                                                     >  2:2535055/1‑100 (MQ=255)
                                                   ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTc                                                  >  1:2145373/1‑100 (MQ=255)
                                                     ggCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAg                                                >  2:2382704/1‑100 (MQ=255)
                                                       cAGGAATGGAATCTGAGCACTTTCAATAGATTCACGAAGCTGTTCATCGGGGTGTGATTACTCCGCGCTTTTGTTGAGGGTAATTAATGTCCGtgtag                                                >  1:3875978/1‑94 (MQ=255)
                                                         ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                            >  1:1101542/1‑100 (MQ=255)
                                                         ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                            <  2:2703205/100‑1 (MQ=255)
                                                          gaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                           >  1:1151379/1‑100 (MQ=255)
                                                            atggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGc                                         >  1:524511/1‑100 (MQ=255)
                                                            atggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTGGc                                         >  1:2497496/1‑100 (MQ=255)
                                                              ggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTa                                       <  2:1857574/100‑1 (MQ=255)
                                                              ggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTa                                       >  1:3469909/1‑100 (MQ=255)
                                                              gcactCTGAGCACTCTCAATAAATTCACGAAGCTAATCATCAGCTTCTGAATACGCCCCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCGTTTTCt                                        >  2:1984621/5‑99 (MQ=255)
                                                                aatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaa                                      >  2:857227/1‑99 (MQ=255)
                                                                 atCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                    <  1:1815313/100‑1 (MQ=255)
                                                                    tGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGCTAATTAATGGCCGTTCAGCTTTTGCTAAAAGCg                                 >  2:295012/1‑100 (MQ=255)
                                                                     gAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                >  1:1298322/1‑100 (MQ=255)
                                                                     gAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                >  2:2084338/1‑100 (MQ=255)
                                                                      aGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGaa                               >  2:51848/1‑100 (MQ=255)
                                                                      aGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGaa                               >  1:2017644/1‑100 (MQ=255)
                                                                        cACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAAt                             <  1:2584658/100‑1 (MQ=255)
                                                                          cTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTg                           >  1:3331018/1‑100 (MQ=255)
                                                                          cTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTg                           >  1:3512158/1‑100 (MQ=255)
                                                                           tttCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                          >  1:248510/1‑100 (MQ=255)
                                                                             tCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgcgc                        <  1:1471475/100‑1 (MQ=255)
                                                                             tCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgcgc                        <  2:3693353/100‑1 (MQ=255)
                                                                             gcggcgtaTTCACAAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgcgc                        >  1:3960351/8‑100 (MQ=21)
                                                                              cggcgtaTTCACAAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTCAAAGCGAAATTGCGCt                       >  2:620839/7‑100 (MQ=21)
                                                                              cggcgtaTTCACAAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCt                       >  2:2681360/7‑100 (MQ=21)
                                                                               aaTAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTa                      >  1:1647974/1‑100 (MQ=255)
                                                                               aaTAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTa                      >  1:3752923/1‑100 (MQ=255)
                                                                                gcgtaTTCACAAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAg                     >  1:216007/5‑100 (MQ=25)
                                                                                 cgtaTTCACAAGCTGATGAACAGCTTCGTGAATACGCTGAGCCTTTGCCAAGGATAATTAATGGACTTTCAGCCTTTGCTAAAGGCGAAATTACGCTAGt                    >  2:2934940/4‑100 (MQ=12)
                                                                                 cgtaTTCACAAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGt                    >  1:1562409/4‑100 (MQ=25)
                                                                                  aaaTTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtg                   >  1:1996635/1‑100 (MQ=255)
                                                                                  aaaTTCACGAAGCTGCTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtg                   <  2:3324445/100‑1 (MQ=255)
                                                                                    aTTCACAAGCTGATGAACAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact                >  1:1520043/1‑100 (MQ=255)
                                                                                     ttCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact                >  1:5571/1‑100 (MQ=255)
                                                                                       cACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactga              >  1:3308578/1‑100 (MQ=255)
                                                                                        aCGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactgac             <  1:456675/100‑1 (MQ=255)
                                                                                         cGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTgactgact            >  1:370791/1‑100 (MQ=255)
                                                                                            aGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCg         <  1:3701495/100‑1 (MQ=255)
                                                                                             gCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg        <  1:578573/100‑1 (MQ=255)
                                                                                              cTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCggg       <  1:3603743/100‑1 (MQ=255)
                                                                                               tGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGc      >  1:4069309/1‑100 (MQ=255)
                                                                                                gTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCa     >  2:1497594/1‑100 (MQ=255)
                                                                                                   catcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATgg  >  1:3642014/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
CTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGG  >  NC_000913/2491167‑2491365

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