Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,348,580 G→A 5.6% G45G (GGG→GGA)  cysT_1 → Sulfate transport system permease protein CysT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,348,5800GA 5.6% 218.3 / 2.6 71G45G (GGG→GGA) cysT_1Sulfate transport system permease protein CysT
Reads supporting (aligned to +/- strand):  ref base G (28/39);  new base A (2/2);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTT  >  KT2440/4348447‑4348714
                                                                                                                                     |                                                                                                                                      
tGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGc                                                                                                                                       >  4:36288/1‑135 (MQ=255)
 gCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGcc                                                                                                                                      <  4:12495/135‑1 (MQ=255)
        gATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCggtgggt                                                                                                                               >  4:52660/1‑135 (MQ=255)
          tGCCAGAGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCAACCTGCTACAGAGGGGCATCGCCACACCCGGCGCCTGGTGGCACCCGCCCCCCCGCTCATGGCAGCGCGGGCCGGGGGGTGc                                                                                                                             >  8:124232/1‑135 (MQ=255)
                  aCCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGgtggtg                                                                                                                      >  8:269064/1‑134 (MQ=255)
                  aCCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGgtggtg                                                                                                                      >  3:159189/1‑134 (MQ=255)
                   ccTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGc                                                                                                                    <  6:93517/135‑1 (MQ=255)
                   ccTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGc                                                                                                                    >  1:233145/1‑135 (MQ=255)
                        gcgaCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGAGCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCGCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTg                                                                                                                <  1:196807/134‑1 (MQ=255)
                               cTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCCGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCCCTGCCGCTGGt                                                                                                        >  3:241522/1‑135 (MQ=255)
                                   cTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGt                                                                                                        <  7:269064/131‑1 (MQ=255)
                                         gTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCg                                                                                              >  4:146851/1‑135 (MQ=255)
                                              aCTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGTGCCGTTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGt                                                                                         <  5:246665/135‑1 (MQ=255)
                                              aCTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGt                                                                                         <  7:21916/135‑1 (MQ=255)
                                               cTGGCCAGCCTGACCACCCTGAGCCTGCTGATCATCGGCACACCCGTCGCCGGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTCCCACCGACGGTg                                                                                        >  4:309391/1‑135 (MQ=255)
                                                     aGCCTGACCACCCTGATGCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGt                                                                                  <  6:68287/135‑1 (MQ=255)
                                                     aGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGt                                                                                  >  7:141645/1‑135 (MQ=255)
                                                     aGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGt                                                                                  >  4:230879/1‑135 (MQ=255)
                                                     aGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGt                                                                                  >  3:283468/1‑135 (MQ=255)
                                                     aGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGt                                                                                  >  1:184231/1‑135 (MQ=255)
                                                        cTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGtttt                                                                               <  4:146356/135‑1 (MQ=255)
                                                        cTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGtttt                                                                               <  7:303066/135‑1 (MQ=255)
                                                          gACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTa                                                                             >  3:150600/1‑135 (MQ=255)
                                                               cccTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGGTTTTACCTGt                                                                        >  8:295815/1‑135 (MQ=255)
                                                                ccTGATCCTGCTGATCATCGGCACACCGGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGtt                                                                       <  3:230879/135‑1 (MQ=255)
                                                                ccTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGtt                                                                       <  8:141645/135‑1 (MQ=255)
                                                                ccTGATCCTGCTGATCATCGGCACACCCGGCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGtt                                                                       <  1:321117/135‑1 (MQ=255)
                                                                       ctgctgATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGcc                                                                                                                                      <  8:205112/65‑1 (MQ=255)
                                                                       ctgctgATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGcc                                                                                                                                      >  7:205112/1‑65 (MQ=255)
                                                                            gATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGtt                                                                       >  1:119629/1‑123 (MQ=255)
                                                                            gATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGtt                                                                       <  2:119629/123‑1 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTTCTGCCACCGACGGt                                                                                         >  6:261764/1‑104 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTTCTGCCACCGACGGt                                                                                         <  5:261764/104‑1 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGc                                                          >  8:24720/1‑135 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGc                                                          >  1:24731/1‑135 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGc                                                          >  1:86954/1‑135 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGc                                                          >  4:62402/1‑135 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTAACTGTTGATCGCCCTTGGc                                                          >  3:224865/1‑135 (MQ=255)
                                                                             atcatcGGCACACCCGTCGCCTGGGGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGc                                                          <  7:290578/135‑1 (MQ=255)
                                                                              tcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGc                                                          <  4:164880/134‑1 (MQ=255)
                                                                              tcatcGGCACACCCGGCGCCTGGTGGCTGGCGCGCACGCGCTCTTGGCGGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTAACTGTTGATCGCCCTTGGc                                                          <  4:224865/134‑1 (MQ=255)
                                                                                  cGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCa                                                     >  3:217327/1‑135 (MQ=255)
                                                                                  cGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCa                                                     <  4:241522/135‑1 (MQ=255)
                                                                                     cacaCCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGACCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACgg                                                  >  5:26782/1‑135 (MQ=255)
                                                                                     cacaCCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGACCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCGTGGCCCGCACgg                                                  >  5:22517/1‑135 (MQ=255)
                                                                                                ccTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGACCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTc                                       <  6:26782/135‑1 (MQ=255)
                                                                                                ccTGGGGGCTCGCGCGCACGCGCTCTGGGCTGCGCGGACCGGGGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTc                                       <  6:22517/135‑1 (MQ=255)
                                                                                                     tggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCAACGACGGGGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCTGGcc                                  >  1:261585/1‑135 (MQ=255)
                                                                                                      ggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCa                                 <  4:255868/135‑1 (MQ=255)
                                                                                                       gCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGcc                                                                >  6:249451/1‑103 (MQ=255)
                                                                                                       gCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGcc                                                                <  5:249451/103‑1 (MQ=255)
                                                                                                       gCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCAc                                <  6:117333/135‑1 (MQ=255)
                                                                                                       gCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGGGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCAc                                <  1:54095/135‑1 (MQ=255)
                                                                                                            cgcgcACGCGCTCTTGGGTGCGCGGGCCGGTGGGGGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAgg                           <  3:12297/135‑1 (MQ=255)
                                                                                                            cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAgg                           <  2:233145/135‑1 (MQ=255)
                                                                                                            cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAgg                           <  3:142052/135‑1 (MQ=255)
                                                                                                            cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAgg                           <  7:265353/135‑1 (MQ=255)
                                                                                                            cgcgcACGCGCGCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGGGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCGGGCTTGGTCAGGCCACCCAgg                           <  3:52660/135‑1 (MQ=255)
                                                                                                                 aggcgcTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                      <  4:283699/133‑1 (MQ=255)
                                                                                                                 aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                      <  6:137193/135‑1 (MQ=255)
                                                                                                                 aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                      <  4:217327/135‑1 (MQ=255)
                                                                                                                      ctctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCTCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCt                 <  7:146518/135‑1 (MQ=255)
                                                                                                                      ctctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCt                 <  2:184231/135‑1 (MQ=255)
                                                                                                                        ctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg               <  7:24720/135‑1 (MQ=255)
                                                                                                                        ctTGGCTGCGCGGGCCGGTGGGGGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg               <  3:65502/135‑1 (MQ=255)
                                                                                                                          tGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg               <  3:146851/133‑1 (MQ=255)
                                                                                                                          tGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg               <  3:36288/133‑1 (MQ=255)
                                                                                                                          tGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg               <  3:158821/133‑1 (MQ=255)
                                                                                                                                  cgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCt                 >  4:15643/1‑123 (MQ=255)
                                                                                                                                  cgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCt                 <  3:15643/123‑1 (MQ=255)
                                                                                                                                  cgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTCGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggt     >  3:271768/1‑135 (MQ=255)
                                                                                                                                     gCCGGTGGGGGCGGTGGTGGCACTGCGGCGGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGtt  <  1:321373/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTT  >  KT2440/4348447‑4348714

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: ‑