Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,894,256 C→G 5.8% L102L (CTG→CTC)  sbp_2 ← Sulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,894,2560CG 5.8% 244.7 / 6.4 69L102L (CTG→CTC) sbp_2Sulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base G (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

ATGACCGACACGTCCTTCTTGATCAGGTCGCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCGACCCCCCATGGGATTGCTGGATTTTCACCT  >  ACB143/5894122‑5894389
                                                                                                                                      |                                                                                                                                     
atgaCCGACACGTCCTTCTTGATCAGGTCGCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCa                                                                                                                                      <  6:326719/136‑1 (MQ=255)
      gACACGTCCTTCTTGATCAGGTCGCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGGCGAGGTGTACGGGGGGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCt                                                                                                                                <  7:321461/136‑1 (MQ=255)
             ccTTCTTGATCAGGTCGCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAg                                                                                                                         >  1:155545/1‑136 (MQ=255)
                 cttGATCAGGTCGCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCtt                                                                                                                      <  4:378921/135‑1 (MQ=255)
                    gATCAGGTCGCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTgg                                                                                                                    <  1:29589/134‑1 (MQ=255)
                             gCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGGAg                                                                                                                                     >  1:278192/1‑108 (MQ=255)
                             gCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGGAg                                                                                                                                     >  1:283846/1‑108 (MQ=255)
                             gCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGGAg                                                                                                                                     <  2:283846/108‑1 (MQ=255)
                             gCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGGAg                                                                                                                                     <  2:278192/108‑1 (MQ=255)
                             gCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTtcgtc                                                                                                          >  3:156831/1‑135 (MQ=255)
                             gCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTtcgtc                                                                                                          >  3:242205/1‑135 (MQ=255)
                               cccAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGAt                                                                                                       >  6:180423/1‑136 (MQ=255)
                                ccAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATg                                                                                                      <  2:37290/136‑1 (MQ=255)
                                 cAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGt                                                                                                     <  4:366691/136‑1 (MQ=255)
                                 cAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGt                                                                                                     <  1:79901/136‑1 (MQ=255)
                                 cAGTCCTTGATGCCTTTCGGGGTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGTCAGTGTCTTGCCAAGCTTGGCCACTTCGTCGATGt                                                                                                     <  2:262242/136‑1 (MQ=255)
                                      cTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCg                                                                                                <  4:242205/136‑1 (MQ=255)
                                            gCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCGAATGTTGGCCTCTTCGTCGATGTCACCggccagg                                                                                          >  7:142096/1‑136 (MQ=255)
                                            gCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAATCTTGGCCACTTCGTCGATGTCACCggccagg                                                                                          >  5:14905/1‑136 (MQ=255)
                                            gCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccagg                                                                                          >  4:394697/1‑136 (MQ=255)
                                            gCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccagg                                                                                          >  5:386035/1‑136 (MQ=255)
                                                 tCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccagggccag                                                                                     >  4:242462/1‑136 (MQ=255)
                                                 tCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccagggccag                                                                                     >  1:214615/1‑136 (MQ=255)
                                                 tCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTACCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccagggcc                                                                                       <  6:196096/134‑1 (MQ=255)
                                                  cGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccagggccag                                                                                     <  1:299654/135‑1 (MQ=255)
                                                   gggTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGccagggccaggg                                                                                   >  8:327395/1‑136 (MQ=255)
                                                                 cACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGgcg                                                                     >  6:16019/1‑136 (MQ=255)
                                                                 cACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGgcg                                                                     >  4:156596/1‑136 (MQ=255)
                                                                    cAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGGCAGTCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGgcgcgc                                                                  <  2:295684/136‑1 (MQ=255)
                                                                       gAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCCgg                                                               >  6:218337/1‑136 (MQ=255)
                                                                        atcacGCTGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGc                                                              <  2:350909/134‑1 (MQ=255)
                                                                            cGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTc                                                          <  3:49139/136‑1 (MQ=255)
                                                                            cGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTc                                                          <  1:222991/136‑1 (MQ=255)
                                                                            cGATGGGCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTc                                                          <  3:54500/136‑1 (MQ=255)
                                                                              aTGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGa                                                        <  3:87676/136‑1 (MQ=255)
                                                                                 gTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATc                                                      >  6:64830/1‑135 (MQ=255)
                                                                                 gTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCa                                                     >  7:122807/1‑136 (MQ=255)
                                                                                    gAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGTCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCg                                                  >  2:300705/1‑136 (MQ=255)
                                                                                    gAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCAcc                                                   <  5:64830/135‑1 (MQ=255)
                                                                                    gAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCAcc                                                   <  7:176772/135‑1 (MQ=255)
                                                                                    gAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCCATCACCg                                                  >  2:123622/1‑136 (MQ=255)
                                                                                     aGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCCGGCCGTCGATCACCGc                                                 >  3:119951/1‑136 (MQ=255)
                                                                                       gtggACGGTGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGTTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACTTCGGCGCGCAGGCCGTCGATCACCGCAc                                               <  6:257096/136‑1 (MQ=255)
                                                                                          tACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACggg                                            >  1:220316/1‑136 (MQ=255)
                                                                                               ggTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGc                                       >  5:67861/1‑136 (MQ=255)
                                                                                                  gCTGGCGGCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTg                                    <  1:367465/136‑1 (MQ=255)
                                                                                                    tGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGcc                                  >  4:315939/1‑136 (MQ=255)
                                                                                                    gggcggcCGGCCGGCGGGTCGGCGAGCTCTCGGGCAGGGGCGTGCCCAGCTTGGCCCCTTCGTCGAGGGCAGCGTCCCGGGCGAGGGTCACCACGTCGCCGCGCCGGCCGTCGATCACCGCACGGGCCGGCTGGcc                                  <  2:92697/135‑1 (MQ=255)
                                                                                                     ggCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGt                                                                                                     >  4:115922/1‑68 (MQ=255)
                                                                                                     ggCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGt                                                                                                     <  3:115922/68‑1 (MQ=255)
                                                                                                     ggCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccag                                                                                           >  8:18891/1‑78 (MQ=255)
                                                                                                     ggCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccag                                                                                           <  7:18891/78‑1 (MQ=255)
                                                                                                     ggCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCCGGCCGTCGCTCACCGCACCGGCCGGCTTGccc                                 >  1:203970/1‑136 (MQ=255)
                                                                                                     ggCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGAACACCGCACGGGCCTGCTTGccc                                 >  5:153669/1‑136 (MQ=255)
                                                                                                      gCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGcc                                  <  2:220316/134‑1 (MQ=255)
                                                                                                      gCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGcc                                  <  6:327539/134‑1 (MQ=255)
                                                                                                      gCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCg                                <  4:274384/136‑1 (MQ=255)
                                                                                                       cGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCg                                                                                                <  2:80136/71‑1 (MQ=255)
                                                                                                       cGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCg                                                                                                >  1:80136/1‑71 (MQ=255)
                                                                                                       cGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCGa                               >  7:360727/1‑136 (MQ=255)
                                                                                                              caggcGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCGGCTTGCCCGcccccc                          >  3:153631/1‑134 (MQ=255)
                                                                                                              caggcGGGACTGCCAGATCTCCGGCAGGGACTTGCCAAGCTAGGCCACTTCGTCGATGTCACCCGCCAGCGCCACGGTCACCACGTAGGCCCGCAGGCCGTCCATCACACCACCCGCCAGCATGCCCCAcccccc                         >  2:168498/1‑135 (MQ=255)
                                                                                                                ggcGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccag                                                                                           >  6:287488/1‑67 (MQ=255)
                                                                                                                ggcGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCggccag                                                                                           <  5:287488/67‑1 (MQ=255)
                                                                                                                     gTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCGACCCCCCATGGGAtt                 <  7:22529/136‑1 (MQ=255)
                                                                                                                     gTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCGACCCCCCATGGGAtt                 <  7:686/136‑1 (MQ=255)
                                                                                                                           cAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCGACCCCCCATGGGATTGCTGGa           <  6:12216/136‑1 (MQ=255)
                                                                                                                                  cGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGCCGCGCACGCCGTCCATCACCGCACCGGCCTGCTTGCCCGACCCCCCATGGGATTGCTGGATTTTCAc    >  1:270217/1‑136 (MQ=255)
                                                                                                                                    ggCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCGACCCCCCATGGGATTGCTGGATTTTCACCt  <  2:393686/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
ATGACCGACACGTCCTTCTTGATCAGGTCGCCCCAGTCCTTGATGCCTTTCGGGTTGCCCTTGCGCACCAGGAACACGATGGTCGAGGTGTACGGGGTGCTGGCGTCAGGCAGGCGGGTCTGCCAGTTCTCGGGCAGGGTCTTGCCAAGCTTGGCCACTTCGTCGATGTCACCGGCCAGGGCCAGGGTCACCACGTCGGCGCGCAGGCCGTCGATCACCGCACGGGCCTGCTTGCCCGACCCCCCATGGGATTGCTGGATTTTCACCT  >  ACB143/5894122‑5894389

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