Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 381,772 T→G 23.5% L198R (CTG→CGG)  soxB → Sarcosine oxidase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122381,7720TG23.5% 86.3 / 2.9 60L198R (CTG→CGG) soxBSarcosine oxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base G (12/3);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 3.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGC  >  ACB122/381642‑381906
                                                                                                                                  |                                                                                                                                        
gCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCAGGGc                                                                                                                                      >  3:97252/1‑135 (MQ=255)
  cGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGGCCGTCACGATGCCGTGGCCTGGGGGGTGGCCCCCGGAGCCGCCCCCCCGGGCgg                                                                                                                                    >  5:356063/1‑134 (MQ=255)
  cGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGCTGCCGTGGCCTGGGGCTTGGCCCGCGCTGCCCCCGCCCGGGGCgg                                                                                                                                    >  1:129661/1‑134 (MQ=255)
  cGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACCCCCGGGGCGt                                                                                                                                    >  2:146525/1‑135 (MQ=255)
       tCCCGTACCTGGCCTGCTCGATGAACACCCGCTCCCCGATCATCGGCGCCCCCGTACCGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTGCGCCCGCGCTGCCGCCGCCCGGGGCGTCGAccg                                                                                                                              >  7:70609/1‑135 (MQ=255)
       tCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCCCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGCTGCCGTGGCCTGGGGCGTCGCCCGCGCTGCCGCCGCCCGGGGCGGCGGcc                                                                                                                               >  6:114910/1‑135 (MQ=255)
       tCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGTCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGccc                                                                                                                               >  4:311894/1‑135 (MQ=255)
           gTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGAt                                                                                                                          <  2:308844/135‑1 (MQ=255)
                          gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGt                                                                                                           >  8:110951/1‑135 (MQ=255)
                          gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCAGCCGACGCCCTGGGCGTCCACC‑TGATCCAGCAAACCGAAg                                                                                                            >  1:6411/1‑134 (MQ=255)
                          gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCCCTGCCGCCGCCCTGGGCGTAGCCC‑GGATGCCGCAAACCgcatgg                                                                                                          >  8:77033/1‑131 (MQ=255)
                          gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTAGCCCCCGCTGCCCACCCCCAGGGGCGCGACC‑GGGACCAGCCAAccccaa                                                                                                            >  4:338472/1‑130 (MQ=255)
                                    cGCTACCCGATCCTCGGCGCCACCGTCCAGCGCCGTGGTGGGGTTGCCCGTCACGATGCCGGGGCCTGGGGCTGCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCtgc                                                                                                 <  6:296444/135‑3 (MQ=255)
                                    cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCCCGCTGCCGACGCCCTGGGCGTCGACC‑AGATCCAGCAAACCGAAGTGATCGGcagac                                                                                                >  7:322760/1‑132 (MQ=255)
                                    cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTAGCCCGCGCTGCCCCCCCCCAGGGCGGGGCCC‑CGATACCGCAAAACCAAAAGATCGGcata                                                                                                 >  3:55364/1‑132 (MQ=255)
                                    cGCTACCCGATCCTCGGCGCAACCGGACCGCGCGGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTc                                                                                                 <  3:49521/135‑1 (MQ=255)
                                        aCCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTAGCCCGCGCTGCCGACCCCCAGGGCGGGGGCC‑TGGAACAGCAAAACGCAATGATGGGCGTACGc                                                                                              >  2:247274/1‑134 (MQ=255)
                                                    gcgcCACCGTCCAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACgg                                                                                   <  5:85907/133‑1 (MQ=255)
                                                    gcgcCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCGGGGCGGCGCCC‑GGCTCCCGCAAACCGAAGTGCTCGGCTTCCGCCAGGCAAACGgcgc                                                                                >  8:66919/1‑136 (MQ=255)
                                                     cgcCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGGGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGgcgc                                                                                <  3:176702/135‑1 (MQ=255)
                                                      gcCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCGGGGCGGCGGCC‑GGATGCAGCCAACCGAAGTGATCCGCTTCCGCCAGGAAAACGGCGCgg                                                                              >  8:19792/1‑136 (MQ=255)
                                                       ccACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCGGGGCGGCGCCC‑GGCTCCCGCAAACCCCAGTGATCGGCGTCCGCCAGGACACCGGCCCggg                                                                             >  7:303465/1‑135 (MQ=255)
                                                        cACCGTACAGCGCCGTGGGGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCTACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGt                                                                             <  1:247274/135‑1 (MQ=255)
                                                         aCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTc                                                                            <  4:142351/135‑1 (MQ=255)
                                                         aCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGAGGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCa                                                                           <  7:332458/136‑1 (MQ=255)
                                                           cGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGGTTCGCCCGCGCTGCCGCCCCCCGGGGCGGCGGCC‑GGGTCCCGCCAACCCAAGTGATCCGGGTGCGCCACGAAAACGGCGCggggat                                                                          >  1:143740/1‑131 (MQ=255)
                                                                  cgcCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCTGGGGGTCGCCC‑TGGTCCAGCCAACCGAAGTGATCGGCTTGCGCCAGGGAAAACGCGCGGGCATGCGtgt‑g                                                                  >  1:47011/1‑135 (MQ=255)
                                                                   gcCGTGGTGGCGTTGCCCGGCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑gg                                                                 <  4:246110/135‑1 (MQ=255)
                                                                     cGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCCGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGaaa                                                              >  7:309554/1‑136 (MQ=255)
                                                                        ggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCGGGGCGGCGGCC‑GGGTGCCGCCAACCGAAGTGGTCGGCGTGCGCCAGGAAAACGGCGCGGGCATCGGTGT‑GGGAAcc                                                            >  8:269855/1‑135 (MQ=255)
                                                                         gtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACca                                                           <  3:144495/135‑1 (MQ=255)
                                                                         gtggCGTTGCCCGGCACGATGCCGTGGCCTGGGGCTTCGCCCGGGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAACCCGAAGTGAGCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACca                                                           <  3:138293/135‑1 (MQ=255)
                                                                         ggggggTTGCCCGTCACGTTGCGGGGGCCTGGGGCTGCGCCCGCGGTGCCTACGCCCTGGGCGTCGACC‑TCACCCCCAAAACCTAAGTTATCGGCTTCCGAAAGTAAAACGGCGCGGTCATCGGTGT‑GGAAACca                                                           <  2:390335/130‑1 (MQ=255)
                                                                          tggCGTTGCCCGTCACGATGCCGTGGCCTGGGGGTTGGCCCGCGCTTCCGACCCCCGGGGCGTGGCCCGTGGACCACCAAACC‑AATTGATAAGCTTCCGGCACGAAAACGGCGCGGTGACTGGGGTGGGAAACca                                                           >  6:335628/1‑135 (MQ=255)
                                                                          tggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGGCGCCCGGGGCGGGGCCC‑CGCTCCACCCAACCGCACTGTTCGGCGTACGCCAGGGAAAACGCGCGGGGATCGGGGG‑GGGAACca                                                           >  4:391371/1‑134 (MQ=255)
                                                                          tggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCGGGGCGTGGACC‑GGGTCCCCCCAACCGAAGTGATCGGCGTCCGCCAGGAAAACGGCGCGGTCCTCGGTGT‑GGGAACca                                                           >  5:109486/1‑134 (MQ=255)
                                                                          tggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGGCC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGGAAACGGCGCGGTCATCGGTGT‑GGaaaacca                                                          >  8:249980/1‑131 (MQ=255)
                                                                          gggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCGGCGCTGCCGACCCCCTAGGCTTCTACC‑TGATCCAGCAAACCGAAGTGAGGGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACca                                                           <  3:282297/133‑1 (MQ=255)
                                                                               ttGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGCCC‑TGCTCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGGAACCAACAAAg                                                     >  6:194994/1‑135 (MQ=255)
                                                                                      tCACGATGCCGTGGCCTGGGGCTTCGCCAGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATc                                             >  8:381665/1‑136 (MQ=255)
                                                                                         cGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCTGGGGGTCGCCC‑TGATCCAGGAAACCGAAGTGCTAGGCGTCCGCAAGGAAAACGGCGCGGTCCTCGGTGT‑GGAAACCCACAAAGGCTTCATCgg                                           >  1:17822/1‑135 (MQ=255)
                                                                                         cGATGCCGTGGCCTGGGGCTTAGCCCGCGCTGCCGACGCCCTGGGCGTCGCCC‑TGCTCCCGGAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCCACAAAGGCGTCAACGgc                                          >  8:376293/1‑136 (MQ=255)
                                                                                               cGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCGACC‑GGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAAAGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATAGGCGCCAAg                                    >  7:51231/1‑136 (MQ=255)
                                                                                                    ccTGGGGCTTCGCCCGCGCTGCCGACCCCCTGGGCGTCGACC‑AGATCCAGCAAACCGAAGTGATCGGCATACGCAAGGAAAACGGCGCGGTCATACGTGT‑GGAAACCAACAAAGGCGTCAACGGCGCCAcgcgcg                                >  3:118827/1‑135 (MQ=255)
                                                                                                         ggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGCCC‑TGATGCAGCCAACCGCAGTGATCGGCGTACGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGc                           >  2:84520/1‑135 (MQ=255)
                                                                                                         ggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑GGATCCAGCAAACCCAAGTGATCGGCGTCCGCAAGGAAAACGGCGCGGTCCTCGGGGT‑GGAAACCAACAAAGGCTTCATCGGCGCCACGCGCCTCGGc                           >  1:248776/1‑135 (MQ=255)
                                                                                                              cGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTc                     <  7:269855/136‑1 (MQ=255)
                                                                                                               gCCCGCGGCCCCGCCGCCCGGCGCGGCGACC‑TGAGCCACCAACCCCCAGAGATCGGCCGCCCCCCACAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTc                     <  1:175169/135‑1 (MQ=255)
                                                                                                               gCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTgg                       <  1:237514/133‑1 (MQ=255)
                                                                                                               gCCCGCGCGGCCGACGCCCGGGGCGTCGACC‑TGATCCACCAAACCGAAGTGATCGGCTTCCGCACGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTc                     <  6:71792/135‑1 (MQ=255)
                                                                                                               gCCCGCGCGGCCGACGCCCGGGGCGTCGACC‑TGACCCAGCAAACCGAAGTGAGCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTc                     <  3:338472/135‑1 (MQ=255)
                                                                                                                 ccGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACaaa                                                      <  5:308674/100‑1 (MQ=255)
                                                                                                                 ccGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACaaa                                                      >  6:308674/1‑100 (MQ=255)
                                                                                                                       tGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCgg             >  7:254728/1‑135 (MQ=255)
                                                                                                                           gACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGAGCGGCTTCCGCAAGTAAAACGGCGCGGTCATGGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAAc         <  1:146525/135‑1 (MQ=255)
                                                                                                                            aCGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACt        <  4:97252/135‑1 (MQ=255)
                                                                                                                                ccTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCgg    <  4:55364/135‑1 (MQ=255)
                                                                                                                                ccTGGGCGACGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCgg    <  5:104040/135‑1 (MQ=255)
                                                                                                                                 gtGGGCGTCTACT‑TTATCCACCAACCCTCAGTTATCGGCTTCCGAAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAATGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGGCACCGCCGGTAACTCCggg   <  6:356063/134‑1 (MQ=255)
                                                                                                                                 cTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCggg   <  3:168619/135‑1 (MQ=255)
                                                                                                                                 cTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCggg   <  2:46369/135‑1 (MQ=255)
                                                                                                                                  tGGGCGTCGACC‑TGATCCTGCAAACCGAAGTGAGCGGCTTCCGCAAGTAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGc  <  1:115398/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
GCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACC‑TGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGT‑GGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGC  >  ACB122/381642‑381906

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