Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,544,643 T→G 22.3% Q412P (CAG→CCG)  pykA ← Pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,544,6430TG22.3% 80.4 / 3.1 74Q412P (CAG→CCG) pykAPyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base T (13/41);  new base G (13/5);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 4.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTGGCAG  >  KT2440/1544512‑1544769
                                                                                                                                   |                                                                                                                              
cccAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGt                                                                                                                            >  8:14605/1‑136 (MQ=255)
  cAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTgg                                                                                                                          <  8:32144/136‑1 (MQ=255)
  cAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGGTCGTCGACCGCCGCCTGGCTCCCCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTgg                                                                                                                          <  7:310478/136‑1 (MQ=255)
   aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGgc                                                                                                                         >  8:268580/1‑136 (MQ=255)
   aGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGgggggg                                                                                                                          >  5:210723/1‑135 (MQ=255)
   aGTCACCTTGCTCCACCAGGCCACGTTTCAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTgg                                                                                                                          <  1:310892/135‑1 (MQ=255)
      cACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGcgcg                                                                                                                      <  8:380417/136‑1 (MQ=255)
        ccTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGCAGGCAATCGGGTACACGCCCCGGAACATCGAGGCGCGCGCCGGGGGGTCGCggg                                                                                                                    >  7:173344/1‑135 (MQ=255)
        ccTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGCGGCGCGCGGCAGGGTGGCGCgg                                                                                                                     >  1:358479/1‑135 (MQ=255)
        ccTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCgg                                                                                                                     >  4:29815/1‑135 (MQ=255)
        ccTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCgg                                                                                                                     >  2:7052/1‑135 (MQ=255)
             cTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTgcgg                                                                                                                <  8:99855/135‑1 (MQ=255)
              tCCCCCCGGCCCCGTTTGCGCAGGTCGGCGACCGCCGCCTGGCTCCCCTTATCGGCCGGCAGCGCGGGCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGCGGCGCGCGCCTGGGTGGCGCGGTgcggc                                                                                                               <  3:295586/135‑1 (MQ=255)
              tCCACCAGGCCACGTTTGAGCAGGTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGTAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTgcgg                                                                                                                <  3:71609/134‑1 (MQ=255)
                  ccaGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCCGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGACGCGCGCGCCCGGGTGGCGCGGGggggcgcga                                                                                                           >  3:321817/1‑132 (MQ=255)
                  ccaGGCCACGTTTGAGCAGGTCGGCGCCCGCCGCCGGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACa                                                                                                           <  1:316310/135‑1 (MQ=255)
                    aGGCCACGTTTGCGCCGTTCGTCCCCCGCCGGCTGGCGCACCTTATCGGCCGGCCGCGCGGCCGGGGCGAAGGCAATCGGGAACACGCCACGGAACACCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAgg                                                                                                         <  5:440367/135‑1 (MQ=255)
                     ggCCACGTTTGAGCAGTTCGTCGCCCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTCCACGCCACGGAACATCGCGGCGCGCGCCTGGGGGGCGCGGGGCGGCGCCCGGGc                                                                                                       >  8:134245/1‑136 (MQ=255)
                     ggCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCCCCGCGGTGGCGCGGTGCGGCCCCAGGGc                                                                                                       >  8:423412/1‑136 (MQ=255)
                      gCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGc                                                                                                       >  6:169218/1‑135 (MQ=255)
                      gCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGCGGTGGCGCGGTGCGGCGCCAGGGc                                                                                                       >  2:120227/1‑135 (MQ=255)
                        cACGTTTGAGCAGTTCGTCGACCCCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGCCGAACGCAATCGGGTTCACAGCACGGAACATCGACGCCACAGCCGGGCTGGCCCGGTGCGGCGGCCAGGCGa                                                                                                     >  6:82608/1‑135 (MQ=255)
                        cACGTTTGAGCAGTTCGGCGACCGCCGCCTGGCTCAGCTTCTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTCCACGCCACGGAACATCGAGGCTCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                    <  7:257923/136‑1 (MQ=255)
                         aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                    <  1:242735/135‑1 (MQ=255)
                         aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGCAg                                                                                                   >  8:352383/1‑136 (MQ=255)
                         aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGGCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                    <  2:337691/135‑1 (MQ=255)
                         aCGTTTGAGCAGTTCGTCGACCGCCGCCTGGCGCACCTTCTCGGCCGGCAGCGGGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGaa                                                                                                    <  6:383458/135‑1 (MQ=255)
                                        gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGCAGGCAATCGGGTACACGCCCCGGAACATCGCGGCGCGCGCCTGGGTGGTGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGa                                                                                     >  5:105802/1‑135 (MQ=255)
                                        gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCCCGGAACATCGAGGCGCGCGGCGGGGTGGCGCGGGGCGGCGGCAGGGCGAAGATCGGTACATGCgg                                                                                     >  4:76039/1‑134 (MQ=255)
                                        gtcgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGCCAGGGCGCAGCTCGGTACATGCGa                                                                                     >  4:199611/1‑135 (MQ=255)
                                          cgACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAAc                                                                                   <  4:85473/135‑1 (MQ=255)
                                             ccgccgccgGGCTCACGTTCTGGGCCGGCCGGGCGGGCGGGTCGATGGCAATCGGGAACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCa                                                                                <  1:398437/135‑1 (MQ=255)
                                             ccgccgCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCa                                                                                <  6:210723/135‑1 (MQ=255)
                                              cgccgcCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAg                                                                               <  1:2728/135‑1 (MQ=255)
                                                 cgccTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGGACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGgcgc                                                                           <  7:347394/136‑1 (MQ=255)
                                                  gccTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGTAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGgcgc                                                                           <  3:406195/135‑1 (MQ=255)
                                                          ccgcTTAGCGGCCGGCCGCGCGGCCGGGTGTAAGGCAATCGGGTCCCCGCCACGTAACATGGAGGCGCGCGCCTGGGGGGCGCGGTGCGGCGACAGGGCGAAGAGCGGTACATGCGAACGCAGTCGCGACATGAt                                                                   <  4:436916/132‑1 (MQ=255)
                                                           aCCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACTTGATCa                                                                 <  7:423412/136‑1 (MQ=255)
                                                            cgttAGCGGCGGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                 <  2:201819/133‑1 (MQ=255)
                                                            ccTTCTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGGGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                 <  2:334841/135‑1 (MQ=255)
                                                            ccTTATCGGCCGGCTGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGCGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                 <  4:20000/135‑1 (MQ=255)
                                                            ccTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                 <  4:368462/135‑1 (MQ=255)
                                                            ccTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCa                                                                 <  4:230184/135‑1 (MQ=255)
                                                              ttATCGGCCGGCAGCGCGGCGGGGTCGAAGGCAATCCGGTACACGCCACGGACCATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCCCCCGGGCGAAGATCGGTACCTGCGAACGCCGGCCCGACCTGATCAGc                                                               >  6:103004/1‑135 (MQ=255)
                                                                aTCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGGGGCGCCGTGCCGCGGCAGGGCGCAGATAGGTTCAAGCGAACCCCGGCGCCCCATGGTAAGccg                                                             >  5:160146/1‑133 (MQ=255)
                                                                aTCGGCCGGCAGCGCGGCCGGGGCGAAGGCAATCGGGTACACGCCACGTAACCTCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCgg                                                             <  8:419601/135‑1 (MQ=255)
                                                                    gccggcAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGCCAGGGCGAAGATCGGTACATGCGCACGCAGGCGCGACATGATCAGCGGGgtgt                                                        >  7:331879/1‑136 (MQ=255)
                                                                          aGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGCTCGGTACATGCGCACGCAGGCGCGACATGCTCAGCGGGGTGTATCCg                                                   >  4:142685/1‑135 (MQ=255)
                                                                           gcggggcgggGGCGACGGCAATCGGCTACACCCCCCGCAACTTGGCGGCGCGCGCCTGGGGGGCGCGGGGGGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGc                                                  <  3:423168/135‑1 (MQ=255)
                                                                                      tCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAgg                                       <  3:348250/135‑1 (MQ=255)
                                                                                       cGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAggg                                      <  3:136250/135‑1 (MQ=255)
                                                                                           ggCTATCGGGTACACGCCACGGCACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAATATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCgat                                  <  5:59168/135‑1 (MQ=255)
                                                                                              aaTCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGGGGCGACAGGGCGAAGCTCGGTACATGCGCACGCCGGCGCGCCATGGTCAGCGGGGTGTAGCCGCGTTGCGTCAGGGCgatgat                               >  5:302183/1‑135 (MQ=255)
                                                                                              aaTCGGGTACACGCCACGGAACATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATc                              >  8:360381/1‑136 (MQ=255)
                                                                                                  gggTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGCCAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCt                          >  8:239440/1‑136 (MQ=255)
                                                                                                  gggTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACCTGATCAGCGGGGTGGTGCCGCTTTCGGTCATGGCGATGATCGcc                           >  4:136250/1‑135 (MQ=255)
                                                                                                   ggTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCt                          <  3:76039/135‑1 (MQ=255)
                                                                                                             aCGGAACATCGCGGCGCGCGCCTGGGGGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCgg                <  3:29815/135‑1 (MQ=255)
                                                                                                                 aaCATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt           <  7:285748/136‑1 (MQ=255)
                                                                                                                 aaCATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGCCCGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATAAGCGGGGTGTAGCCGCTTTAGGTGAGGGCGATGAAAGCCGTTTCGCCCGGGAAGt           >  8:241518/1‑136 (MQ=255)
                                                                                                                  ccATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt           <  1:120227/134‑1 (MQ=255)
                                                                                                                  agaTCGTGGGGCGGGCCTGGGTGGGGCGGTGCTGCGACTGGGCGTAGAGCGGTAGATGCGATCGCAGGCGCGACATGATGAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt           <  5:82608/133‑1 (MQ=255)
                                                                                                                  aCATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt           <  1:242087/135‑1 (MQ=255)
                                                                                                                  aCATCGAGGCGCGCGCCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGt           <  5:169218/135‑1 (MQ=255)
                                                                                                                     tCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt        <  3:142685/135‑1 (MQ=255)
                                                                                                                     gcgcggCGCGCGCCGGGGGGGCGCGGTGCGGCGACAGGGCGAAGAGCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt        <  2:195298/131‑1 (MQ=255)
                                                                                                                       gtggggcgcgcCGGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGt        <  2:420856/128‑1 (MQ=255)
                                                                                                                       gcggCGCGCGCCTGGGTGGCGCGGGGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg      <  6:20896/133‑1 (MQ=255)
                                                                                                                       gAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg      <  7:268580/135‑1 (MQ=255)
                                                                                                                       gAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg      <  8:30429/135‑1 (MQ=255)
                                                                                                                        aGGCGCGCGCGTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg      <  1:62942/134‑1 (MQ=255)
                                                                                                                        aGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTg      <  3:199611/134‑1 (MQ=255)
                                                                                                                         ggCGCGCGCCTGGGTGGCGCGGTGCGGCGCCAGGGCGCAGCTAGGGTCATGCGCACGCAGGCGCGACATGTTCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGAATATCGCCTTTACGCCCGGGAAGGGGTTg      >  4:188770/1‑133 (MQ=255)
                                                                                                                           cgcgcgcCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATAGCCTTTACGCCCGGGAAGTGGTTGGCAg  >  6:194552/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                              
CCCAGTCACCTTGCTCCACCAGGCCACGTTTGAGCAGTTCGTCGACCGCCGCCTGGCTCACCTTATCGGCCGGCAGCGCGGCCGGGTCGAAGGCAATCGGGTACACGCCACGGAACATCGAGGCGCGCGCCTGGGTGGCGCGGTGCGGCGACAGGGCGAAGATCGGTACATGCGAACGCAGGCGCGACATGATCAGCGGGGTGTAGCCGCTTTCGGTCAGGGCGATGATCGCCTTTACGCCCGGGAAGTGGTTGGCAG  >  KT2440/1544512‑1544769

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: