Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,000,955 A→T 6.9% L189Q (CTA→CAA)  ggt_1 ← Glutathione hydrolase proenzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,000,9550AT 6.9% 211.8 / 2.6 67L189Q (CTA→CAA) ggt_1Glutathione hydrolase proenzyme
Reads supporting (aligned to +/- strand):  ref base A (36/26);  new base T (3/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTGTGCCAACTCAGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCACCCCGGGCA  >  KT2440/4000837‑4001084
                                                                                                                      |                                                                                                                                 
ggTCTGTGCCAACTCAGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTc                                                                                                                   >  4:89388/1‑135 (MQ=255)
ggTCTGTGCCAACTCAGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTc                                                                                                                   <  4:95088/135‑1 (MQ=255)
       gCCAACTCAGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGaaa                                                                                                              <  2:29652/133‑1 (MQ=255)
       gCCAACTCAGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAccc                                                                                                            >  2:121111/1‑135 (MQ=255)
       gCCAACTCAGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAAcc                                                                                                            >  3:253122/1‑135 (MQ=255)
       gCCAACTCAGGGTTGCGCAGGATTATACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAAcc                                                                                                            >  7:283728/1‑135 (MQ=255)
       gCCAACTCAGGGTGGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGaaa                                                                                                              <  7:161602/133‑1 (MQ=255)
       gCCAACTCAGGGTGGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAAcc                                                                                                            <  2:307233/135‑1 (MQ=255)
               aGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGcc                                                                                                    >  2:297376/1‑135 (MQ=255)
                   ttGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGACGCCGTTg                                                                                                >  1:139703/1‑135 (MQ=255)
                   ttGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTg                                                                                                >  5:3281/1‑135 (MQ=255)
                      cgcAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGt                                                                                             <  5:297408/135‑1 (MQ=255)
                       gcAGGATTGTACCCACTGCCAGCGGCCTGCCCTGTTGATCCAGGAAGTAACGGGCCATTGCGGGTGAGCTTGCAATGTAGGGGTCCCGGGCCAGTAGTGTGTGCAGGGGCTCTGAAACCGCGAAGCCGTTGCGt                                                                                             <  8:30715/134‑1 (MQ=255)
                         aGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGcc                                                                                          >  2:210108/1‑135 (MQ=255)
                         aGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGc                                                                                           >  7:51237/1‑134 (MQ=255)
                         aGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGc                                                                                           >  1:211121/1‑134 (MQ=255)
                         aGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGc                                                                                           >  2:77003/1‑134 (MQ=255)
                         aGGATTGTACCCACGGCCAGCGGCCTGCACTGTTCATCCAGGAAGTAAAGGTCCATTGCGGTTGAGCGGGCAATGTAGGCGTCCCCGCCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGc                                                                                           >  2:331122/1‑134 (MQ=255)
                                   ccACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTTGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATgg                                                                                >  5:135919/1‑135 (MQ=255)
                                   ccACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTTGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATgg                                                                                >  5:135915/1‑135 (MQ=255)
                                        gCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATg                                                                                 >  5:301653/1‑129 (MQ=255)
                                        gCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATg                                                                                 <  6:301653/129‑1 (MQ=255)
                                        gCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCggg                                                                           >  7:176816/1‑135 (MQ=255)
                                        gCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGACAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCggg                                                                           >  3:6223/1‑135 (MQ=255)
                                          cAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGc                                                                         <  5:173574/135‑1 (MQ=255)
                                          cAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGc                                                                         >  3:294382/1‑135 (MQ=255)
                                          cAGAGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCGGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCTTTGCGTGCCAAGGCGATGGCCGGGGc                                                                         <  3:278613/135‑1 (MQ=255)
                                             cGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCgca                                                                      >  1:258369/1‑133 (MQ=255)
                                             cGGCCTGCCCTGGTCATCCAGGAAGTAACGGGCCAGTGCGGGTGCGCGGGCAATGTAGGGGTCCCCGGACAGTAGTGTGTGCCGGAGCTCGGACACCGCGAAGCCGTTGCGTGCCACGGCGATGGCCGGGGCGaa                                                                      >  4:84268/1‑135 (MQ=255)
                                               gCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGGCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACa                                                                    <  1:102636/135‑1 (MQ=255)
                                                            aTCCAGGAAGTCACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAgg                                                        <  8:176546/134‑1 (MQ=255)
                                                            aTCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAgg                                                        >  1:89769/1‑134 (MQ=255)
                                                                aGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAAc                                                   <  2:83795/135‑1 (MQ=255)
                                                                      tAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCg                                             >  4:98899/1‑135 (MQ=255)
                                                                       caCGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGt                                            <  6:197485/134‑1 (MQ=255)
                                                                         cGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCAAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGtg                                           <  2:89769/134‑1 (MQ=255)
                                                                                  gCGGGTGAGCGGGCAATGTAGGGGGCCCCGGCCAGTAGTGTGGGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTg                                  <  5:209019/134‑1 (MQ=255)
                                                                                   cGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAGCAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTggg                                >  7:261240/1‑135 (MQ=255)
                                                                                    gggTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCTGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGc                               >  7:269265/1‑135 (MQ=255)
                                                                                    gggTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGc                               >  5:260344/1‑135 (MQ=255)
                                                                                        gAGCGGGCAATGTAGGGGTCCCCGGCAAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGc                           >  5:126958/1‑135 (MQ=255)
                                                                                         aGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGATCAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCt                          <  8:261240/135‑1 (MQ=255)
                                                                                         aGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCt                          <  6:327080/135‑1 (MQ=255)
                                                                                         aGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCt                          >  1:278294/1‑135 (MQ=255)
                                                                                         aGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTACt                          <  1:121111/135‑1 (MQ=255)
                                                                                           cGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTc                        <  4:253122/135‑1 (MQ=255)
                                                                                             ggCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAg                      >  4:49510/1‑135 (MQ=255)
                                                                                             ggCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAg                      <  1:297376/135‑1 (MQ=255)
                                                                                                aaTGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGc                   >  4:40755/1‑135 (MQ=255)
                                                                                                 aTGTAGGGGTCCCCGGCCAGTAGTGTGTGCTGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCa                  <  8:269265/135‑1 (MQ=255)
                                                                                                  tGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCAc                 <  1:128779/135‑1 (MQ=255)
                                                                                                   gTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCAc                 >  3:93494/1‑134 (MQ=255)
                                                                                                   gTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGCAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGGGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACg                <  6:320944/135‑1 (MQ=255)
                                                                                                     aGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACgc               <  8:214590/134‑1 (MQ=255)
                                                                                                          tCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGcc                                              >  2:266622/1‑98 (MQ=255)
                                                                                                          tCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGcc                                              <  1:266622/98‑1 (MQ=255)
                                                                                                             ccGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCACCCCg      >  4:66268/1‑135 (MQ=255)
                                                                                                              cGGCCAGTTGTGGGTGCAGGCGCTCGGAAACCGGGAGGCCGTTGCGTGCCAAGGCGAGGGCGGGGGCGAACTGGGCGGGCCAGGGCAACTTGCCGGGCTGCTCGGGGGCCTGCTTCTGGGCACGCAGCACCCCgg     <  3:98899/135‑1 (MQ=255)
                                                                                                               ggCCAGTTGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAAc                                                                     <  3:112308/70‑1 (MQ=255)
                                                                                                               ggCCAGTTGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAAc                                                                     >  4:112308/1‑70 (MQ=255)
                                                                                                               ggCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCCCCCCggg    >  8:9501/1‑135 (MQ=255)
                                                                                                               ggCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCACCCCggg    >  8:8433/1‑135 (MQ=255)
                                                                                                               ggCCAGTAGTGTGTGCAGGCGCTCGGAAACAGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCACCACggg    >  1:105301/1‑135 (MQ=255)
                                                                                                                 ccAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAg                                                                                       <  4:198200/50‑1 (MQ=255)
                                                                                                                 ccAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAg                                                                                       >  3:198200/1‑50 (MQ=255)
                                                                                                                 ccAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCACCCCGGGCa  >  4:161473/1‑135 (MQ=255)
                                                                                                                 ccAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCACCCCGGGCa  >  3:61561/1‑135 (MQ=255)
                                                                                                                      |                                                                                                                                 
GGTCTGTGCCAACTCAGGGTTGCGCAGGATTGTACCCACGGCCAGCGGCCTGCCCTGTTCATCCAGGAAGTAACGGGCCATTGCGGGTGAGCGGGCAATGTAGGGGTCCCCGGCCAGTAGTGTGTGCAGGCGCTCGGAAACCGCGAAGCCGTTGCGTGCCAAGGCGATGGCCGGGGCGAACAGGGCGTGCCAGGGCAACTTGCCGTGCTGCTCGTGGGCCTGCTTCAGGGCACGCAGCACCCCGGGCA  >  KT2440/4000837‑4001084

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