Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 477,777 A→T 6.8% A109A (GCT→GCA)  glpE ← Thiosulfate sulfurtransferase GlpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440477,7770AT 6.8% 183.5 / 4.4 59A109A (GCT→GCA) glpEThiosulfate sulfurtransferase GlpE
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base T (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTCCGTGATA  >  KT2440/477643‑477908
                                                                                                                                      |                                                                                                                                   
ggCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                                                     <  4:194890/135‑1 (MQ=255)
 gCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAg                                                                                                                                    >  2:268665/1‑135 (MQ=255)
     gAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGtt                                                                                                                                >  6:168942/1‑135 (MQ=255)
      aGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTg                                                                                                                               >  1:60358/1‑135 (MQ=255)
           aaGCGTAGTTGGGTATCTGCAAGGTAAAGGGGTGCGAAAGGCAGAAGTTCCGAACGACCGCTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCt                                                                                                                           <  4:161718/134‑1 (MQ=255)
            aGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCt                                                                                                                           <  5:195526/133‑1 (MQ=255)
              cGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                                                     >  8:125854/1‑121 (MQ=255)
              cGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                                                     <  7:125854/121‑1 (MQ=255)
              cGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTCGGc                                                                                                                       >  1:80766/1‑135 (MQ=255)
               gTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCgg                                                                                                                     >  4:77955/1‑136 (MQ=255)
               gTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCgg                                                                                                                     >  4:295999/1‑136 (MQ=255)
                   ttGGTTATCTGCAAGGTAAAGGGGGGACAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGtc                                                                                                                   <  7:100759/134‑1 (MQ=255)
                     ggTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCg                                                                                                                >  1:112432/1‑135 (MQ=255)
                                aGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTc                                                                                                    <  4:162631/136‑1 (MQ=255)
                                 ggTAACGGGGTGATAAAGGAATAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTc                                                                                                    <  1:187217/135‑1 (MQ=255)
                                    aaaGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTcgcg                                                                                                 >  5:176952/1‑135 (MQ=255)
                                         ggTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGcc                                                                                             >  8:318449/1‑134 (MQ=255)
                                         ggTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCca                                                                                            >  4:360671/1‑135 (MQ=255)
                                              gaAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCg                                                                                                            <  4:97055/114‑1 (MQ=255)
                                              gaAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCg                                                                                                            >  3:97055/1‑114 (MQ=255)
                                              gaAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGtt                                                                                       <  2:60358/135‑1 (MQ=255)
                                              gaAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGtt                                                                                       >  7:40481/1‑135 (MQ=255)
                                              gaAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGtt                                                                                       >  6:15290/1‑135 (MQ=255)
                                                aaGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTc                                                                                      >  6:132134/1‑134 (MQ=255)
                                                 aGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                                   <  3:252442/136‑1 (MQ=255)
                                                 aGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGa                                                                                    <  5:86395/135‑1 (MQ=255)
                                                  ggCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                                   <  1:148751/135‑1 (MQ=255)
                                                  ggCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                                   <  7:364978/135‑1 (MQ=255)
                                                  ggCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                                   <  8:40481/135‑1 (MQ=255)
                                                      gAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                                   <  7:47293/131‑1 (MQ=255)
                                                      gAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                                   >  8:47293/1‑131 (MQ=255)
                                                       aaGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCTGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTc                                                                              >  8:127200/1‑135 (MQ=255)
                                                       aaGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCTGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTc                                                                              >  8:127202/1‑135 (MQ=255)
                                                       aaGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTc                                                                              <  2:241986/135‑1 (MQ=255)
                                                       aaGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTc                                                                              <  2:112432/135‑1 (MQ=255)
                                                          ttACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCTCGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGt                                                                           <  5:347410/135‑1 (MQ=255)
                                                          ttACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCTGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGt                                                                           <  7:127200/135‑1 (MQ=255)
                                                          ttACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCTGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGt                                                                           <  7:127202/135‑1 (MQ=255)
                                                          ttACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGt                                                                           <  6:7922/135‑1 (MQ=255)
                                                             cgaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTcc                                                                         <  1:198680/134‑1 (MQ=255)
                                                               aacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACg                                                                      >  8:67471/1‑135 (MQ=255)
                                                               aacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACg                                                                      >  8:61654/1‑135 (MQ=255)
                                                               aacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACg                                                                      >  7:182276/1‑135 (MQ=255)
                                                                acgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCt                                                                    <  3:295999/136‑1 (MQ=255)
                                                                 cgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCt                                                                    <  5:311428/135‑1 (MQ=255)
                                                                    ccGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                                                     >  4:184171/1‑67 (MQ=255)
                                                                    ccGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                                                     <  3:184171/67‑1 (MQ=255)
                                                                      gTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                                                     <  1:273760/65‑1 (MQ=255)
                                                                      gTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                                                     >  2:273760/1‑65 (MQ=255)
                                                                       ttCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTacac                                                              >  6:196090/1‑135 (MQ=255)
                                                                                 aaGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGaa                                                     <  5:261821/134‑1 (MQ=255)
                                                                                   gCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAATAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAAcc                                                  <  5:156890/135‑1 (MQ=255)
                                                                                     cGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCt                                                <  8:276386/135‑1 (MQ=255)
                                                                                         aGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGcc                                            <  2:244593/135‑1 (MQ=255)
                                                                                         aGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCt                                           <  3:280503/136‑1 (MQ=255)
                                                                                                  ggcgCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTa                                   >  5:323484/1‑135 (MQ=255)
                                                                                                                        aTAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGaa             >  2:319511/1‑135 (MQ=255)
                                                                                                                        aTAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGaa             >  1:164899/1‑135 (MQ=255)
                                                                                                                        aTAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAAt            >  3:92893/1‑136 (MQ=255)
                                                                                                                          aaaaaaTTATTCAGCGTAGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACATGTCGGAGAAACCCTGGCCTACCAGCTAGGCGGCAGCACTTTGGCTGGAAtt           >  5:210289/1‑135 (MQ=255)
                                                                                                                            aaaaTTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTc         >  1:81866/1‑135 (MQ=255)
                                                                                                                            aaaaTTATTCAGCGTTGCCATGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTCGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTc         >  1:11191/1‑135 (MQ=255)
                                                                                                                                   ttCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTCCGTGATa  <  3:77955/135‑1 (MQ=255)
                                                                                                                                   ttCAGCGTTGCCCTGGGCGGGCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTCCGTGATa  <  3:54935/135‑1 (MQ=255)
                                                                                                                                    tCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTCCGTGATa  <  7:318449/134‑1 (MQ=255)
                                                                                                                                    tCAGCGTTGCCCTGGGCGGGCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTCCGTGATa  <  8:182276/134‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGGAATTTCCGTGATA  >  KT2440/477643‑477908

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