Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,277,0430TG31.3% 92.4 / 15.4 68G278G (GGT→GGG) hbd3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (13/33);  new base G (21/0);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 9.89e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGC  >  ACB122/4276907‑4277173
                                                                                                                                        |                                                                                                                                   
aCCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGCAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGt                                                                                                                                     <  2:20918/137‑1 (MQ=255)
   tGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGgttt                                                                                                                                  >  6:325578/1‑136 (MQ=255)
     aTCGGCCTGTACACCCTGCTGGCAATCATGGAGGCCTGCCACGACGGCTTAAACGACAGAAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGCCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt                                                                                                                                <  8:179088/137‑1 (MQ=255)
     aTCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt                                                                                                                                <  6:67304/137‑1 (MQ=255)
     aTCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt                                                                                                                                <  2:179492/137‑1 (MQ=255)
           cTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTa                                                                                                                          <  2:48046/137‑1 (MQ=255)
            gggACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTAAACGACAGAAAATACCGCCCTGCTCCACTGCTCAAGCAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTAc                                                                                                                         <  1:179531/136‑1 (MQ=255)
             ggACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCGCCGAcg                                                                                                                        >  8:7794/1‑136 (MQ=255)
                      gctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCaa                                                                                                               >  4:10839/1‑137 (MQ=255)
                       ctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGTGGCGCAAGAGCGGTCGCGGTTTAATCACCTACTGAGGAGCaac                                                                                                              >  6:331028/1‑136 (MQ=255)
                       ctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCATCACCTACTGAGGGGCAAg                                                                                                              >  6:331023/1‑137 (MQ=255)
                       ctgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTCACCAACTGAAGACCaac                                                                                                              >  7:207413/1‑136 (MQ=255)
                        tgGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAgc                                                                                                             >  5:245125/1‑137 (MQ=255)
                         ggCAATCCTGGCGGCCGTCCACGACGGCTTCAACGACAGCAAATCCCGCCCGGCTCCCCTGATCAAGCAAATGGCCGCGGCCGGCTGGCTTGGTTGCAATTGCGGGGGCTGTTTCTTGAGCTACTGAGGAGCATgcg                                                                                                            <  4:150249/137‑1 (MQ=255)
                         ggCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAgcg                                                                                                            <  2:45046/137‑1 (MQ=255)
                                aTGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCAACTACTGAGGAGCAAGCGCCATGcc                                                                                                     >  6:330345/1‑137 (MQ=255)
                                aTGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTAACCTACTGAGGAGCAAGCGcccagcc                                                                                                     >  2:174434/1‑132 (MQ=255)
                                 tGGAGGCCTTCCACGAAGGCTTCAACGACACCAAATCCCGCCCTGCTCCACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGccc                                                                                                    <  4:31050/137‑1 (MQ=255)
                                  ggAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCCCCTACTGAGGAGCAAGCGCCATGcccc                                                                                                   >  4:75277/1‑137 (MQ=255)
                                             aCGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTGTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGcc                                                                                        >  4:270234/1‑137 (MQ=255)
                                             aCGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGGTTGTTAACAAAAGAGGAGCAAGCGCCAATCCCCCCGTCaaatcc                                                                                        >  5:40037/1‑133 (MQ=255)
                                              cGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGCGGCGCAAGCGCGGTCGCGGGTTCTTCACCGACAGAGGAGCAAGCGCCCAGCCCCCCGTCAAAGCCg                                                                                       >  6:113920/1‑137 (MQ=255)
                                                   gCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATg                                                                                  <  6:329359/137‑1 (MQ=255)
                                                        aactccAGCAAATACCGCCCTGCTCCACTGCTCACGGAATTGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCtg                                                                             <  1:134424/132‑1 (MQ=255)
                                                        aaCGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTGTTGACCTAATGAGGAGCAAGCCCCAAGCCCCCGGTCAATGCCGCAATTCGCtg                                                                             >  5:305359/1‑137 (MQ=255)
                                                        aaCGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCACCAAAGTCGGGGCAAGCGCCCAGCCCCCCCGCAAAGCCGCCAAGCGCtg                                                                             >  4:186384/1‑137 (MQ=255)
                                                        aaCGACAGCAAATACCGCCCTGCGCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCtg                                                                             <  5:303235/137‑1 (MQ=255)
                                                         aCGACAGCAAATCCCGCCCTGCTCCACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCtgt                                                                            <  1:174434/137‑1 (MQ=255)
                                                         aCGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAAAGCCGCAATGCGCtgt                                                                            >  8:230538/1‑137 (MQ=255)
                                                                cAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTGGTGCCCTACTGAGGAGCAAGCGCCAAGCCCCCGGGCCATGCCGCCATGCGCggtggccacg                                                                     >  1:175883/1‑131 (MQ=255)
                                                                cAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGGTCACCGACGGATGAGCAAGCGCCATGCCCCCGGTCACTGCCGCCATGCGCggtggccacg                                                                     >  6:183066/1‑131 (MQ=255)
                                                                  aaTACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTAACCGACGGAGGAGCCAGCCCCATGCCCCCGGGCAATGCCCCCATGCGCAGTGCCAACttt                                                                   >  6:270645/1‑137 (MQ=255)
                                                                     cccGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTgag                                                                <  3:164061/136‑1 (MQ=255)
                                                                          ccTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGGTGCCCTACTGAGGAGCAAGCGCCATGCCCCCGGTGAATGCCGCAATGCCCTGGGCCAACTTTGAGGAGCg                                                           >  8:140147/1‑137 (MQ=255)
                                                                           cTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCGACTGAGGCGCAAGCGCCCTGCCCCCGGTCAATGGCGCAATTCGCGGTGCCCACGTTGTGGAGCgc                                                          >  6:263406/1‑136 (MQ=255)
                                                                            tGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGAACGACGGAGGAGCAAACCCCCAGCCCCCCGGAAAAGCCGCAAAGCGCTGTGGCCACTTTTAGGGGCgcc                                                         >  7:166638/1‑135 (MQ=255)
                                                                              ccccccTGCACAAGTAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                       <  2:159656/132‑1 (MQ=255)
                                                                              cTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                       <  7:223788/137‑1 (MQ=255)
                                                                              cTCCACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                       <  5:325578/137‑1 (MQ=255)
                                                                              cTCCACTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                       <  2:287463/137‑1 (MQ=255)
                                                                              cTCCACTCCTAAAGTAATTGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCTCTGTGCCAACTTTGAGGAGCGACGt                                                       <  2:287467/137‑1 (MQ=255)
                                                                                   cTGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCACTGTGCCAACTTTGAGGAGCGACGTG‑ACCg                                                 <  3:149005/137‑1 (MQ=255)
                                                                                    tGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCgg                                                <  6:118782/137‑1 (MQ=255)
                                                                                    tGCTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCgg                                                <  4:71441/137‑1 (MQ=255)
                                                                                         aaGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTGTTCGCCGACGGAGTCACCACCGCCCACCCCCCCGTGAAATCCCCAATTCCGTGTGGCAACATTGTGGAGCGACGGG‑GCCGGGCCAt                                           >  6:281400/1‑137 (MQ=255)
                                                                                         aaGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCTACTGGGGAGCCAGCCCCCAGCCCCCCGGCAAAGCCGCCATGCGCTGTGCCCAAATTTTGGAGCGGCGGT‑CCCGGGccca                                           >  8:116405/1‑135 (MQ=255)
                                                                                         aaGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGACCTACTGAGGGGCAAGCGCCATGCCCCCCGGCAATTCCGCAAAGCGCTGTGCCAACATTGAGGAGCGCCGTG‑ACCGGGCCaa                                           >  6:175315/1‑136 (MQ=255)
                                                                                           ggAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTGTTCACCTTCTTATGATCAAGCGCCCTGCCCCCCGTCAATGCCGCAATGGGGTGTTGCAACTTTGTGGTGCGACGTG‑GAGGGGCGATgg                                         >  4:203875/1‑137 (MQ=255)
                                                                                           ggAAATGGTAGCGGCCGGCTGGCTCGGGCGCAAGAGCGGTCGCGGATTCTTAACCTACAGAGGAGCAACAGCCAACGCCCCCGTCAATCCCCCAAAGCGCTAGCGCACCAATGAGGAGCGACCTG‑TCCGCGccgagg                                         >  4:199891/1‑133 (MQ=255)
                                                                                              aaTGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCAt                                      <  6:40037/137‑1 (MQ=255)
                                                                                                   tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCAACTGAGGTGCAAGCGCCCTGGCCCCGGGGAATGACCCAATGAGCGTGGGCAACGTTGATGAGCGCCGAG‑ACCGGGCCCTGGCATTGttt                                 >  2:89484/1‑137 (MQ=255)
                                                                                                   tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGAACGACTGAGGAGCAACCCCCAAGCCCCCCGGCAAAGCCGCCAAGCGCGGGGCGAACATTTAGGGGGCACCTG‑CCCCGGGCCTGGCAGTGttt                                 >  6:221343/1‑137 (MQ=255)
                                                                                                   tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTCACCGACGGAGGGGCAAGCGCCATGCCCCCCGGCAATGCCGCAATGCGCTGTGCCCACGTTGAGGAGCGCCGTG‑CCCGGGCCCTGGGAGTGttt                                 >  6:73682/1‑137 (MQ=255)
                                                                                                   tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCCCCTACGGAGGAGCAAGCGCCCTGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCAAGGCAGTGgtt                                 >  8:243642/1‑137 (MQ=255)
                                                                                                   tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTCACCTAATTAGGAGCAAGCGCCAAGCCCCCCGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGCCGTG‑ACCGGGCCATGGCATTGttt                                 >  7:254948/1‑137 (MQ=255)
                                                                                                   tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTAGTGACCTAAAGAGGAGCAAGCGCCATGCCCCCGGTCAATTCCGCAATGCGCTGTGCCCAAGTTGAGGAGCGCCCTG‑ACCGGGCCAAGGGATTGttt                                 >  7:315530/1‑137 (MQ=255)
                                                                                                        gccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCTCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGa                            <  4:209674/137‑1 (MQ=255)
                                                                                                        gccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGa                            <  1:78/137‑1 (MQ=255)
                                                                                                        gccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGa                            <  6:168144/137‑1 (MQ=255)
                                                                                                             ctggctggGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACGGGGCCTTGGCATTGTTTGCCGAAAAgg                       <  7:274413/137‑1 (MQ=255)
                                                                                                                  tggGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCttt                  <  3:10839/137‑1 (MQ=255)
                                                                                                                  gggggCGCAAGAGCGGGCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCttt                  <  8:271221/136‑1 (MQ=255)
                                                                                                                       cgcAAGAGCGGTCGCGTTTTCTTCACCTACTGAGGAGCAAGAGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGATGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCa             <  7:157042/137‑1 (MQ=255)
                                                                                                                       cgcAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCa             <  6:315583/137‑1 (MQ=255)
                                                                                                                        gcAAGAGCGGTCGCGGGTTATTCACCAACTGAGGAGCAAGCGCCATGCCCCCCGGGCATGCCCCAATGCGCTGTGCCAACATTTAAGAGCGCCGTGCCCCGGCCCATGG‑ATTGTTTGCCGAAAATGGCGTTGGcccc            >  4:241237/1‑135 (MQ=255)
                                                                                                                         aaaGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGCG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAgg           <  8:58935/136‑1 (MQ=255)
                                                                                                                                ggTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCAt    <  5:262468/137‑1 (MQ=255)
                                                                                                                                  tCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATTCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGGCATGGCATTGTTTGGCGAAAAGGGCTTTGGGCAGGTCAGCATgc  >  2:11572/1‑137 (MQ=255)
                                                                                                                                  tCGCGGGTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑CCCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATgc  >  8:173039/1‑137 (MQ=255)
                                                                                                                                  gcgcgGTTGCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCGGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATgc  <  7:243642/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
ACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTG‑ACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGC  >  ACB122/4276907‑4277173

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