Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,135,590 A→C 27.0% Q32P (CAG→CCG)  tgnC → (Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,135,5900AC27.0% 51.4 / 5.6 61Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base C (14/4);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 5.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCT  >  ACB143/4135455‑4135713
                                                                                                                                       |                                                                                                                           
cacgCCCGCACCCTCACCGATTCCGCCAAGGGGTTTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGTTGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                             <  4:167990/133‑1 (MQ=255)
cAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                             <  3:336659/136‑1 (MQ=255)
cAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                             <  8:11053/136‑1 (MQ=255)
     ccccACCCTCACCTAATCCGCCATGGGGATTTTTGCTTGACCTTCGTGGTTTGCCAGATGTGCATCGACGGCCAATGGCGCTATGCCCAGAGCTCCAAGACCTTCGACAGCGTCAACCCGGCCACCTCCCAGGGCt                                                                                                                        <  6:34427/133‑1 (MQ=255)
      cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGccag                                                                                                                       >  6:299212/1‑134 (MQ=255)
      cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGccag                                                                                                                       >  7:126501/1‑134 (MQ=255)
      cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCTg                                                                                                                       >  1:156599/1‑136 (MQ=255)
          cccTCACCGAATCCGCCAAGGGTTTTGTTGCATGACCCTCGTTCGTTTCCAGTTGTGCATCGACGGCCAATGGCGCTATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGgcg                                                                                                                   <  3:207714/136‑1 (MQ=255)
           ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGgcgc                                                                                                                  >  1:155813/1‑136 (MQ=255)
           ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCAGGgcgc                                                                                                                  >  7:247496/1‑136 (MQ=255)
                 cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCACCTGc                                                                                                            >  2:247952/1‑136 (MQ=255)
                 cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCAGGGCGCAACTGc                                                                                                            >  2:284888/1‑136 (MQ=255)
                  gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGcc                                                                                                           >  8:184475/1‑136 (MQ=255)
                        gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACAGcccccccc                                                                                                     >  6:299455/1‑131 (MQ=255)
                        gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGc                                                                                                     >  3:59288/1‑136 (MQ=255)
                           aaGGGGTTTGTTGCATGACCCTCGTTCGTTTCCAATTGTCCAGCGCCGGCCAAGGGCGCGAGGCCCAGAGCGGAAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGa                                                                                                  <  3:188147/136‑1 (MQ=255)
                              gggATTTTTGCATGACCCGCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                               <  5:385034/136‑1 (MQ=255)
                              gggATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                               <  7:94637/136‑1 (MQ=255)
                                           gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGCTGAAGCCGATGTGGAg                                                                                  >  4:132545/1‑136 (MQ=255)
                                           gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCAGGGCGCAACTGCCCGACCCCGAAGAAGCCGAAGTGGAg                                                                                  >  8:351262/1‑136 (MQ=255)
                                             cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCt                                                                                >  5:378271/1‑136 (MQ=255)
                                             cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGCCACAGCCCCCAGCCGATGAAGCCGATGTGGAg                                                                                  >  1:367163/1‑134 (MQ=255)
                                             cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCAGGGCGCAACAGCCCGCCGCCGCTGAAGCCGAAGTGGAGCt                                                                                >  2:320249/1‑136 (MQ=255)
                                             cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCACCCCTGGGGGCAACAGCCCGCCGCCCAACAACCCCATGTGGAGCt                                                                                >  1:144641/1‑136 (MQ=255)
                                                tcgttcttTTCCAGTTGTGAATCGACGCCCAAGGGCGCGATGCCCATAGCGCAAACACATTCGACCGCCACACCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGc                                                                             <  4:363135/136‑1 (MQ=255)
                                                tcgttcgtTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCAGGGCGCCACTGCCCGCCGCCGCTGAAGCCCAGGGGGAGCTGGc                                                                             >  7:135287/1‑136 (MQ=255)
                                                tcgttcgtTTCCAGATGAGCAACGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCAGGGCGCACCAGACCGCCCCCCATGAAGACGATGTGGAGCTGGc                                                                             >  7:79077/1‑136 (MQ=255)
                                                  gttcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCAGGGCGCCACAGCCCCCCGCCGAAGAACCCGAAGTGGACCAGGccc                                                                           >  5:127634/1‑135 (MQ=255)
                                                    tcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCAGGGCGCCACTGCCCCCCCCCCCTGAAGCCGATGTGGAGCTGGCCGtt                                                                         >  1:228458/1‑136 (MQ=255)
                                                    tcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCAGGGCGCCACAGCCCCCCGCCGATGAAGCCGATGTGGAGCTGGCCGt                                                                          >  1:27743/1‑135 (MQ=255)
                                                    tcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCACGCCGGGGCGCAACTGCCCGACCCCGATGAAGCCGAAGTGGAGCTGGccct                                                                          >  6:127648/1‑133 (MQ=255)
                                                     cttttCCAGATGTGCATCTACGGCCAATGGCGCGAGGCCCAGCTCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCGTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGtt                                                                         <  4:344312/133‑1 (MQ=255)
                                                     cgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCTGGGCGCCACTGCCCGACGCCGAAGAAGCCGATGTGGAGCTGGCCGtt                                                                         >  3:70151/1‑135 (MQ=255)
                                                      gttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTc                                                                        <  6:103628/135‑1 (MQ=255)
                                                         tCCAGTTGTGCTGCTCCGGCCATTGGCGCGATGCCCAGAGCGGCATGACCTTCGACACCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAgg                                                                     <  2:367163/135‑1 (MQ=255)
                                                          ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCAGGGCGCCACTGCCCCCCGCCCCTGAAGCCCATGTGGAGCTGGCCGTTCAGgcc                                                                   >  6:50383/1‑136 (MQ=255)
                                                              aTGTGCATCGACGCCCAAGGGCGCGATGCCCAGAGCGCCAACCCCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGccc                                                               <  4:352064/136‑1 (MQ=255)
                                                               tgtgCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGccc                                                               <  2:162562/135‑1 (MQ=255)
                                                                 gccATCCACGCCCAAGGGCGCGAGGCCCCGAGCGGCACCCCCTTCCCCCGCCTCACCCCGGCCCCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAg                                                             <  4:238904/133‑1 (MQ=255)
                                                                    aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGGCAATGCCCGCCGCCGCTGAAGCCGAGGTGGAGCTGGCCGTGCCGGCCCCCCcgcggg                                                         >  2:393763/1‑134 (MQ=255)
                                                                      cGACGGCCAATGGCGCAATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCGTGcc                                                       >  3:90619/1‑136 (MQ=255)
                                                                             cAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                <  3:30505/136‑1 (MQ=255)
                                                                              aaTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGGACCTGCCCGACGCCGCTGAAGCCGCTGTGGAGCTGGCCGTTCCGGCCGCGCAGCGGGCCTTGGGCAg                                               >  8:375304/1‑136 (MQ=255)
                                                                                        tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCGTGCCTTCGCCAGCAAGGCGTgg                                     >  2:236114/1‑136 (MQ=255)
                                                                                        ggCCCAGCGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                      <  6:127634/134‑1 (MQ=255)
                                                                                           ccAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                  <  4:112552/136‑1 (MQ=255)
                                                                                           cacgagCGCCACCACCTGCGACACCCGAAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                  <  6:232523/133‑1 (MQ=255)
                                                                                             agagCGGCAAGACCTTCGACCGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                <  2:259862/136‑1 (MQ=255)
                                                                                             agagCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                <  2:155813/136‑1 (MQ=255)
                                                                                             agacCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                <  2:209483/136‑1 (MQ=255)
                                                                                                     aaGACCTTCGACCGCCTCACCCCGCCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCg                        <  1:256596/136‑1 (MQ=255)
                                                                                                       gACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGATGTGGTGCTGGCCGTTCAGGCCGCCCCGCGTGCCTTCGCCAGCAACGCGTGGCGCAGCATCACCGc                       >  8:250628/1‑135 (MQ=255)
                                                                                                       gACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCGTCGACAGCACGGCGGGGCGCCGCATCAcccc                       >  6:16272/1‑133 (MQ=255)
                                                                                                       gACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCCATGAAGCCGATGTGGAGCTGGCCGTTAAGGCCGCCCAGCGTGCCTTCGCCAGCAAGGCGTGGCGCCGCATAAcccc                       >  8:338589/1‑133 (MQ=255)
                                                                                                       gACCTTCGACAGCCTCAACCCGGCCACCCCCCAAGCCTGGGCGCCACTGCCCCACGCCCATGAAGACGATGTGGAGCTAGCCCTTAAGGCCGCCCACCGTGCCATAGACAGCCAGGCGTTGCGCCGCATAACCGc                       >  6:336568/1‑135 (MQ=255)
                                                                                                                 cgCCACACCCCGGCCCCCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg            <  7:342060/135‑1 (MQ=255)
                                                                                                                 aGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg            <  8:190487/136‑1 (MQ=255)
                                                                                                                 aGCCTCAACCCGCCCACCGCCATGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg            <  1:320249/136‑1 (MQ=255)
                                                                                                                   ccTCACCCCGGCCACCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGc          <  5:305208/136‑1 (MQ=255)
                                                                                                                      cAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCCGCATCACCGCTACCCCCCCCGGCAAg       >  6:313891/1‑136 (MQ=255)
                                                                                                                        aCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGCCGCCGATGCAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGGGCCTTCGACAGCAAGGCGTGGCGCAGCATCCCCGCTACCGCCCGCGGCCAgct     >  2:288477/1‑136 (MQ=255)
                                                                                                                        aCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTTCCTTCGACAGCAAGGCGTGGCGCAGCATAACCGCTACCGCCCGCGGCAAg       >  2:121366/1‑134 (MQ=255)
                                                                                                                           cgcccACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct  <  6:250723/133‑1 (MQ=255)
                                                                                                                           cGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct  <  2:148565/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                           
CAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCT  >  ACB143/4135455‑4135713

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