Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,491,268 T→A 16.1% D228V (GAT→GTT)  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,2680TA16.1% 345.6 / 48.3 130D228V (GAT→GTT) oxcoxalyl CoA decarboxylase, ThDP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (67/42);  new base A (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 = 9.98e-01

Only 100 of 143 total aligned reads displayed.

CCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCA  >  NC_000913/2491170‑2491367
                                                                                                  |                                                                                                    
ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                     >  2:3732276/1‑100 (MQ=255)
ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                     >  2:4095671/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)
    aGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCt                                                                                                 <  1:3984554/100‑1 (MQ=255)
     gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCtt                                                                                                >  1:2448263/1‑100 (MQ=255)
      aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCGCGAAGCTGTTCATCAGCTtg                                                                                              <  1:2533472/100‑1 (MQ=255)
       aaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTtgt                                                                                             <  2:179415/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)
           ggATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAt                                                                                         <  1:408329/100‑1 (MQ=37)
            gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtt                                                                                        <  2:1577730/100‑2 (MQ=37)
             aTGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAttt                                                                                       <  1:1415391/100‑3 (MQ=25)
              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)
                cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGATGAACAGCTTCGTGAAtttat                                                                                     <  1:2444899/100‑5 (MQ=37)
                cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttatt                                                                                    <  2:801221/100‑6 (MQ=25)
                cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttatt                                                                                    <  2:3990671/100‑6 (MQ=25)
                 gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttattg                                                                                   <  1:2714967/100‑7 (MQ=21)
                 gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTCCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttattg                                                                                   <  2:2874243/100‑7 (MQ=21)
                    tCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttattgaaa                                                                                <  1:2622930/100‑10 (MQ=18)
                     cTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGcc                                                                               >  2:2453859/1‑100 (MQ=255)
                     cTTCAAGGATCCCTTACGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATTAATTCACAAGCTGATGAACAGCTTCGTGAAtttattgaaa                                                                                <  2:1713885/99‑10 (MQ=18)
                      ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCt                                                                              >  2:541231/1‑100 (MQ=255)
                        cAAGGAGCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCtt                                                                             >  2:3973574/1‑99 (MQ=255)
                         aaGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTg                                                                           <  2:726794/100‑1 (MQ=255)
                         aaGGATCCCTTTCGCCAGAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTg                                                                           >  1:1373640/1‑100 (MQ=255)
                            gATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCa                                                                        >  1:2073380/1‑100 (MQ=255)
                             atgccgttCGCCATATACATTGGCAGGAATGGAATCTGAGCACTTTCAAGAAATTCACGAGGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCa                                                                        <  2:3799449/93‑1 (MQ=255)
                              tCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAg                                                                      >  1:3674207/1‑100 (MQ=255)
                               cccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAgg                                                                     <  1:3775364/100‑1 (MQ=255)
                               cccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAgg                                                                     <  1:2193866/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)
                                      gCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaatta                                                              >  2:3278253/1‑100 (MQ=255)
                                       ccATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaa                                                             >  1:1021715/1‑100 (MQ=255)
                                       ccATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaa                                                             >  1:504771/1‑100 (MQ=255)
                                        cATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                            >  1:3064289/1‑100 (MQ=255)
                                        cATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                            >  1:2425887/1‑100 (MQ=255)
                                           aGACATTGGCAGGAATGTAATCTGAGCACTTACAATAAATTGACGAAGCTGTTCATCAGCGTGTGAATGCGCCGCGCCTTTGCCAAGGAtaat                                                                <  2:2556993/93‑1 (MQ=255)
                                             aCATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCACCTTTGCCAAGGATAATTAATGGCCg                                                       >  2:2535055/1‑100 (MQ=255)
                                                ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTc                                                    >  2:1427535/1‑100 (MQ=255)
                                                ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTc                                                    >  1:2145373/1‑100 (MQ=255)
                                                    cAGGAATGGAATCTGAGCACTTTCAATAGATTCACGAAGCTGTTCATCGGGGTGTGATTACTCCGCGCTTTTGTTGAGGGTAATTAATGTCCGtgtag                                                  >  1:3875978/1‑94 (MQ=255)
                                                    cAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCt                                                <  2:248727/100‑1 (MQ=255)
                                                      ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                              >  1:1644120/1‑100 (MQ=255)
                                                      ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                              >  1:1101542/1‑100 (MQ=255)
                                                      ggaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                              <  2:2703205/100‑1 (MQ=255)
                                                       gaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                             >  1:1151379/1‑100 (MQ=255)
                                                       gaatggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                             <  2:2736069/100‑1 (MQ=255)
                                                         atggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGc                                           >  1:524511/1‑100 (MQ=255)
                                                         atggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTGGc                                           >  1:2497496/1‑100 (MQ=255)
                                                           ggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTa                                         >  1:3469909/1‑100 (MQ=255)
                                                           gcactCTGAGCACTCTCAATAAATTCACGAAGCTAATCATCAGCTTCTGAATACGCCCCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCGTTTTCt                                          >  2:1984621/5‑99 (MQ=255)
                                                             aatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaa                                        >  2:857227/1‑99 (MQ=255)
                                                              atCTGAGCATTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                      <  2:449595/100‑1 (MQ=255)
                                                              atCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                      <  1:1815313/100‑1 (MQ=255)
                                                               tCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAg                                     <  1:3572791/100‑1 (MQ=255)
                                                                 tGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGCTAATTAATGGCCGTTCAGCTTTTGCTAAAAGCg                                   >  2:295012/1‑100 (MQ=255)
                                                                 tGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCg                                   >  2:3798642/1‑100 (MQ=255)
                                                                  gAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                  >  1:1298322/1‑100 (MQ=255)
                                                                  gAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                  >  2:2084338/1‑100 (MQ=255)
                                                                   aGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGaa                                 >  1:2017644/1‑100 (MQ=255)
                                                                     cACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAAt                               <  1:2584658/100‑1 (MQ=255)
                                                                      aCTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAAtt                              <  2:2913531/100‑1 (MQ=255)
                                                                       cTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTg                             >  1:3512158/1‑100 (MQ=255)
                                                                       cTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTg                             >  1:3331018/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)
                                                                            aaTAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTa                        >  1:1647974/1‑100 (MQ=255)
                                                                               aaaTTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtg                     >  1:1996635/1‑100 (MQ=255)
                                                                               aaaTTCACGAAGCTGCTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtg                     <  2:3324445/100‑1 (MQ=255)
                                                                                 aTTCACAAGCTGATGAACAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact                  >  1:1520043/1‑100 (MQ=255)
                                                                                  ttCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCCAGtgact                  <  1:2176535/100‑1 (MQ=255)
                                                                                   tCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg                 <  2:2073367/100‑1 (MQ=255)
                                                                                    tcacaaGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact                  >  2:2623645/5‑99 (MQ=37)
                                                                                    cACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactga                >  1:3308578/1‑100 (MQ=255)
                                                                                    cACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactga                >  1:138909/1‑100 (MQ=255)
                                                                                      cGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTgactgact              >  1:370791/1‑100 (MQ=255)
                                                                                      cGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTgactgact              >  1:505589/1‑100 (MQ=255)
                                                                                         aGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCg           <  1:3701495/100‑1 (MQ=255)
                                                                                          gCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg          <  1:578573/100‑1 (MQ=255)
                                                                                             gTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCa       >  2:1497594/1‑100 (MQ=255)
                                                                                              ttCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCAt      >  2:2148949/1‑100 (MQ=255)
                                                                                               tcatcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATg     >  2:4085646/1‑100 (MQ=255)
                                                                                                catcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATgg    >  1:3642014/1‑100 (MQ=255)
                                                                                                 atcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc   <  1:2126670/100‑1 (MQ=255)
                                                                                                 atcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc   >  1:3884381/1‑100 (MQ=255)
                                                                                                 atcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc   <  2:504767/100‑1 (MQ=255)
                                                                                                 atcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc   <  1:1674233/100‑1 (MQ=255)
                                                                                                 atcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCAAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATgg    <  1:683121/99‑1 (MQ=255)
                                                                                                  |                                                                                                    
CCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCA  >  NC_000913/2491170‑2491367

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: