Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,491,264 T→A 16.5% E229D (GAA→GAT) ‡ 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,2640TA16.5% 350.1 / 49.9 135E229D (GAA→GAT) oxcoxalyl CoA decarboxylase, ThDP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (69/42);  new base A (9/13);  total (79/56)
Fisher's exact test for biased strand distribution p-value = 9.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 140 total aligned reads displayed.

AGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCAT  >  NC_000913/2491165‑2491363
                                                                                                   |                                                                                                    
aGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGt                                                                                                      <  1:3791098/100‑1 (MQ=255)
    gCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtca                                                                                                   >  2:2860142/1‑99 (MQ=255)
     ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                 >  2:3732276/1‑100 (MQ=255)
     ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                 >  2:4095671/1‑100 (MQ=255)
      cGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatca                                                                                                <  2:722840/100‑1 (MQ=255)
       gCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGttagatcg                                                                                               >  1:234235/1‑94 (MQ=255)
       gCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAg                                                                                               >  2:4019794/1‑100 (MQ=255)
       gCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAg                                                                                               >  2:1685510/1‑100 (MQ=255)
        cAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGc                                                                                              >  2:2746745/1‑100 (MQ=255)
         aGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCt                                                                                             <  1:3984554/100‑1 (MQ=255)
          gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCtt                                                                                            >  1:2448263/1‑100 (MQ=255)
           aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑g                                                                                          <  1:2532930/100‑1 (MQ=255)
            aaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gt                                                                                         <  2:179415/100‑1 (MQ=255)
             aGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gtg                                                                                        <  1:2746761/100‑1 (MQ=255)
             aGTGGATGCGTATCTTCAAGGATCCATTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gtg                                                                                        <  2:53590/100‑1 (MQ=255)
                ggATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAAt                                                                                     >  1:3998259/1‑100 (MQ=255)
                ggATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAt                                                                                     <  1:408329/100‑1 (MQ=37)
                   tGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcg                                                                                  >  1:2587765/1‑100 (MQ=255)
                   tGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAc                                                                                   >  2:445345/1‑99 (MQ=255)
                    gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgc                                                                                 <  2:3536086/100‑1 (MQ=255)
                     cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgcc                                                                                >  1:2951029/1‑100 (MQ=255)
                     cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgcc                                                                                >  2:2153990/1‑100 (MQ=255)
                     cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgcc                                                                                <  1:996243/100‑1 (MQ=255)
                     cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgc                                                                                 <  2:42613/99‑1 (MQ=255)
                     cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGATGAACAGCTTCGTGAAtttat                                                                                 <  1:2444899/100‑5 (MQ=37)
                     cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttatt                                                                                <  2:801221/100‑6 (MQ=25)
                     cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttatt                                                                                <  2:3990671/100‑6 (MQ=25)
                      gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttattg                                                                               <  1:2714967/100‑7 (MQ=21)
                      gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTCCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttattg                                                                               <  2:2874243/100‑7 (MQ=21)
                       tATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTC‑CCAAGCTGATGAACAGCTTCGTGAAtttattg                                                                               <  2:1680776/99‑7 (MQ=25)
                          cTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGcc                                                                           >  2:2453859/1‑100 (MQ=255)
                            tCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCtt                                                                         <  1:2124614/100‑1 (MQ=255)
                             cAAGGAGCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCtt                                                                         >  2:3973574/1‑99 (MQ=255)
                              aaGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTg                                                                       <  2:726794/100‑1 (MQ=255)
                              aaGGATCCCTTTCGCCAGAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTg                                                                       >  1:1373640/1‑100 (MQ=255)
                                ggATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGcc                                                                     >  1:323855/1‑100 (MQ=255)
                                 gATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCa                                                                    >  1:2073380/1‑100 (MQ=255)
                                  atgccgttCGCCATATACATTGGCAGGAATGGAATCTGAGCACTTTCAAGAAATTCACGAGGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCa                                                                    <  2:3799449/93‑1 (MQ=255)
                                   tCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAg                                                                  >  1:3674207/1‑100 (MQ=255)
                                    cccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAgg                                                                 >  2:2534491/1‑100 (MQ=255)
                                    cccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAgg                                                                 <  1:2193866/100‑1 (MQ=255)
                                     ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGa                                                                >  2:3397341/1‑100 (MQ=255)
                                     ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGa                                                                <  1:1669331/100‑1 (MQ=255)
                                      cTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAt                                                               >  1:8612/1‑100 (MQ=255)
                                          cGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaatt                                                           >  2:1815300/1‑100 (MQ=255)
                                          cGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaatt                                                           <  1:1391578/100‑1 (MQ=255)
                                           gCCATAGACATTGGCAGGAATGGAATCTTAGCACGTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaatta                                                          <  2:1219657/100‑1 (MQ=255)
                                           gCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaatta                                                          >  2:1727045/1‑100 (MQ=255)
                                            ccATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaattaa                                                         >  1:1021715/1‑100 (MQ=255)
                                               tAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGCTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATg                                                       >  2:2126212/1‑99 (MQ=255)
                                                aGACATTGGCAGGAATGTAATCTGAGCACTTACAATAAATTGACGAAGCTGTTCATCAGCGT‑GTGAATGCGCCGCGCCTTTGCCAAGGAtaat                                                            <  2:2556993/93‑1 (MQ=255)
                                                 gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                    >  1:2322293/1‑100 (MQ=255)
                                                 gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                    >  1:1399239/1‑100 (MQ=255)
                                                  aCATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCACCTTTGCCAAGGATAATTAATGGCCg                                                   >  2:2535055/1‑100 (MQ=255)
                                                     ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTc                                                >  2:1427535/1‑100 (MQ=255)
                                                       ggCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAg                                              <  1:1949023/100‑1 (MQ=255)
                                                       ggCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAg                                              >  2:2382704/1‑100 (MQ=255)
                                                         cAGGAATGGAATCTGAGCACTTTCAATAGATTCACGAAGCTGTTCATCGGGGT‑GTGATTACTCCGCGCTTTTGTTGAGGGTAATTAATGTCCGtgtag                                              >  1:3875978/1‑94 (MQ=255)
                                                           ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                          >  1:1644120/1‑100 (MQ=255)
                                                           ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                          <  2:2703205/100‑1 (MQ=255)
                                                           ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                          >  1:1101542/1‑100 (MQ=255)
                                                            gaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                         >  1:1151379/1‑100 (MQ=255)
                                                              atggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGc                                       >  1:524511/1‑100 (MQ=255)
                                                              atggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTGGc                                       >  1:2497496/1‑100 (MQ=255)
                                                                ggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTa                                     >  1:3469094/1‑100 (MQ=255)
                                                                ggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTa                                     >  1:3469909/1‑100 (MQ=255)
                                                                gcactCTGAGCACTCTCAATAAATTCACGAAGCTAATCATCAGCTT‑CTGAATACGCCCCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCGTTTTCt                                      >  2:1984621/5‑99 (MQ=255)
                                                                  aatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaa                                    >  2:857227/1‑99 (MQ=255)
                                                                   atCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                  <  1:1815313/100‑1 (MQ=255)
                                                                      tGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGCTAATTAATGGCCGTTCAGCTTTTGCTAAAAGCg                               >  2:295012/1‑100 (MQ=255)
                                                                      tGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCg                               >  2:3798642/1‑100 (MQ=255)
                                                                       gAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                              >  1:1298322/1‑100 (MQ=255)
                                                                        aGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGaa                             >  2:51848/1‑100 (MQ=255)
                                                                        aGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGaa                             >  1:2017644/1‑100 (MQ=255)
                                                                          cACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAAt                           <  1:2584658/100‑1 (MQ=255)
                                                                            cTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTg                         >  1:3512158/1‑100 (MQ=255)
                                                                             tttCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                        >  1:248510/1‑100 (MQ=255)
                                                                               tCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgcgc                      <  1:1471475/100‑1 (MQ=255)
                                                                               tCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgcgc                      <  2:3693353/100‑1 (MQ=255)
                                                                                cggcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTCAAAGCGAAATTGCGCt                     >  2:620839/7‑100 (MQ=21)
                                                                                cggcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCt                     >  2:2681360/7‑100 (MQ=21)
                                                                                 aaTAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTa                    >  1:3752923/1‑100 (MQ=255)
                                                                                 aaTAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTa                    >  1:1647974/1‑100 (MQ=255)
                                                                                  gcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAg                   >  1:216007/5‑100 (MQ=25)
                                                                                   cgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCTGAGCCTTTGCCAAGGATAATTAATGGACTTTCAGCCTTTGCTAAAGGCGAAATTACGCTAGt                  >  2:2934940/4‑100 (MQ=12)
                                                                                    aaaTTCACGAAGCTGCTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtg                 <  2:3324445/100‑1 (MQ=255)
                                                                                      aTTCAC‑AAGCTGATGAACAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact              >  1:1520043/1‑100 (MQ=255)
                                                                                       ttCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact              >  1:5571/1‑100 (MQ=255)
                                                                                       ttCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCCAGtgact              <  1:2176535/100‑1 (MQ=255)
                                                                                        tCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg             <  2:2073367/100‑1 (MQ=255)
                                                                                         cACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactga            >  1:3308578/1‑100 (MQ=255)
                                                                                         cACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactga            >  1:138909/1‑100 (MQ=255)
                                                                                          aCGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactgac           <  1:456675/100‑1 (MQ=255)
                                                                                           cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTgactgact          >  1:505589/1‑100 (MQ=255)
                                                                                           cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTgactgact          >  1:370791/1‑100 (MQ=255)
                                                                                              aGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCg       <  1:3701495/100‑1 (MQ=255)
                                                                                                cTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCggg     <  1:3603743/100‑1 (MQ=255)
                                                                                                 tGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGc    >  1:4069309/1‑100 (MQ=255)
                                                                                                  gTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCa   >  2:1497594/1‑100 (MQ=255)
                                                                                                   ttCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCAt  >  2:2148949/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
AGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCAT  >  NC_000913/2491165‑2491363

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: ‑