Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,135,607 A→C 49.6% 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,6070AC49.6% 6.0 / 34.1 63Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (14/17);  new base C (24/8);  total (38/25)
Fisher's exact test for biased strand distribution p-value = 2.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

AAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTG  >  KT2440/4135473‑4135740
                                                                                                                                      |                                                                                                                                      
aaGCCCGCACCCTCCCCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                        <  2:211306/135‑1 (MQ=255)
         cccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCCGCCACCGCCCAGGCCTGGgcg                                                                                                                             >  4:119615/1‑136 (MQ=255)
         cccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTGCGTTTCCAGATGTGCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGgcg                                                                                                                             <  4:336304/136‑1 (MQ=255)
          ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCTGGgcg                                                                                                                             >  1:291030/1‑135 (MQ=255)
                cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCTACAGCCTCAACCCGGCCACCGCCCCGGCCAGGGCGCCACGGc                                                                                                                      >  4:328386/1‑136 (MQ=255)
                 gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCGGCCGGGGCGCCCCTGcc                                                                                                                     >  3:150777/1‑136 (MQ=255)
                 gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCAGGGCGCCACTGc                                                                                                                      >  8:88563/1‑135 (MQ=255)
                     ccGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTCGGCGCCACTGCCCgcc                                                                                                                 >  4:110492/1‑134 (MQ=255)
                       gCCAAGGGGATTGTTGCATTACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCcgccgc                                                                                                               >  4:334103/1‑136 (MQ=255)
                       gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGCCCGGGGCGCAACTGCCCGACGc                                                                                                               >  3:158473/1‑136 (MQ=255)
                          aaGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGCCGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCGGCCAGGGCGCACCTGCCCGACGcc                                                                                                              >  4:237448/1‑134 (MQ=255)
                                    ttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACAGCCCGCCGCCGCTGAAGCCGa                                                                                                   >  7:116923/1‑135 (MQ=255)
                                     tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACAGCCCGCCGCCGATGAAGCCGag                                                                                                  >  6:337659/1‑134 (MQ=255)
                                     tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCACAGACCGCCGCCGATGAAGCCGag                                                                                                  >  5:70617/1‑134 (MQ=255)
                                      gctccccccTCGTGCTCTACCTGATGTGCATTGACGGCAACGGGGGCTAGCCCCATCGCGTCATGACTATCGACCGCCACAACCCGCCCACCCCCAGCGCCTGGCCGCAACTGCCCGACGCCGATGAAGCCGAtgt                                                                                                <  4:311732/130‑1 (MQ=255)
                                           aCCCTCGTGCGATGCCAGATGGGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACATTCGACACCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAg                                                                                            <  5:165505/135‑1 (MQ=255)
                                            cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACGGCCCGACGCCGATGAAGCCGATGTGGAGc                                                                                           >  2:145233/1‑135 (MQ=255)
                                            cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCGAGCCGCACCAGACCCACCCCGACGAAGCCGCAGTGGAGc                                                                                           >  7:310688/1‑135 (MQ=255)
                                                    cgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGCAGAAGCAGATGTGGAGCTTGCCGtt                                                                                   >  8:346073/1‑135 (MQ=255)
                                                    cgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACGGCCCGCCGCCGAGGAAGCCGATGTGGCGCTGGCCCtt                                                                                   >  6:286721/1‑135 (MQ=255)
                                                    cgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGACCATGCAGCCGATGTGGAGCTTGCCGt                                                                                    >  8:306309/1‑134 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGCAACAGCCCCACCCCGATGAACCCGATGTGGAGCTGGCCGTTCACgc                                                                              >  2:52217/1‑135 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGgcc                                                                             >  4:89329/1‑136 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCACCCCAGGGCGCAACTGCCCGACGCCGATGCAGCCCCTGTGGAGCTGGCCCTTCAGgcc                                                                             >  4:122856/1‑136 (MQ=255)
                                                         ccAGAGGTGCCTCGACGGCCACGGGCGCGATGCCCAGCGCGGCAAGACCTTCGACAGCCTCAACCCGCCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAgg                                                                               <  5:115361/134‑1 (MQ=255)
                                                          cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGgccc                                                                            >  4:206445/1‑135 (MQ=255)
                                                          cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCCACTGCCCGCCGCCGCTGGAGCCGATGTGGAGCTGGGCGTTCAGgcc                                                                             >  8:3638/1‑135 (MQ=255)
                                                                       aCGGCCAATGGCGCGAGGCCCAGAGCGGCAAGACCTTCGCCAGCCTCACCCCCGCCCCCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCtt                                                               <  3:254904/136‑1 (MQ=255)
                                                                        cgcccAAGGGCGCGATCCCCAGCGCGCCAACCCCTTCCCCACCCACACCCCGCCCCCCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCtt                                                               <  7:6003/132‑1 (MQ=255)
                                                                        cGGCCACGGGCGCGAGGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCtt                                                               <  7:301830/135‑1 (MQ=255)
                                                                             ctgggCGCGAGGCCCCGAGCGCCAACACCTCCCACCCCCCCACCCCGCCCCCCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                          <  5:101406/132‑1 (MQ=255)
                                                                             accggCGCGAGGCCCGGCGCGGCAAGACCTTCGACAGCCGCACCCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                          <  5:300508/132‑1 (MQ=255)
                                                                                 gcgcGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                     <  4:114565/136‑1 (MQ=255)
                                                                                 gcgcGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGATGACGCCGATGTGGCGCTGGCCGTTCAGGCCGCCCCGCGTGCCTTCGCCAGCAAg                                                     >  4:31547/1‑136 (MQ=255)
                                                                                    cgATGCCCAGAGCGGCAAGACCTGCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCg                                                  <  4:66373/136‑1 (MQ=255)
                                                                                     gaggCCCAGAGCGCCACGCCCTCCCCCAGCCCCACCCCGCCCCCCCCCCCGGCGTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCg                                                  <  1:5067/132‑1 (MQ=255)
                                                                                      aTGCCCAGAGCGGCAAGACCTTCGACAGCTTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                <  4:109785/136‑1 (MQ=255)
                                                                                       tGCCCAGCGCGGAACGACCTTCCACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                <  6:337687/135‑1 (MQ=255)
                                                                                       tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACTCGGCCACCGCCCCGGCCGGGGCGCACCTGCCCGGCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCGTGCCGTCGCCCTCAAGGCGTg                                                >  8:324976/1‑135 (MQ=255)
                                                                                       tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACAGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTGAGGCCGCCCCGCGTGCCGTCGACACCCAGGCGTgg                                               >  3:225414/1‑136 (MQ=255)
                                                                                       tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCAGGGCGCAACTGCCCGACGCCGATGAAGCCGCTGTGGCGCTGGGCGTTCAGGCCGCCCACCGTGCCGTCGCCAGCAAGGCGTgg                                               >  4:260155/1‑136 (MQ=255)
                                                                                       tGCCCAGAGCGCCAAGACCTGCGCCCGCCCCACCCCGCCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                <  2:193907/135‑1 (MQ=255)
                                                                                          ccAGAGCGGCAAGACCTTCGACCGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                            <  4:59362/136‑1 (MQ=255)
                                                                                           cAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                            <  8:357784/135‑1 (MQ=255)
                                                                                               gccgCAAGACCTGCGACAGCCACCACCCGCCCACCCCCCAGCCGGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcagca                                        <  1:18663/132‑1 (MQ=255)
                                                                                                 ggCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGCGGCAGCCGCTGTGGAGCTGGCCGTTCAGGCCCCCCCGCGTGCCTTCGACACCAAGGCGTGGCGCAGCATc                                      >  1:105563/1‑135 (MQ=255)
                                                                                                     acaCCTGCGACAGCCACAACCCGCCCACCCCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCGTCACCg                                  <  8:12183/133‑1 (MQ=255)
                                                                                                      gACCTTCGACAGCCTCAACCCGGCCCCCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCCGCCTCACCg                                  >  7:280206/1‑134 (MQ=255)
                                                                                                      gACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCAACAGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGGCGGTGAGGCCGCCCAGCGTGCCGTCGACAGCAAGGCGGGGCGCAGCATCACCGc                                 >  4:140974/1‑135 (MQ=255)
                                                                                                      gACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCCCTGTGGAGCTGGCCGTTCAGGCCGCCCCGCCTGCCTTCGCCACCAACGCGTGGCGCAGCCACAccc                                  >  6:180490/1‑133 (MQ=255)
                                                                                                      gACCTTCGACAGCCTCAACCCGGCCACAGCCCCGGCCTCGGCGCAACAGCCCGACGCCCATGAAGCCCATGTAGAGCATGACGATAATGCCCCCCAGCGTGCCATCGACAGCAAGGCGTGTCGCAGCATAAACGc                                 >  8:161813/1‑135 (MQ=255)
                                                                                                          tgcGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                             <  6:64282/133‑1 (MQ=255)
                                                                                                           tCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGCTGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCAGCGTGGCGTCGCCCGCAAGGCGTGGCGCCGCATCACCGCTAccgc                           >  4:302605/1‑136 (MQ=255)
                                                                                                              aCAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCCACTGCCCGCCGCCGAAGAAGCCCATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGGCAGCAAGGCGTGGCGCCGCCTCACCGCTAcccccc                         >  2:319962/1‑135 (MQ=255)
                                                                                                                     cAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCGGGCCGTCGCCACCAACGCCTGGCGCAGCAACACCGCAACCCCCCCCGGccc                  >  6:11169/1‑133 (MQ=255)
                                                                                                                         ccGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg              <  1:116873/135‑1 (MQ=255)
                                                                                                                          cgccccccGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct            <  4:159922/130‑1 (MQ=255)
                                                                                                                           gcccACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct            <  8:54465/133‑1 (MQ=255)
                                                                                                                               aCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGcgt        <  8:95567/135‑1 (MQ=255)
                                                                                                                               aCCGCCCAGGCCTGGGCGCAACGGCCCGACGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCGTGCCGTCGCCAGCAAGGCGTGGCGCAGCATCACCGCGACCGCCCGCGGCAAGCGGCTGcgg        >  1:268598/1‑134 (MQ=255)
                                                                                                                               aCCCCCCAGCCCTGGGCGCAACTGCCCCACGCCCGTGGAGCCGATGTGGCGCTGGGCGTTCAGGCCGCCCAGCGTTCCGTGGACAGCCAGGGGTGGCGCCCCATCACCCCCACCCCCCCCCGCAAgcggcgggggc       >  4:5259/1‑130 (MQ=255)
                                                                                                                                   cccAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCACCGTGCCGTCGCCAGCAAGGCGTGGCGCAGCATCACCCCTACCGCCCCCCGCCAgcggcggcggcggg    >  8:74900/1‑133 (MQ=255)
                                                                                                                                    ccAGGCCTGGGCGCAACTGCCCGACGCCGAAGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCAGCGTGCCGTCGACAGCAACGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCACGCGGCTGCGTCGTCt   >  8:356237/1‑135 (MQ=255)
                                                                                                                                     cAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTg  <  8:63710/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTG  >  KT2440/4135473‑4135740

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