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

CACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCAT  >  KT2440/1728884‑1729143
                                                                                                                               |                                                                                                                                      
cACGCCGACGCCGGCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCgg                                                                                                                               <  8:395323/136‑1 (MQ=255)
 acgccgacgccgTCCGGACCGTACAGCTTGTGGCGGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCgg                                                                                                                               <  4:138135/135‑1 (MQ=255)
   gccgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCCTCGACCACCGTCAGCGCGCCCAGGGC‑GCggggc                                                                                                                            >  8:294484/1‑134 (MQ=255)
   gccgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCa                                                                                                                            <  8:423626/136‑1 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCCTCGACCACCGTCAGCGCCCCCTGGGC‑GCGGGCagg                                                                                                                          >  8:382971/1‑134 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGG‑GCGGGCATg                                                                                                                          >  7:251614/1‑136 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGG‑GCCGGCAtt                                                                                                                          >  8:34711/1‑135 (MQ=255)
     cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCAGGGC‑GCGGGCAt                                                                                                                           >  5:126880/1‑135 (MQ=255)
      gacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCCCCACGCCCTGGGCACCCTCGACCACCGTCAGCGCGCCCAGGGCGGCGGGctgg                                                                                                                          >  8:392596/1‑133 (MQ=255)
          ccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGcc                                                                                                                      <  6:26582/135‑1 (MQ=255)
           cgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATTCCCCACGCCCTGGGCCCCATCGCACACCGTCAGGGCGACCGGGTC‑GGGGTCGTGGGGcag                                                                                                                    >  8:288350/1‑136 (MQ=255)
           cgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGccc                                                                                                                     >  4:189850/1‑134 (MQ=255)
                ggACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCCTGGGCCcgccgc                                                                                                                >  5:429973/1‑135 (MQ=255)
                ggACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCCCCCAGGGC‑GCGGGCATGGGCCcgcagc                                                                                                                >  4:409146/1‑135 (MQ=255)
                  aCCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAgcgc                                                                                                              <  1:275999/135‑1 (MQ=255)
                        cAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAg                                                                                                        <  7:34711/135‑1 (MQ=255)
                         aGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑Cg                                                                                                     >  8:290765/1‑136 (MQ=255)
                          gCTCGTGGCTGGAAAACACTCAGACGTCGCCGCCCAGGTGCGGCATGTCCTGGCGACCATGCACCCCGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑Cg                                                                                                     <  4:256483/135‑1 (MQ=255)
                              gtgGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCCCCCTGGGC‑CCGGGCAGGGGCCAGCAGCGCAGGCAG‑CGGTTg                                                                                                 >  1:236356/1‑135 (MQ=255)
                                ggCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGc                                                                                                <  4:387912/134‑1 (MQ=255)
                                ggCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGGCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGc                                                                                                <  5:34896/134‑1 (MQ=255)
                                 gCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCCTGGGCCCGCCCCGCGGGCCG‑CCGTTGGCa                                                                                              >  1:117055/1‑135 (MQ=255)
                                 gCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGCCCACCGTCAGCGCGCCCTGGGC‑GCGGGCCTGGGCCAGCAGCGCTGGCAG‑AGGTTGccc                                                                                              >  2:96009/1‑134 (MQ=255)
                                 gCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAg                                                                                             >  7:351344/1‑136 (MQ=255)
                                  cTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAg                                                                                             <  3:362643/135‑1 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCCCGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGGTGCCAg                                                                                             >  2:291591/1‑134 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGCCCACCGTCAGCGCGCCCAGGGC‑GCGGGCCTGGGCCAGCAGCGCTGGCCC‑CGGTTGCCAg                                                                                             >  2:345262/1‑134 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCCTGGGCCCGCCGCGCTGGCAG‑CGGGTGCCAgg                                                                                            >  3:151173/1‑135 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCCGGGC‑GCCGGGAAGGGCCCCCAGCGCAGGCCG‑CGGGTGCCac                                                                                             >  2:325019/1‑133 (MQ=255)
                                   tGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCCTCAGCCCGCCCGGGGG‑GCGGGCCTGGGCCCCCAGCGCAGGCCGCCGGTTGcagg                                                                                             >  1:335836/1‑132 (MQ=255)
                                     gAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGACCTGGGGACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTAc                                                                                         <  8:311324/136‑1 (MQ=255)
                                      aaaaCACATAGCAGTCGCAGCCCAGTTGCTGCATGGCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTAc                                                                                         <  3:397068/135‑1 (MQ=255)
                                      aaaaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGt                                                                                           <  4:33329/133‑1 (MQ=255)
                                      aaaaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTAc                                                                                         <  4:60889/135‑1 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTAccc                                                                                       >  4:176134/1‑135 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCCAGGGCCCGCAGCGCGGGCCG‑CGGGTGCCACGTACCCa                                                                                      >  8:75043/1‑136 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCCCGCCCGGGGG‑GGCGGCATGGGCCAGCCGCGCTGGCCG‑CCGTTTCCAGGTAACCa                                                                                      >  8:183545/1‑136 (MQ=255)
                                        aaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGAGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCCACAAGGGC‑GCGCGCAAGGGCCATCCCAGCAGGGAG‑CGGTTGCCAGGAAccc                                                                                       >  1:22761/1‑135 (MQ=255)
                                         acacATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCa                                                                                      <  4:98156/135‑1 (MQ=255)
                                                 aaGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGtt                                                                               <  4:277374/134‑1 (MQ=255)
                                                       cAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCATCGACCCCCGTCAGCGCGCCCGGGGG‑GCGGGGCTGGGCCCGCCGCGCCGGCCG‑CGGGGGGCCGGTGCCCAGCCCCTGGGTCAgc                                                                        >  3:2494/1‑135 (MQ=255)
                                                       cAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCCGCGCTGGCAG‑CGGTTGCCACGTACCCAGCACGTTGGACCgct                                                                       >  8:279629/1‑136 (MQ=255)
                                                       cAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACAACGCCCTGGGCACCAACGACCCCCGTCATCGCGCCCTCGGC‑GGGGGCATGGGCCCCCAGCCCGGGACC‑CGCTTGGCAGGGACCCCGCACGTGGGGCCgct                                                                       >  8:170421/1‑136 (MQ=255)
                                                          cccAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCCCCGTCAGCGCGCCCGGGGC‑GCGGGCCTGGGCCAGCCGCGCGGGGAG‑CGGTTGCCAGGTACCCCGCACGTTGGGCAgcggg                                                                     >  3:108135/1‑135 (MQ=255)
                                                          cccAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGGATGGGCCCGCCGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGGTGGACCgcgggc                                                                    >  7:129330/1‑136 (MQ=255)
                                                               ttGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGAt                                                                <  6:177087/135‑1 (MQ=255)
                                                               ttGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGATGGACAGCTGGCTGAt                                                                <  4:290800/135‑1 (MQ=255)
                                                                 gctgcATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGCCCACCGTCAGCGCGGCCGGGGC‑GCGGGCCTGGGCCCGCCGCGCTGGCCG‑CGGTTTGCAGGTACCCCGCAAGTTGGACAGAGGGGTGATgg                                                              >  2:362497/1‑135 (MQ=255)
                                                                       tGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCCGGGC‑GCGGGCCTGGGCCAGCCGCGCGGGCAG‑CGGTTGCCAGGTAACCCGCACGTTGGCCAGCTGGCGGAGGGCGAGcc                                                        >  2:380814/1‑134 (MQ=255)
                                                                          ccTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCCTGGGCCAGCCGCGGTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGCCAGCTGGCTGATGgcgcgcccgcg                                                    >  8:297002/1‑131 (MQ=255)
                                                                            tGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GGGGGCATGGGCCCGCAGCGCTGGGCG‑CGGTTGCCAGGTACCCAGCACGGTGGACAGCGGGCTGGTGGCGCGCCCGCgg                                                   >  6:366324/1‑135 (MQ=255)
                                                                                       caccacGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCg                                       <  7:75769/136‑1 (MQ=255)
                                                                                       caccacGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGC‑GCGGGCAGGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGcc                                        <  5:267354/135‑1 (MQ=255)
                                                                                        accacGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGcc                                        <  1:267626/134‑1 (MQ=255)
                                                                                              cccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                                <  7:290765/136‑1 (MQ=255)
                                                                                              cccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGt                                 <  6:126880/135‑1 (MQ=255)
                                                                                              cccTGGGCACCATCGACCACCGGCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                                <  7:382971/136‑1 (MQ=255)
                                                                                               ccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                                <  4:184574/135‑1 (MQ=255)
                                                                                                  gggCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCCCGGTGGACCGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGGCGCTCAGTTTCa                             >  2:431098/1‑135 (MQ=255)
                                                                                                   ggCACCCTCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAg                            <  5:123477/135‑1 (MQ=255)
                                                                                                   ggCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCTATCAGTTGCAg                            <  4:108135/135‑1 (MQ=255)
                                                                                                   ggCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  8:351344/136‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  3:409146/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  3:439690/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                           <  5:423450/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGGGCGCGGGCCGATCAGTTGCAGt                           <  4:72894/135‑1 (MQ=255)
                                                                                                    gCACCATCGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCCG‑CGGTTGCCAGGTACCCAGCCCGGTGGACAGCTGGCTGATGGCGCGCCCTCGGGTGCGCGGGCCGCTCAGTTGCAGTg                          >  8:41551/1‑136 (MQ=255)
                                                                                                      aCCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTTCCCAGCACGTTGGACAGCTGGCTGATGGCGCGCACTCGGGTGCGCGGGCCGCTCAGTTGCAGTGc                         >  2:329457/1‑135 (MQ=255)
                                                                                                           cGACCACCGTCAGCGCGCCCGGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCCG‑CGGTTTGCAGGTACACAGCACCTTGGACACCTGGCTGATGGCGAGCACTCGGGGGCGCGGGCCGATCAGTTGCAGTGcctggt                    >  4:23377/1‑135 (MQ=255)
                                                                                                            gACCACCGTCAGCGCGCCCGGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCCG‑CGGTTGCCAGGTACCCAGCACGTTGGTCAGCTGGCTGCTGGCGCGCCATCGGGTGCGCGGGCCGCTCAGTTTCAGTGcctggtc                   >  4:298713/1‑135 (MQ=255)
                                                                                                                 ccGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  7:294484/136‑1 (MQ=255)
                                                                                                                  cggcAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGGTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  1:325019/132‑1 (MQ=255)
                                                                                                                  cGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  3:189850/135‑1 (MQ=255)
                                                                                                                  cGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc             <  1:345262/135‑1 (MQ=255)
                                                                                                                     cAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATa         >  7:336684/1‑136 (MQ=255)
                                                                                                                            cggtGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGGGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCAt  <  8:400541/133‑1 (MQ=255)
                                                                                                                            cccGGGGC‑GCGGGCAGGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCa   <  4:28167/135‑1 (MQ=255)
                                                                                                                               |                                                                                                                                      
CACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGC‑GCGGGCATGGGCCAGCAGCGCTGGCAG‑CGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGACCAT  >  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: ‑