Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,274,993 T→A 12.7% V153E (GTA→GAA)  nsaD → 2‑hydroxychromene‑2‑carboxylate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,274,9930TA12.7% 177.1 / 5.4 68V153E (GTA→GAA) nsaD2‑hydroxychromene‑2‑carboxylate isomerase
Reads supporting (aligned to +/- strand):  ref base T (42/17);  new base A (2/7);  total (44/24)
Fisher's exact test for biased strand distribution p-value = 7.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTATCGGTGCTGTTCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGACTTTGTCGAGGAGGCGCTACGT  >  ACB122/3274861‑3275123
                                                                                                                                    |                                                                                                                                  
tGCTATCGGTGCTGTTCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTaa                                                                                                                                  <  3:159805/135‑1 (MQ=255)
  cTATCGGTGCTGTTCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAAAACTCGAGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTaaaa                                                                                                                                >  5:48462/1‑135 (MQ=255)
   tATCGGTGCTGTTCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTaaaa                                                                                                                                >  1:3703/1‑134 (MQ=255)
   tATCGGTGCTGTTCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAag                                                                                                                               >  8:143540/1‑135 (MQ=255)
            tGTTCACTGGGATTTGGGTACAGGGCCGCAACCTGAGCTACAGCGCTGTGGTGAATTAAACACTCGCGCTGGCTGGCTGCGACCCACAGTTTTTTAATGCCCTGGCCGCTGACAGCGATGTAAAAGCAGCGCTaaa                                                                                                                     <  8:73084/136‑1 (MQ=255)
               tCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGCCAGCGAAGTAAAAGCAGCGCTAAAGc                                                                                                                   >  4:42805/1‑135 (MQ=255)
               tCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCa                                                                                                                  >  7:294116/1‑136 (MQ=255)
                 aCTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGAAGCGCTAAAGCAc                                                                                                                 >  6:278111/1‑135 (MQ=255)
                      gCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGAcc                                                                                                           >  7:121039/1‑136 (MQ=255)
                      gCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGAc                                                                                                            >  3:325258/1‑135 (MQ=255)
                        tATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGAc                                                                                                            >  3:206590/1‑133 (MQ=255)
                           gggTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAgg                                                                                                       >  3:4186/1‑135 (MQ=255)
                             gTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGcc                                                                                                     >  5:8472/1‑135 (MQ=255)
                             gTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGcc                                                                                                     >  6:39873/1‑135 (MQ=255)
                             gTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGcc                                                                                                     >  1:68318/1‑135 (MQ=255)
                             gTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCg                                                                                                    >  8:68210/1‑136 (MQ=255)
                               aCAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGc                                                                                                   >  3:257400/1‑135 (MQ=255)
                                cAGCGCCGCAACCTGCGCGACAGCGCGGTGCTTAATCAAACACTCGCGCCGGCAGGCTCCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGAAAACCCAGCCAAAACGCCCGCGACCGAGGCCGCt                                                                                                  <  4:156929/135‑1 (MQ=255)
                                  gcgcCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCtgt                                                                                                >  4:340729/1‑135 (MQ=255)
                                  acgccgcAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCtgt                                                                                                >  1:46349/2‑135 (MQ=255)
                                    gccgcAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgc                                                                                             >  8:43167/1‑136 (MQ=255)
                                    gccgcAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgc                                                                                             >  7:383841/1‑136 (MQ=255)
                                      cgcAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgcg                                                                                            >  5:69721/1‑135 (MQ=255)
                                      cccAACCGCAGCGACAGCGCTGGGCTGAATCAAACCCGCGCGCAGGCAGGCTGCCCCCCACAGGTTTTTCACGCCCGGGCCGCTCACAGCGAAGTAAACGCAGCGCTAAACCCCGCGACCGAGGCCGCTGTGgcg                                                                                            <  4:257400/133‑1 (MQ=255)
                                       gcAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgcgc                                                                                           >  2:109614/1‑135 (MQ=255)
                                        cAACCTGCGCGACAGCGCTGTGCTGAATGACACACGCCCGCCGGCAGGCTGCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAACGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgcgcg                                                                                          <  5:135541/135‑1 (MQ=255)
                                         aaCCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgcgcgc                                                                                         >  1:130884/1‑135 (MQ=255)
                                         aaCCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgcgcgc                                                                                         >  5:377115/1‑135 (MQ=255)
                                         aaCCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgcgcgc                                                                                         >  4:320767/1‑135 (MQ=255)
                                         aaCCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGgcgcgc                                                                                         >  3:156929/1‑135 (MQ=255)
                                           ccTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCggg                                                                                      >  8:199242/1‑136 (MQ=255)
                                                ccGACAGCGGTGTGCGTAATGAACCCCCCCCCCAGGCAGGCTGCCACCCACTGGTTTTTCATCCCCTGGCCGCTCAAAGCGAAGTAAAAGCAGCCCTAACCCCCGCGACCGAGGCCGCTGTGGCGCGCGGGgtgt                                                                                  <  6:223099/134‑1 (MQ=255)
                                                  gACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGGGttt                                                                                >  2:161364/1‑135 (MQ=255)
                                                     aGCGGTTTGCTAAATCAAACCCTCGCGCAGGATGGCTTCTACCCACAGTTTTTTCATCCCGTGGCCGCTGACAGCGAATAAAACCCCCCGCAAACCCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGgc                                                                             <  6:377115/135‑1 (MQ=255)
                                                          ggggCTTAATGAAACACTCGCGCCGGCAGGCTGCCACCCAAACGTTTTTCATCCCCTGGCCGAAGACCGCTATTAAAACTCACCGCAAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGcc                                                                        <  6:73602/132‑1 (MQ=255)
                                                            tgCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCa                                                                      >  2:266922/1‑135 (MQ=255)
                                                            tgCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCa                                                                      >  1:352334/1‑135 (MQ=255)
                                                            tgCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCa                                                                      >  4:230653/1‑135 (MQ=255)
                                                                   ttaaaCACGCGCGCAGGCTGGCTCCGACCCACTGGTTTTTCATGCCCTGGCCGCTGACAGCGAATAAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGttt                                                               <  2:352334/133‑1 (MQ=255)
                                                                    gaaACACTCGCGCGGGCAGGCTTCGACCCACAGGATTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGCGGTGTTTGGCGCGCCCACCTGtttt                                                              >  4:384399/1‑135 (MQ=255)
                                                                    gaaACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGtttt                                                              >  5:139308/1‑135 (MQ=255)
                                                                    gaaACACTCGCGCAGGCAGGATTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGtttt                                                              >  3:215951/1‑135 (MQ=255)
                                                                         acTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTaaaa                                                                                                                                >  3:284170/1‑64 (MQ=255)
                                                                         acTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTaaaa                                                                                                                                <  4:284170/64‑1 (MQ=255)
                                                                         acTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCgg                                                         >  2:207207/1‑135 (MQ=255)
                                                                         acTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGt                                                        >  8:296032/1‑136 (MQ=255)
                                                                         acTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGGTTTGTCgg                                                         >  1:134592/1‑135 (MQ=255)
                                                                         acTCGCGCAGGCAGGCTTAGACCATCAGGTATTACATGCCCTGGCCTCTGACATCGAACAAAAAGCAGCGCTAAAGCACTCGACCGAGGCCGCTCTGGCGCGCGGGGTGTATGCCGCGCCCACCAGTTTTCTCGGt                                                        >  7:12016/1‑136 (MQ=255)
                                                                                caggcaggCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACggg                                                  >  6:83159/1‑135 (MQ=255)
                                                                                   gccggcTGCGACCCACAGGTTTTTCATGCCCGGCCCGCTGACAGCGAAGTAAAAGCCGCGCTAAAGCCCGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGt                                              <  8:260549/133‑1 (MQ=255)
                                                                                   gcaggcTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGt                                              <  8:375987/136‑1 (MQ=255)
                                                                                   gaagggTTCGACCCACAGGTTTTTCTTTCCCTGGCCGCTGACTGTGATGTAAAAGCTGGGCTAATGCTCGCGACCGAGGCCGGTGTGGCGCGCGGGGTGTTTGGCGCGCGCACCTTTTTTGTCGGTGACGGGATGt                                              <  8:12016/134‑1 (MQ=255)
                                                                                    ccggattCGCCCCACTGTTTTTTCAGGCCGGGGCCGCTTACAGCTAAGAAAACCCCGCGCTAACCCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGt                                              <  4:4186/130‑1 (MQ=255)
                                                                                    acggcTGCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACCGCGAAGTAAAAGCCGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGt                                              <  4:279310/133‑1 (MQ=255)
                                                                                        cTTCGACCCACAGGTATTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCCGGCTGTTTGGCGCGCCCACCGGTTATGTCGGTGACGGGATGttctt                                          >  5:171807/1‑135 (MQ=255)
                                                                                         tccgcccccCAGGTTTTTCCGCCCCGGCCCGCTCACAGCGAAGAAAACCCCCCCCAAACGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGttcttc                                         <  5:336058/130‑1 (MQ=255)
                                                                                          tCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCgg                                       >  7:12701/1‑136 (MQ=255)
                                                                                            gACCCACAGGTTATTCCTGCCCGGGCCGCTGACAGCGAATTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGcc                                     <  7:242133/136‑1 (MQ=255)
                                                                                             ccccACTGGTTTTTCATGCCCTGCCCGCTGACAGCGAAGTAAAAGCAGCGATAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGcc                                     <  2:46349/134‑1 (MQ=255)
                                                                                             aCCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAACGCAGCGCTAACGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGcc                                     <  3:250677/135‑1 (MQ=255)
                                                                                                       ttttCCTGCGGTGGCGGCTGACAGCGTTGTAATAGCCGGGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCTCCCAGCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGccg                           <  1:316778/135‑2 (MQ=255)
                                                                                                          tcctccccGGGCTGGTGACAGCTAATAAAACGCAGCGTTAACCCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTTTTTTGTCGTTGACGGGATGTTCTTCGGCCAGGACCGCCTGGa                        <  7:200588/130‑1 (MQ=255)
                                                                                                           cATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGa                        <  4:325258/134‑1 (MQ=255)
                                                                                                                       gCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGACTTTGTCgaggag           >  3:247621/1‑135 (MQ=255)
                                                                                                                       gCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGACTTTGTCgaggag           >  3:12708/1‑135 (MQ=255)
                                                                                                                              gCGAAGAAAAAGCAGCGCTAAAGCACGCGACCGACCCCGCTGTGCCGCGCCGGGTCATTGGCGCCCCCACCTGTTATGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGACTTTGTCGCGGAGGCGCTa     >  5:181621/1‑134 (MQ=255)
                                                                                                                               ctaaGTAAACCCAGCGATAAAGCCCGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGACTTTGTCGAGGAGGCGCTACGt  <  8:12701/134‑1 (MQ=255)
                                                                                                                                gAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGACTTTGTCGAGGAGGCGCTACGt  <  2:68318/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
TGCTATCGGTGCTGTTCACTGGGCTATGGGTACAGCGCCGCAACCTGAGCGACAGCGCTGTGCTGAATGAAACACTCGCGCAGGCAGGCTTCGACCCACAGGTTTTTCATGCCCTGGCCGCTGACAGCGAAGTAAAAGCAGCGCTAAAGCACGCGACCGAGGCCGCTGTGGCGCGCGGGGTGTTTGGCGCGCCCACCTGTTTTGTCGGTGACGGGATGTTCTTCGGCCAGGACCGCCTGGACTTTGTCGAGGAGGCGCTACGT  >  ACB122/3274861‑3275123

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: