Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,874,7890TC52.6% 77.7 / inf 545L21L (TTA→TTGcreDprotein CreD
Reads supporting (aligned to +/- strand):  ref base T (112/146);  new base C (125/161);  total (238/307)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

Only 100 of 549 total aligned reads displayed.

TGCACCGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCG  >  DH1/3874640‑3874937
                                                                                                                                                     |                                                                                                                                                    
 gCACCGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTa                                                                                                                                                     >  4:1397019/1‑150 (MQ=255)
  cACCGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTaac                                                                                                                                                   >  8:1229927/1‑151 (MQ=255)
  cACCGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCaac                                                                                                                                                   >  6:1046343/1‑151 (MQ=255)
   aCCGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGACTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCCTCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCaac                                                                                                                                                   >  5:1667362/1‑150 (MQ=255)
    ccGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTaacaa                                                                                                                                                 <  7:1905218/151‑1 (MQ=255)
    ccGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTaacaa                                                                                                                                                 <  4:528995/151‑1 (MQ=255)
     cGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAAACACCTGCCGAATCAGCATTATCGGAATTaacaac                                                                                                                                                >  1:153349/1‑151 (MQ=255)
      gTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTaacaac                                                                                                                                                >  8:538918/1‑150 (MQ=255)
      gTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCaacaac                                                                                                                                                >  7:209132/1‑150 (MQ=255)
       tataAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTGCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAacaacac                                                                                                                                              >  2:1078839/1‑150 (MQ=255)
             gCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATAc                                                                                                                                        <  6:691105/151‑1 (MQ=255)
                cGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGc                                                                                                                                     >  7:303269/1‑151 (MQ=255)
                 ggTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGcc                                                                                                                                    >  7:1146421/1‑150 (MQ=255)
                    cacaGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCa                                                                                                                                 <  8:1272155/151‑1 (MQ=255)
                     acaGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCaa                                                                                                                                <  5:401180/151‑1 (MQ=255)
                     acaGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCa                                                                                                                                 >  4:147946/1‑150 (MQ=255)
                       aGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTACACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTTTCGGAATCAACAACAGCAATACTGCACCAAAc                                                                                                                              >  2:1436069/1‑151 (MQ=255)
                        ggAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACa                                                                                                                             <  5:305211/151‑1 (MQ=255)
                         gAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAg                                                                                                                            <  2:1155646/151‑1 (MQ=255)
                          aaTAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACTACAGCAATACTGCACCAAACAgg                                                                                                                           >  1:616102/1‑151 (MQ=255)
                             aGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTa                                                                                                                        <  2:683659/151‑1 (MQ=255)
                             aGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTa                                                                                                                        <  4:950933/151‑1 (MQ=255)
                              gCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGATCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAg                                                                                                                       <  2:348548/151‑1 (MQ=255)
                                     gCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCAtttt                                                                                                                <  4:810906/151‑1 (MQ=255)
                                        gCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCa                                                                                                             <  3:683698/151‑1 (MQ=255)
                                        gCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCa                                                                                                             <  5:1789920/151‑1 (MQ=255)
                                          ccAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGa                                                                                                           <  5:1149139/151‑1 (MQ=255)
                                            aaCGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAAc                                                                                                         >  3:1331056/1‑151 (MQ=255)
                                               gAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAgg                                                                                                      <  7:1093309/151‑1 (MQ=255)
                                                    tttGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAgggggg                                                                                                  <  5:615042/150‑1 (MQ=255)
                                                       gCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGAttt                                                                                              <  3:204335/151‑1 (MQ=255)
                                                       gCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGAtt                                                                                               <  8:1084308/150‑1 (MQ=255)
                                                         ggCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTc                                                                                             >  4:1672742/1‑150 (MQ=255)
                                                           cccGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAAc                                                                                          >  4:246207/1‑151 (MQ=255)
                                                            ccGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACa                                                                                         <  3:11985/151‑1 (MQ=255)
                                                             cGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACAt                                                                                        <  6:818116/151‑1 (MQ=255)
                                                                 ggTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGtc                                                                                    <  5:1325396/151‑1 (MQ=255)
                                                                 ggTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGtc                                                                                    >  5:228507/1‑151 (MQ=255)
                                                                    aCTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGtcttc                                                                                 >  7:1747351/1‑151 (MQ=255)
                                                                    aCTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGtctt                                                                                  <  8:190994/150‑1 (MQ=255)
                                                                       ttGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTtctc                                                                               <  4:1264210/150‑1 (MQ=255)
                                                                         ggCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCtt                                                                            <  8:44065/151‑1 (MQ=255)
                                                                          gCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCttt                                                                           >  1:855018/1‑151 (MQ=255)
                                                                           cGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGCCAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCGAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTg                                                                          >  2:1682573/1‑151 (MQ=255)
                                                                             aaTCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGc                                                                         >  1:1666490/1‑150 (MQ=255)
                                                                                cGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATg                                                                     >  6:460928/1‑151 (MQ=255)
                                                                                cGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAAt                                                                      <  6:1139775/150‑1 (MQ=255)
                                                                                       tCCACATCGCTACGGTAATCAGCACGGTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGgcagc                                                                <  5:1560330/149‑1 (MQ=255)
                                                                                        ccaaaacggtACCGTAAACAACACGTTCCAAAATCACCTTCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGcagcag                                                              <  8:1492977/141‑1 (MQ=255)
                                                                                        ccACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGcagcag                                                              <  8:643316/150‑1 (MQ=255)
                                                                                         cacaTCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAgg                                                             >  2:601601/1‑150 (MQ=255)
                                                                                         cacaTCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATAAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCATCAgg                                                             >  2:884859/1‑150 (MQ=255)
                                                                                          acaTCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGAt                                                           <  6:316321/151‑1 (MQ=255)
                                                                                           caTCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGAt                                                           >  3:1104773/1‑150 (MQ=255)
                                                                                              cGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGa                                                       <  3:27868/151‑1 (MQ=255)
                                                                                              cGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGAGTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACg                                                        <  3:855323/150‑1 (MQ=255)
                                                                                               gCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCATCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGaa                                                      <  8:1659634/151‑1 (MQ=255)
                                                                                               gCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGaa                                                      <  4:1009648/151‑1 (MQ=255)
                                                                                                  aCGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGa                                                    <  1:1384935/150‑1 (MQ=255)
                                                                                                  aCGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGAc                                                   >  1:499501/1‑151 (MQ=255)
                                                                                                    ggTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTa                                                 <  3:480469/151‑1 (MQ=255)
                                                                                                     gTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACACGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTAt                                                <  4:75867/151‑1 (MQ=255)
                                                                                                     gTAATCACCACCGACGACAATCCCCTGCCGAATCAGCATTATCGGAATCAACAACATCAAAACTGCACCAAACACGCTAGTCATTTTCCAGAACATTTGCGATTTCAACATAGTCTTATCCTTTGCAATGGCAGCAGGATACGAAGACTAt                                                <  6:1615729/151‑1 (MQ=255)
                                                                                                      tAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTAtg                                               >  4:646325/1‑151 (MQ=255)
                                                                                                       aaTCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTAtgtg                                             <  8:307297/151‑1 (MQ=255)
                                                                                                       aaTCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTAtgt                                              <  1:1249550/151‑1 (MQ=255)
                                                                                                       aaTCACCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTACTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTAtgt                                              <  8:947489/151‑1 (MQ=255)
                                                                                                        aTCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTAtgtg                                             <  8:710134/151‑1 (MQ=255)
                                                                                                           aGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTggg                                           >  3:1546390/1‑150 (MQ=255)
                                                                                                            gCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGga                                          >  8:1006900/1‑150 (MQ=255)
                                                                                                              aCGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGga                                          >  1:1281648/1‑148 (MQ=255)
                                                                                                               cGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGgaagaa                                      <  7:461568/151‑1 (MQ=255)
                                                                                                               cGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGgaaga                                       >  4:400513/1‑150 (MQ=255)
                                                                                                                gTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAAt                                     >  8:869786/1‑151 (MQ=255)
                                                                                                                  tCGACAATCACCTGCCGAATCAGCATTATCGGTATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAAtt                                    <  6:1721480/150‑1 (MQ=255)
                                                                                                                   cGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTg                                  <  7:553536/151‑1 (MQ=255)
                                                                                                                      cAATCACCTGCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTNa                                 >  1:1793851/1‑149 (MQ=255)
                                                                                                                      cAATCACCTGCCGAATCAGCATTATAGGAATTAACAACAGCAATACTGCACCAAACAGGATAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAg                               <  5:1093023/151‑1 (MQ=255)
                                                                                                                          cACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTa                            <  5:86903/150‑1 (MQ=255)
                                                                                                                               gCCGAATCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTAtgtg                        <  7:1610101/149‑1 (MQ=255)
                                                                                                                               gCCGAATCAGCATTATAGGAATCAACCACAGCAATACTGCACCCAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGa                       <  6:1254971/150‑1 (MQ=255)
                                                                                                                                ccGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAg                     <  7:1348815/151‑1 (MQ=255)
                                                                                                                                ccGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAg                     <  2:183299/151‑1 (MQ=255)
                                                                                                                                   aaTCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGa                  <  2:278799/151‑1 (MQ=255)
                                                                                                                                   aaTCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGa                  <  6:1340775/151‑1 (MQ=255)
                                                                                                                                   aaTCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGa                  <  2:126308/151‑1 (MQ=255)
                                                                                                                                    aTCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGAc                 >  3:1794612/1‑151 (MQ=255)
                                                                                                                                     tCAGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACg                >  4:175913/1‑151 (MQ=255)
                                                                                                                                      cAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACgg               >  7:748685/1‑151 (MQ=255)
                                                                                                                                       aGCATTATCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGt              <  1:1120119/151‑1 (MQ=255)
                                                                                                                                         cATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTg             >  4:1103297/1‑150 (MQ=255)
                                                                                                                                            tATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAg          <  3:1862166/150‑1 (MQ=255)
                                                                                                                                             aTCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTc        >  6:1812389/1‑151 (MQ=255)
                                                                                                                                             aTCGGAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCACCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTc        <  7:139238/151‑1 (MQ=255)
                                                                                                                                                ggAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGCCGGTGAAGTCGa      >  3:65965/1‑150 (MQ=255)
                                                                                                                                                ggAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGaa     <  1:1000734/151‑1 (MQ=255)
                                                                                                                                                ggAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGa      >  1:291767/1‑150 (MQ=255)
                                                                                                                                                ggAATCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGaa     <  1:1054261/151‑1 (MQ=255)
                                                                                                                                                  aaTCAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGc   <  2:1565232/151‑1 (MQ=255)
                                                                                                                                                   aTTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATCTGGGAAGAATTTGACGCTATGTGAAGTGACGGTGAAGTCGAAGCg  >  1:39688/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCACCGTATAAAGCTCGGTCACAGGAATAGCGATGAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCATTATCGGAATTAACAACAGCAATACTGCACCAAACAGGCTAGTCATTTTCCAGAACAGGGGGGATTTCAACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCG  >  DH1/3874640‑3874937

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