Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,135,607 A→C 26.5% Q32P (CAG→CCG)  tgnC → (Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,135,6070AC26.5% 64.7 / 8.6 72Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base C (16/5);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 3.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  KT2440/4135478‑4135741
                                                                                                                                 |                                                                                                                                      
cgaccgcTCACCGAATCCGCCACGGGGATTGTTGATTGACCGTCGTTCGTTTCCAGTTGTGCATCGTCGGCCAATGGCGCGATGCCCAGAGTGGCAAGACCTTCGAAAGCCTCAACCCGTCCACCTCCCAGGCCt                                                                                                                                   <  8:357915/129‑1 (MQ=255)
cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCCCCGCCCCGGccg                                                                                                                                   >  1:119112/1‑134 (MQ=255)
cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGccg                                                                                                                                   >  7:122408/1‑134 (MQ=255)
cGCACCCTCACCGAATCCGACAAGGTGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCAACGACGGCCAATGGCGCGATGCCCAGAGCGGCACCCCATTCGACAGACTCACCCCGGCCACCGCCCCGGcca                                                                                                                                   >  1:337754/1‑134 (MQ=255)
cGCACCCTCACCGAATACGCCAAGGGGATTGTTGCATGACCCTCTTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCAGccca                                                                                                                                   >  6:236111/1‑134 (MQ=255)
  cACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGATCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGACTgg                                                                                                                                 >  1:54885/1‑135 (MQ=255)
   aCCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGCCCTTCGACAGCCTCAACCCGGCCACCGCCCAGcccgcgg                                                                                                                                >  3:323596/1‑131 (MQ=255)
     ccTCACCGAATCCGCCAAGGGTTTTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGgcg                                                                                                                              <  8:7222/135‑1 (MQ=255)
     ccTCACCGAATCCGCCAAGGGGATTGTTGCATGCCCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACCGCCTCAACCCGGCCACCGCCCAGGCCTGGgcg                                                                                                                              >  6:189634/1‑135 (MQ=255)
     ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCAGGgcg                                                                                                                              >  2:184824/1‑135 (MQ=255)
      cTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGCCCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGgcgc                                                                                                                             >  3:238297/1‑135 (MQ=255)
           cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCCCCCCGGCCGGGGCGCCACTg                                                                                                                        >  2:11432/1‑135 (MQ=255)
           cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCACCTg                                                                                                                        >  5:270498/1‑135 (MQ=255)
            gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGgcgccccggc                                                                                                                       >  3:348121/1‑129 (MQ=255)
               tCCGCCAGGGGGATGGTTGCATGACCGTCGTTCGTTTCCATAGGTGCATCGACGGCCAATGCCGCGATGCCCAGAGCGGCAAGACCTGCGCCCGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCg                                                                                                                    <  5:316003/135‑1 (MQ=255)
                  gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCCACTGCCCGACg                                                                                                                 >  2:228066/1‑135 (MQ=255)
                  gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGACg                                                                                                                 >  2:269867/1‑135 (MQ=255)
                  gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCccgccg                                                                                                                 >  6:94380/1‑135 (MQ=255)
                  gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACg                                                                                                                 >  4:110808/1‑135 (MQ=255)
                    cAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGcc                                                                                                               >  4:202572/1‑135 (MQ=255)
                         ggATTGTTGCATCACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGAGGCCCAGAGCGCCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                          <  3:77294/135‑1 (MQ=255)
                           aTTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAg                                                                                                        <  4:252537/135‑1 (MQ=255)
                               tggCATGACCCTGGTTGGTTTCCAGATTGGCATCGAGGGCCAATGGGGCTATGCCCAGAGCGGAAATACCTTCTACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCg                                                                                                     <  4:35561/132‑1 (MQ=255)
                                tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGAt                                                                                                   >  8:118469/1‑135 (MQ=255)
                                      aCCCTGGTTCGTTCCCAGAGGTCCATCGCCGCCCAAGGGGGCGAGGCCCAGAGCGGCAAGACCTCCGCCAGCCCCACCCCGCCCCCCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAg                                                                                             <  7:308873/135‑1 (MQ=255)
                                       cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCATGCCTTGGCGCAACAGCCCCACGCCGATGAAGCCGAGTTGGAga                                                                                            >  2:66874/1‑134 (MQ=255)
                                       cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGGGCAACTGCCCGCCGCCGCTCAAGCCGCTGTGGGGc                                                                                            >  6:263171/1‑135 (MQ=255)
                                       cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGc                                                                                            <  6:80500/135‑1 (MQ=255)
                                             ttcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCTACTGCCCCACCCCGATCAAGCCGATGTGGAGCTGGCCg                                                                                      >  5:287575/1‑135 (MQ=255)
                                             ttcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCACGCCAGGGCGCCACAGCCCCCCCCCGCTGAAGCCCATGTGGCGCGGGCCg                                                                                      >  6:167155/1‑135 (MQ=255)
                                              tcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCAACTGCCCGCCGCCGCAGCAGCCGAAGTGGCGCGGGGCGt                                                                                     >  5:87941/1‑135 (MQ=255)
                                               cgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGCCAAGGCCCGCCGCCGATGACGCCGATGTGGAGCTGGccct                                                                                     >  6:44534/1‑132 (MQ=255)
                                                gttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTc                                                                                   <  4:15538/135‑1 (MQ=255)
                                                    ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCCATGAAGCCGATGTGGAGCTGGCCGTTCAGgc                                                                               >  8:308873/1‑135 (MQ=255)
                                                     cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACATCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGgcc                                                                              >  8:358052/1‑135 (MQ=255)
                                                         tgggCAACGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGccc                                                                          <  4:50080/132‑1 (MQ=255)
                                                         ggggCATCGACGGCCATGGGCGCGAGGCCCAGCGCGGCAACACCTGCGACAGCCTCAACCCGCCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGccc                                                                          <  4:187559/132‑1 (MQ=255)
                                                           ttcATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGACTGGGCGCAACAGCCCGACGCCCATGAAGGCGAAGTGGCGCGGGCCGTTCAGGACCCCCAg                                                                        >  2:304536/3‑135 (MQ=255)
                                                            gaagccaCGGCCAAGGGCGCGATGCCCAGAGCGGCAACACCTTCGACCGCCGCACCCGGGCCCCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGc                                                                       <  8:342288/129‑1 (MQ=255)
                                                            gCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTGGACAGCCTCAACCCGGCCACCGCCCAGGCGTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGc                                                                       <  4:293953/135‑1 (MQ=255)
                                                             cATAGACGGCCAATGGCGCGATCCCCAGAGCGGCAAGACCTTCGACAGCCTCAACACGCCCACCCCCCAGGCAGGGGCGCAAAAGGCCGCCACCCCAGAAGACGCAGATGAGCTGGCCGTTAAGGCCGCCCaccc                                                                      >  6:349291/1‑132 (MQ=255)
                                                               tCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCACCTGCCCGACGCCGCTGAAGCCGCTGTGGCGCTGGCCGTTGAGGCCGCCCCGCGTg                                                                    >  8:200417/1‑135 (MQ=255)
                                                                cGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCACTGCCCGCCGCCGCGGAAGCCGCTGTGGGGCTGGCCGGTGAGGCCGCCCCGCGGGc                                                                   >  5:167915/1‑135 (MQ=255)
                                                                        aaTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                           <  2:301669/135‑1 (MQ=255)
                                                                        aaTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGACCCCGATGGGGAGCTGGGCGTTCAGGCCGCCCACCGTGCCGTCGAcc                                                           >  2:227964/1‑134 (MQ=255)
                                                                               cgATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCAGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCAGCGGGCCGTCGCCAGCCAGGc                                                    >  2:129057/1‑135 (MQ=255)
                                                                                  tGCCCAGCGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                 <  8:282016/135‑1 (MQ=255)
                                                                                  ggCCCAGCGGGCCAAGACCTTCCACCGCCTCAACCCGCCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                 <  3:161196/134‑1 (MQ=255)
                                                                                  ggCCCAGAGCGGCAAGACCTTCGCCAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                 <  7:288356/134‑1 (MQ=255)
                                                                                    cccAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGCTGCAGCCGATGTGGAGCTGGCCGTTAAGGCCGCCCCGCGTGCCTTCGCCAGCCAGGCGTGgc                                               >  6:149196/1‑135 (MQ=255)
                                                                                       agagCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCACGCCAGGCCGCAACTGCCCGCCGCCGCTGAAGCCGAAGTGGAGCTGGCCGGTGACGCCGCCCAGCCTGCCGTCGGCAGCCAGGCGTGGCgca                                            >  2:109961/1‑135 (MQ=255)
                                                                                        gcgCGGCACGACCTTCGACAGCCTCACCCCGCCCCCCCCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                           <  2:215528/133‑1 (MQ=255)
                                                                                          gCGGCAACCCCTCCCCCACCCTCACCCCGCCCCCCCCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCtcagca                                         <  5:334521/135‑1 (MQ=255)
                                                                                                aGACCTTCGACAGCCGCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCg                                   <  2:322237/135‑1 (MQ=255)
                                                                                                    ctccGACAGCCGCACCCCGCCCCCCGCCCCGGCGGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAc                               <  1:249794/132‑1 (MQ=255)
                                                                                                     ttCGACAGCCTCAACCCGGCCACGGCCCTGGCGTGGGCGCAACTGCCCGACGCCGATGAAGCCGAGGGGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  8:193943/135‑1 (MQ=255)
                                                                                                     ttCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  8:122408/135‑1 (MQ=255)
                                                                                                        gccccccTCACCCCGCCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgc                            <  3:24299/129‑1 (MQ=255)
                                                                                                         aCAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAATCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCCTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgccc                          >  4:26552/1‑135 (MQ=255)
                                                                                                         aCAGCCTCAACCCGGCCACCGCCCACGCCTGGGCGCAACAGCCCCCCGCCCAGGAAGCCGATGTGGAGCTGGCCCTTCATGCCGCCCACCGTGCCGTATACACCAAAGCCTGGCCCACCAAAACAGCAAcccccc                          >  4:364297/1‑135 (MQ=255)
                                                                                                            gCCTCCACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  7:375047/135‑1 (MQ=255)
                                                                                                            gCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  2:6697/135‑1 (MQ=255)
                                                                                                            gCCCCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  1:45808/135‑1 (MQ=255)
                                                                                                                cAACCCGGCCACCGCCCAGGCCTGTGCGCAACTGccc                                                                                                                     >  4:356668/1‑37 (MQ=255)
                                                                                                                cAACCCGGCCACCGCCCAGGCCTGTGCGCAACTGccc                                                                                                                     <  3:356668/37‑1 (MQ=255)
                                                                                                                cAACCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGACGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCCCGGCac                   >  7:273975/1‑134 (MQ=255)
                                                                                                                    ccGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg               <  6:270498/135‑1 (MQ=255)
                                                                                                                    ccGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg               <  8:109756/135‑1 (MQ=255)
                                                                                                                    ccGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg               <  2:119112/135‑1 (MQ=255)
                                                                                                                      ggCCACCGCCCTGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct             <  5:236111/135‑1 (MQ=255)
                                                                                                                      ggCCACCGCCCAGGCGTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct             <  6:180592/135‑1 (MQ=255)
                                                                                                                          aCCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTAAGGCCCCCCCGCGTGCCGTCGGCAGCAAGGCGTGGCGCAGCCTCACCGCTACCCCCCCCGGCCAGCTGCTGcgg         >  4:253465/1‑134 (MQ=255)
                                                                                                                                 aGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCAGCTGCGTCGTCTgg  >  7:288930/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                      
CGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  KT2440/4135478‑4135741

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