Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,906,8320CA31.3% 137.2 / 11.9 98intergenic (+3/‑32)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (33/24);  new base A (26/0);  total (73/25)
Fisher's exact test for biased strand distribution p-value = 1.83e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCT  >  DH1/3906689‑3906981
                                                                                                                                                 |                                                                                                                                                     
gcgcgccacgccccccccTCCCACCTAAAACCACTTTTAAATTGCTGCGAGTTTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGa                                                                                                                                                <  1:1126881/141‑1 (MQ=14)
cccgccaacgcccgccccTCGCACCAACAACTCCTTTTTAATTCCTGTCCGTGGGCGAAGGTGGCGTAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGCCCCCCCCTCGCACCAAACGAGGCGa                                                                                                                                                <  5:1183044/141‑1 (MQ=2)
 gcgcccattCCCCCCCCGCGCTCCAAAAAACACTTTGAAATTGATCCCCGTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTGCAAGTCCCCCCCCTCGCACCAAACGAGGCGat                                                                                                                                               <  3:1564092/143‑1 (MQ=17)
    gcAAGGCCCCCCCCTCGCACCAAAAACCACGTTGATATTGGTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATc                                                                                                                                            <  7:51108/150‑1 (MQ=255)
       acgccccccccTCCCACCAACAACCACTTTGAAATACCTGACTCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaaa                                                                                                                                         <  8:1459421/148‑1 (MQ=18)
         gccgcccccTCCCACCAACCAACACTTTTATATAGCTGCCTGTCTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaaaaa                                                                                                                                       <  3:257814/147‑1 (MQ=14)
               cccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGCACGGGAGATCATATAAAGttt                                                                                                                                 >  6:142606/1‑149 (MQ=18)
               cccTCGCACCAAAAACCACGTTGAAATTGCTCGCTCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTaa                                                                                                                                 <  3:531226/151‑1 (MQ=255)
                ccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAATAACGGTACATCAATAAAAGtttc                                                                                                                                >  1:1045657/1‑148 (MQ=18)
                  tCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAACGAGGCGAGATCAAATAATGttacat                                                                                                                              >  1:759395/1‑146 (MQ=21)
                  gcccACCAACCACCACGTTTATATAGCTCCCCGTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAACGAGGCGATATCAAAAAAAGTAAGAt                                                                                                                              <  5:906139/148‑1 (MQ=18)
                   cgcacctaccaactcGTTGATAGAGCTCGCAGTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTGCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATg                                                                                                                             <  5:537852/137‑1 (MQ=18)
                   cGCACCAAAAACCTCGTTGAAATTCCTCTCTGTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATg                                                                                                                             <  4:463157/151‑1 (MQ=18)
                   cGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTGAGGCTAGTTCAAATAAAGTTTGATg                                                                                                                             >  1:1117189/1‑151 (MQ=18)
                   cGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCCCCAAAGAAGGGGAGATAAAAAAAAGTatggt                                                                                                                              >  8:579592/1‑146 (MQ=18)
                    gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTACCCCCCCTCGCACCAAAAGACGGTAAATCAAAAAAAGTAAgcaa                                                                                                                             >  4:819344/1‑147 (MQ=18)
                      acaaaacaCCACTTTGATATTCCGCGCAGTGGTCGAATGTGGAGGACTTGGTAGACCCGATAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACt                                                                                                                          <  7:585702/148‑1 (MQ=18)
                      aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATAAAAAAAAGTAAGATGACt                                                                                                                          >  2:1011722/1‑151 (MQ=255)
                      aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGCCGAAGACCAAAAATGTAACAtcc                                                                                                                            >  7:1647517/1‑147 (MQ=18)
                      aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAACTCCCCCCCCTAGCCCCAAAAATCGTTATATCAAAAAAAGtttcttacat                                                                                                                          >  8:1634033/1‑142 (MQ=18)
                       ccAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAACTCGGGATATCAAAAAAAGTTAGATGAAtg                                                                                                                         >  3:294041/1‑151 (MQ=18)
                        cAAAAACCACGTTGATATTGCTCGCACTGTTGTTAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTGTTACGGTCGTGGGGGGTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACtgt                                                                                                                        <  1:504150/151‑1 (MQ=21)
                         aaacacgtcGTTTATATTGCTCGGTCTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGGCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACtgtg                                                                                                                       <  1:1164885/143‑1 (MQ=18)
                         aaaaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACtgtg                                                                                                                       >  1:954061/1‑151 (MQ=255)
                          aaaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAAAGGTTAAATCAAATAAAGTATAATGCCtgtgt                                                                                                                      >  7:1386430/1‑150 (MQ=18)
                          aaaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAACGAGGAGATATCAATTAATGTTACATGAGtgtgt                                                                                                                      >  1:1317786/1‑150 (MQ=18)
                           aaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAATAAGGGGATATCAATTAATGTTTGTTGACtgtttt                                                                                                                     >  1:1097415/1‑148 (MQ=17)
                           aaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGGATATCAAATATATTTAGATGCAtgtggt                                                                                                                     >  8:598291/1‑149 (MQ=18)
                           aaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTACCCCCCCACGCACCAAAACAAGGGAAAACAAAAAAAGGTAGTAGACtgtggt                                                                                                                     >  2:1368977/1‑149 (MQ=17)
                                aCGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGt                                                                                                                <  1:1010336/151‑1 (MQ=255)
                                 cGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTg                                                                                                               <  2:954527/151‑1 (MQ=255)
                                  gTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAGGAGGTTAGATCCATTAAAGTTTGATGCCtgttggttggtcg                                                                                                              >  3:749096/1‑141 (MQ=17)
                                  gTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGACGTTAGATCAAAAAATGTTAGATgcgtgtttgttgctgg                                                                                                              >  3:741187/1‑136 (MQ=17)
                                  gTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGTCGAGATCATAAAAAGTTTGATGCCtgttgtttgctgg                                                                                                              >  3:803592/1‑141 (MQ=17)
                                  gTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGTTATTTCAAATAAAGTTTGTTGACtgatttatcatgc                                                                                                              >  1:1473948/1‑140 (MQ=17)
                                    tGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGACGGGATATCAAATATTGTTAGATGCGtgtgttttggtggcg                                                                                                            >  3:369037/1‑140 (MQ=17)
                                      taattagtcattgtGGGCGAAGGTGTCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTATTACTGACGTGGGGGTGCAAGTCCCCCCCCTCGCACGAAACGAGTCGATATCTATAAAAGTAAGATGAGTGTGCGAAGGTGGCgg                                                                                                           <  6:1189678/137‑1 (MQ=18)
                                       tatTGCTCGCACTGGGCGAAGGTAGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAACAGGGGAATTCAATAAAAGTTTCATGACtgtgttatgttgccggat                                                                                                         >  3:878639/1‑137 (MQ=14)
                                       aatttgttcctctctgCGAAGGGTTCGTAATTGGTAGAGGAGCTAGCTTTAGGTGTTAGTGTTATTACGGACGTGGTGGTTCAAGTCCCCCCCCTCGCACTAAACGATCCGATATCAATAAAAGTAATATGTCTGTGCGAAGGTGGCGGaa                                                                                                         <  4:403560/136‑1 (MQ=12)
                                        atTTCTAGCAGTAGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTGCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGaa                                                                                                         <  3:928117/150‑1 (MQ=255)
                                        atTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAAAGAGGGGAGATCAAAAAATGTTAGTTGCCTGTTGGTTGGTGGCGGatt                                                                                                        >  7:111141/1‑149 (MQ=17)
                                          tGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTCGCACCAAAAATGGGGAGATCAAAAAAAGTTATATGCCTGTTGTTTGGTGGCGGAAttt                                                                                                      >  1:550118/1‑150 (MQ=14)
                                          tGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCCAAATACGGGATATCAAAAAAAGTTAGATGAGTGTGGGTTGGTGGCGGAAtgg                                                                                                      >  3:1634210/1‑149 (MQ=17)
                                          tGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCTACATCAAAAAATGTAAGATGAGTGTTGGTTGCTGGCGGAATTg                                                                                                      >  5:498256/1‑151 (MQ=17)
                                              ggccgtgcgCGACGGGGGCGGAATTGGTAGACGCGCTAGCGGCCGGTGTTAGGGTTCGTCGGGACGTGGGGGTTCAAGTCCCCCCCCTCCCCCCACCCGCGGCGATTTCAAAAAAAGTAAGCTGACTGTGCGAAGGTGGCGGAATTGGTa                                                                                                   <  3:1375130/142‑1 (MQ=18)
                                               gCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTACCCCCCCTCGCACCAAACGACGGTAAATCAAAAAAAGTTAGATGAATGTTATTAGGTGGCGGAAttaattt                                                                                                  >  7:436661/1‑145 (MQ=17)
                                                cACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGCGATATCAAAAAAAGTTAGATGCCTGTGGTAAGGTGGCGGCATTGGTATAc                                                                                                >  3:487343/1‑151 (MQ=18)
                                                    gggCGAAGGGGGCGGACCCGGCAGCCGCGCTAGCTTCAGTTTTTAGTGTGCTTACGGAGGGTGCGGTTCAAGGCGCAGGGCTCGGCAGTTCAGACGCGATATCAAAAAATGTAATCTGACTGTGCGAAGGTGGCGGAATTGGTAGAcgc                                                                                              <  7:1494547/149‑1 (MQ=255)
                                                       cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAAAGGTTAAATAAAAAAAAAGAAGAAGCCTGTGCCAACGAGGCCCCATTAGGAGAcgcggttg                                                                                         >  5:931883/1‑147 (MQ=2)
                                                            gTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAAAGTAAGATGCCTGTTGTTTGGTGGCGGATTTGGTataatagctatcatc                                                                                     >  1:559801/1‑136 (MQ=17)
                                                            gTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGGTCAATTACCCCCCCACGCCACCAAATAGGGGAAATCAAAAAAAGTAAGATGAATTTGAGAAGGTTGGGGAATTGTTataaacgaaaaaatc                                                                                     >  8:1443377/1‑136 (MQ=12)
                                                                         tAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTACCCCCCCTCGCACCAAAAAAGGTTAAATCAAAAAATGTTTGATGATTGTTTTTTGGTGCCGGAATTAGTATAATTGATATCTGCAGGTTTTAgtgaa                                                                        >  2:779348/1‑148 (MQ=2)
                                                                             cgcgCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACTAGGCGATATCAAAAAAAGTAAGATGCATGTGAGAAGGTGGCGGAATTGGTATACGCGCTAGCATCAGGTGTTAGTGTACtttc                                                                   >  5:247684/1‑149 (MQ=25)
                                                                               cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGGTGAAATCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAATTACAAGCCTGTTCGAAGTTGGCGGAATTGGTAGACGCGCTACGGTGAGGTGTTACTGTGCTTAAgg                                                                 >  5:927654/1‑151 (MQ=18)
                                                                                 cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGt                                                                                                                <  1:149179/102‑1 (MQ=255)
                                                                                  tAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAATAGGGGATTTCAAAAAAAGTTTGATGCCTGTGGGTAGCTGGCGGAATTGGTTTCCGGGCTAGCATCAGGTTTTACTGTCGTTACGGACg                                                              >  3:1340957/1‑151 (MQ=18)
                                                                                         aGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGCGATATCAAAAAATGTTACTTGAATGTGATTTCGTGGCGGATTTTGTATAAGAGCTATCATCAGGTGTTACTCACCATACGGACGGggcggg                                                       >  1:1557/1‑150 (MQ=12)
                                                                                         aGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGTCGAGTTCATAAAAAGTTAGATGAGTGTTGGTTCGTGGCGGAATTTGATGAAGTTGTAGCTTCAGGTTTTATTGTCGATACGGCAGGGGGGGt                                                       >  5:805386/1‑151 (MQ=9)
                                                                                         aGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGGATATCAATAAAAGTTAGATGCCTGTGGGATGGTGGCGGAATTGGTATACGCTGAATCATCAGGTTTTAGTGTACTTACGGCAGTGGGGGt                                                       >  5:455288/1‑151 (MQ=21)
                                                                                           gtgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCTAAATCAAAAAAAGTAACATGACTGTGCGTACCTGGCTGAATTGGTTGACGCGGGAGCTTCAGGGGGGAGTCTCATTTCGCCCGAGGGGGttt                                                     >  6:311169/1‑150 (MQ=16)
                                                                                            tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAATGTTAGTTGCCTGTTGTATGGTGGCGGAATTGGTATAAGGTATATCTTCAGGATTTAGTGTCATTACGGACGGGGGGGttg                                                     >  3:37827/1‑149 (MQ=16)
                                                                                             gtTAGTGTTCTTACGGACGTGGGGGTTGAAGTGCCCCCCCGCGCCCCCAAAGAGGGGAGATGAAAAAAAGTTAGATGAAGGGTCGAAGGGGGGGGGATTGGGAGACCCCCTAGCGTGAGGGGTTAGTGGGCGTTCGGGCGTtgggggtaac                                                   >  5:360893/1‑146 (MQ=255)
                                                                                             gtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAACAGGGGATATCAAAAATAGTTTGTTGAGTGTGAGTTGGTGGCGGAATTGATATAAGGGGTATCGGCAGGTGTTAGTGTAATAAGGGCCGGGGGGGttggt                                                   >  5:918941/1‑148 (MQ=12)
                                                                                                  tgtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCACCAAAAGAGGCGATATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGGTAGACGGGGTAGCGTCAGGTGGTAGTGTAATTAAGGACGTggggggtggactaac                                              >  8:946467/1‑141 (MQ=25)
                                                                                                     tCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAAAGTTAGATGACTGTGCGTAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTGAAGTcacccc                                           >  5:232448/1‑151 (MQ=25)
                                                                                                           ggACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAACGAGGTTATATCAAAAAAAGTATCATGCGTGAGCTAAGGTGGCGGAATTGGGATACGGGCAAGCATCAGGTTTTACTGACCTTAAGGCCGGGGCGGGTGAAggaaacccccaagc                                     >  7:1551656/1‑138 (MQ=255)
                                                                                                            gACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGTGGATAGATCAAAAAATGTAAGATGACTGATATAACGTGGCGGAATTGATTTACTCGGTATCTTCATGAGTTAGTCTAATTCAGGACGTCGGGGTGCAAggaacaaccaatgcc                                    >  3:1521894/1‑137 (MQ=255)
                                                                                                            gACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGGCCCCCCCCTCGcc                                    >  3:538419/1‑150 (MQ=255)
                                                                                                            gACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGAGATCAAAAAAAGTATGATGCCTGTGGGAAGGTGGCGGATTTGGTTTCCTCGCTAGCTTCAGGTGTTAGTGTCGTTTCGGCCGTGGGGGTTGTAGGCCCCCGCCACGcc                                    >  7:983538/1‑150 (MQ=21)
                                                                                                                       tCAAGTCCCCCCCCTCGCACCAAAGGAGGCGATATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGGTATACGCGGTAGCTTCAGGTGTTAGTGTACTTACGGCCGTGGGGGTGGTAGGCCCCGCCCACGCCCCaagtttttt                         >  7:710298/1‑144 (MQ=21)
                                                                                                                            tCCCCCCCCTCGCACCAAAGGAGGCGATATCAAAAAAAGTAACATGAATGTGAGAAGGTGGCGGATTAGGTTGACGAGCTATCTTCAGGATTTAGTGAAATTAAGGACGTGGGGATTCAAGAACACACCATGGCACCAAGTATAttttcgt                    >  2:1200390/1‑147 (MQ=255)
                                                                                                                            tCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCCAGTCCCCCCCCGCGCACCAATTATCTTTAcgt                    >  5:468774/1‑149 (MQ=255)
                                                                                                                            tCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGGCCCCCCCCCCCCCCCAATTATCTTTACtt                    >  1:1372927/1‑151 (MQ=25)
                                                                                                                              cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGCCGTGGGGGTTCAAGTCCCCCCCCGAGCCCCAAATATATTTTCttgc                  >  6:511029/1‑149 (MQ=21)
                                                                                                                                 ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCttt               >  4:639534/1‑151 (MQ=255)
                                                                                                                                   ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCt             >  5:685155/1‑151 (MQ=255)
                                                                                                                                    cTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCGCCCCCCAAATATCTTTTCTTGCTTTGtt            >  7:1472720/1‑151 (MQ=21)
                                                                                                                                     tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCACGGCCCCCCCCGCGCCCGCATTATCGTTGCTTCCGTTCTTg           >  8:977089/1‑151 (MQ=21)
                                                                                                                                      cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTAGCACCAATTATCTTTTCTTCCTTTCtttt          >  8:70251/1‑149 (MQ=255)
                                                                                                                                      cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAATAACTTTACTTCCTTTCTTGt          >  1:1555829/1‑151 (MQ=255)
                                                                                                                                      cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGTCCCCCCCCGCGCCCCAATTATCTTTACTTCCTTTCTTGt          >  2:1379097/1‑151 (MQ=255)
                                                                                                                                       ccccACCAAACGAGTCCATATCTTATAATTTATGATGACTGTGCGAAGGTGTCGTAATTTGTAGACGCGCTTGCTTCAGTTGTTATTGTTCTTACGTACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTg           <  8:941488/148‑1 (MQ=17)
                                                                                                                                        gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCTTGtt         >  5:1016014/1‑151 (MQ=255)
                                                                                                                                        gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGtt         >  5:1099639/1‑151 (MQ=255)
                                                                                                                                        gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGtt         >  3:1004656/1‑151 (MQ=255)
                                                                                                                                        gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCCGCACCAATTATCTTTACTTCCGTTCTTGtt         >  3:275068/1‑151 (MQ=255)
                                                                                                                                          cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGttt        >  2:376749/1‑151 (MQ=255)
                                                                                                                                          cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCCAGTCCCCCCCCGCGCACCAATTATGTTTACTTCCGTTCTTGttt        >  5:1190096/1‑151 (MQ=255)
                                                                                                                                           cccAAACGATGCTACATCAAAATACGCACTATGCCTCTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttc       <  8:657487/150‑1 (MQ=17)
                                                                                                                                           cccAAACGAGGCGATATCAAAATATTTAATATTAATGTGCGAAGGTGGCTGAATTCCAATACGCGCTAGCTTCTTGTCTTAGACGACTTTCTTACGTGGTGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttc       <  8:1101530/150‑1 (MQ=16)
                                                                                                                                           acgaaaCTATGCTATCTCAAAAAAAGTAATATGACTTTGCGCAGGTTGCGTATTTGTTAGACTCGCTATTGTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttc       <  6:360055/148‑1 (MQ=18)
                                                                                                                                             cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTACCCCCCCGAGGCCCAATTATCTTTACTTACGTTGTTGTttctt     >  8:1164916/1‑151 (MQ=21)
                                                                                                                                              aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttcttc    >  1:663602/1‑151 (MQ=255)
                                                                                                                                              aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttcttc    >  3:911207/1‑151 (MQ=255)
                                                                                                                                              aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttcttc    >  4:1221391/1‑151 (MQ=255)
                                                                                                                                                aCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCTCGCCCCAATTATCTTTACTTGCTTTGTTGTTTCTTCCt  >  6:1188591/1‑151 (MQ=255)
                                                                                                                                                aCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCt  >  5:1109118/1‑151 (MQ=255)
                                                                                                                                                aCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACCCACCAATTATCTTTAATTCCTTTCTTGGttattgc   >  8:583822/1‑148 (MQ=21)
                                                                                                                                                 |                                                                                                                                                     
GGGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCT  >  DH1/3906689‑3906981

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