Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,729,011 T→G 25.8% Q190P (CAG→CCG)  csdA ← Cysteine desulfurase CsdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,729,0110TG25.8% 89.0 / 4.3 76Q190P (CAG→CCG) csdACysteine desulfurase CsdA
Reads supporting (aligned to +/- strand):  ref base T (14/34);  new base G (17/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 3.37e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCAT  >  KT2440/1728884‑1729143
                                                                                                                               |                                                                                                                                      
cACGCCGACGCCGGCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCgg                                                                                                                               <  8:395323/136‑1 (MQ=255)
 acgccgacgccgTCCGGACCGTACAGCTTGTGGCGGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCgg                                                                                                                               <  4:138135/135‑1 (MQ=255)
   gccgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCCTCGACCACCGTCAGCGCGCCCAGGGCGCggggc                                                                                                                            >  8:294484/1‑134 (MQ=255)
   gccgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCa                                                                                                                            <  8:423626/136‑1 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCCTCGACCACCGTCAGCGCCCCCTGGGCGCGGGCagg                                                                                                                          >  8:382971/1‑134 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGGGCGGGCATg                                                                                                                          >  7:251614/1‑136 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGGGCCGGCAtt                                                                                                                          >  8:34711/1‑135 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCAGGGCGCGGGCAt                                                                                                                           >  5:126880/1‑135 (MQ=255)
      gacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCCCCACGCCCTGGGCACCCTCGACCACCGTCAGCGCGCCCAGGGCGGCGGGctgg                                                                                                                          >  8:392596/1‑133 (MQ=255)
          ccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGcc                                                                                                                      <  6:26582/135‑1 (MQ=255)
           cgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATTCCCCACGCCCTGGGCCCCATCGCACACCGTCAGGGCGACCGGGTCGGGGTCGTGGGGcag                                                                                                                    >  8:288350/1‑136 (MQ=255)
           cgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGccc                                                                                                                     >  4:189850/1‑134 (MQ=255)
                ggACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCCTGGGCCcgccgc                                                                                                                >  5:429973/1‑135 (MQ=255)
                ggACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCCCCCAGGGCGCGGGCATGGGCCcgcagc                                                                                                                >  4:409146/1‑135 (MQ=255)
                  aCCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAgcgc                                                                                                              <  1:275999/135‑1 (MQ=255)
                        cAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAg                                                                                                        <  7:34711/135‑1 (MQ=255)
                         aGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCg                                                                                                     >  8:290765/1‑136 (MQ=255)
                          gCTCGTGGCTGGAAAACACTCAGACGTCGCCGCCCAGGTGCGGCATGTCCTGGCGACCATGCACCCCGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCg                                                                                                     <  4:256483/135‑1 (MQ=255)
                              gtgGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCCCCCTGGGCCCGGGCAGGGGCCAGCAGCGCAGGCAGCGGTTg                                                                                                 >  1:236356/1‑135 (MQ=255)
                                ggCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGc                                                                                                <  4:387912/134‑1 (MQ=255)
                                ggCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGGCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGc                                                                                                <  5:34896/134‑1 (MQ=255)
                                 gCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCCTGGGCCCGCCCCGCGGGCCGCCGTTGGCa                                                                                              >  1:117055/1‑135 (MQ=255)
                                 gCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGCCCACCGTCAGCGCGCCCTGGGCGCGGGCCTGGGCCAGCAGCGCTGGCAGAGGTTGccc                                                                                              >  2:96009/1‑134 (MQ=255)
                                 gCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAg                                                                                             >  7:351344/1‑136 (MQ=255)
                                  cTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAg                                                                                             <  3:362643/135‑1 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCCCGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGGTGCCAg                                                                                             >  2:291591/1‑134 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGCCCACCGTCAGCGCGCCCAGGGCGCGGGCCTGGGCCAGCAGCGCTGGCCCCGGTTGCCAg                                                                                             >  2:345262/1‑134 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCCTGGGCCCGCCGCGCTGGCAGCGGGTGCCAgg                                                                                            >  3:151173/1‑135 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCCGGGCGCCGGGAAGGGCCCCCAGCGCAGGCCGCGGGTGCCac                                                                                             >  2:325019/1‑133 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCCTCAGCCCGCCCGGGGGGCGGGCCTGGGCCCCCAGCGCAGGCCGCCGGTTGcagg                                                                                             >  1:335836/1‑132 (MQ=255)
                                     gAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGACCTGGGGACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAc                                                                                         <  8:311324/136‑1 (MQ=255)
                                      aaaaCACATAGCAGTCGCAGCCCAGTTGCTGCATGGCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAc                                                                                         <  3:397068/135‑1 (MQ=255)
                                      aaaaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGt                                                                                           <  4:33329/133‑1 (MQ=255)
                                      aaaaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAc                                                                                         <  4:60889/135‑1 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAccc                                                                                       >  4:176134/1‑135 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCCAGGGCCCGCAGCGCGGGCCGCGGGTGCCACGTACCCa                                                                                      >  8:75043/1‑136 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCCCGCCCGGGGGGGCGGCATGGGCCAGCCGCGCTGGCCGCCGTTTCCAGGTAACCa                                                                                      >  8:183545/1‑136 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGAGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCCACAAGGGCGCGCGCAAGGGCCATCCCAGCAGGGAGCGGTTGCCAGGAAccc                                                                                       >  1:22761/1‑135 (MQ=255)
                                         acacATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCa                                                                                      <  4:98156/135‑1 (MQ=255)
                                                 aaGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGtt                                                                               <  4:277374/134‑1 (MQ=255)
                                                       cAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCATCGACCCCCGTCAGCGCGCCCGGGGGGCGGGGCTGGGCCCGCCGCGCCGGCCGCGGGGGGCCGGTGCCCAGCCCCTGGGTCAgc                                                                        >  3:2494/1‑135 (MQ=255)
                                                       cAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCCGCGCTGGCAGCGGTTGCCACGTACCCAGCACGTTGGACCgct                                                                       >  8:279629/1‑136 (MQ=255)
                                                       cAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACAACGCCCTGGGCACCAACGACCCCCGTCATCGCGCCCTCGGCGGGGGCATGGGCCCCCAGCCCGGGACCCGCTTGGCAGGGACCCCGCACGTGGGGCCgct                                                                       >  8:170421/1‑136 (MQ=255)
                                                          cccAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCCCCGTCAGCGCGCCCGGGGCGCGGGCCTGGGCCAGCCGCGCGGGGAGCGGTTGCCAGGTACCCCGCACGTTGGGCAgcggg                                                                     >  3:108135/1‑135 (MQ=255)
                                                          cccAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGGATGGGCCCGCCGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGGTGGACCgcgggc                                                                    >  7:129330/1‑136 (MQ=255)
                                                               ttGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGAt                                                                <  6:177087/135‑1 (MQ=255)
                                                               ttGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGATGGACAGCTGGCTGAt                                                                <  4:290800/135‑1 (MQ=255)
                                                                 gctgcATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGCCCACCGTCAGCGCGGCCGGGGCGCGGGCCTGGGCCCGCCGCGCTGGCCGCGGTTTGCAGGTACCCCGCAAGTTGGACAGAGGGGTGATgg                                                              >  2:362497/1‑135 (MQ=255)
                                                                       tGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCCGGGCGCGGGCCTGGGCCAGCCGCGCGGGCAGCGGTTGCCAGGTAACCCGCACGTTGGCCAGCTGGCGGAGGGCGAGcc                                                        >  2:380814/1‑134 (MQ=255)
                                                                          ccTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCCTGGGCCAGCCGCGGTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGCCAGCTGGCTGATGgcgcgcccgcg                                                    >  8:297002/1‑131 (MQ=255)
                                                                            tGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGGGGGCATGGGCCCGCAGCGCTGGGCGCGGTTGCCAGGTACCCAGCACGGTGGACAGCGGGCTGGTGGCGCGCCCGCgg                                                   >  6:366324/1‑135 (MQ=255)
                                                                                       caccacGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCg                                       <  7:75769/136‑1 (MQ=255)
                                                                                       caccacGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGCGCGGGCAGGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGcc                                        <  5:267354/135‑1 (MQ=255)
                                                                                        accacGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGcc                                        <  1:267626/134‑1 (MQ=255)
                                                                                              cccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                                <  7:290765/136‑1 (MQ=255)
                                                                                              cccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGt                                 <  6:126880/135‑1 (MQ=255)
                                                                                              cccTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                                <  7:382971/136‑1 (MQ=255)
                                                                                               ccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                                <  4:184574/135‑1 (MQ=255)
                                                                                                  gggCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCCCGGTGGACCGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGGCGCTCAGTTTCa                             >  2:431098/1‑135 (MQ=255)
                                                                                                   ggCACCCTCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAg                            <  5:123477/135‑1 (MQ=255)
                                                                                                   ggCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCTATCAGTTGCAg                            <  4:108135/135‑1 (MQ=255)
                                                                                                   ggCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  8:351344/136‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  3:409146/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  3:439690/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  5:423450/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGGGCGCGGGCCGATCAGTTGCAGt                           <  4:72894/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCCGCGGTTGCCAGGTACCCAGCCCGGTGGACAGCTGGCTGATGGCGCGCCCTCGGGTGCGCGGGCCGCTCAGTTGCAGTg                          >  8:41551/1‑136 (MQ=255)
                                                                                                      aCCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTTCCCAGCACGTTGGACAGCTGGCTGATGGCGCGCACTCGGGTGCGCGGGCCGCTCAGTTGCAGTGc                         >  2:329457/1‑135 (MQ=255)
                                                                                                           cGACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCCGCGGTTTGCAGGTACACAGCACCTTGGACACCTGGCTGATGGCGAGCACTCGGGGGCGCGGGCCGATCAGTTGCAGTGcctggt                    >  4:23377/1‑135 (MQ=255)
                                                                                                            gACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCCGCGGTTGCCAGGTACCCAGCACGTTGGTCAGCTGGCTGCTGGCGCGCCATCGGGTGCGCGGGCCGCTCAGTTTCAGTGcctggtc                   >  4:298713/1‑135 (MQ=255)
                                                                                                                 ccGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  7:294484/136‑1 (MQ=255)
                                                                                                                  cggcAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGGTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  1:325019/132‑1 (MQ=255)
                                                                                                                  cGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  3:189850/135‑1 (MQ=255)
                                                                                                                  cGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  1:345262/135‑1 (MQ=255)
                                                                                                                     cAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATa         >  7:336684/1‑136 (MQ=255)
                                                                                                                            cggtGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGGGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCAt  <  8:400541/133‑1 (MQ=255)
                                                                                                                            cccGGGGCGCGGGCAGGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCa   <  4:28167/135‑1 (MQ=255)
                                                                                                                               |                                                                                                                                      
CACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCAT  >  KT2440/1728884‑1729143

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: