Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 3,894,823 T→A intergenic (‑23/+104) deoA ← / ← deoC thymidine phosphorylase/2‑deoxyribose‑5‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,894,8230TA98.4% 1914.5 / 8.9 508intergenic (‑23/+104)deoA/deoCthymidine phosphorylase/2‑deoxyribose‑5‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base A (259/241);  total (262/246)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 511 total aligned reads displayed.

GAGGGCGGCAATCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGTATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGACCCAGCGCTT  >  DH1/3894675‑3894970
                                                                                                                                                    |                                                                                                                                                   
gAGGGCGGCAATCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCTAGAAACAAGGTAACCTCACGAGGAAAAGAGAAt                                                                                                                                                   >  1:1370831/1‑151 (MQ=255)
gAGGGCGGCAATCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAAt                                                                                                                                                   >  6:765293/1‑151 (MQ=255)
    gcggcAATCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAgg                                                                                                                                               <  8:1698102/151‑1 (MQ=255)
      ggcAATCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAggcg                                                                                                                                             <  6:1158683/151‑1 (MQ=255)
       gcAATCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAggcgg                                                                                                                                            <  1:686086/151‑1 (MQ=255)
         aaTCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCt                                                                                                                                          >  7:639034/1‑151 (MQ=255)
              gCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCttt                                                                                                                                        <  8:1331735/148‑1 (MQ=255)
              gCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGc                                                                                                                                     >  7:1504201/1‑151 (MQ=255)
                  ttCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCgg                                                                                                                                  >  4:1459750/1‑150 (MQ=255)
                  ttCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCgg                                                                                                                                  >  1:333609/1‑150 (MQ=255)
                    cGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCAAGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGATGAAAAGAGAATGAGGAGTCTTTTGCCGGAt                                                                                                                                <  6:1350125/150‑1 (MQ=255)
                     ggAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGc                                                                                                                              >  1:1307020/1‑151 (MQ=255)
                     ggAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGc                                                                                                                              <  2:79758/151‑1 (MQ=255)
                       agaTAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGt                                                                                                                            >  2:1730727/1‑151 (MQ=255)
                         aTAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCt                                                                                                                          >  7:110806/1‑151 (MQ=255)
                            gtgtTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACg                                                                                                                       >  7:235842/1‑151 (MQ=255)
                             tgttgtCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGc                                                                                                                      >  5:1283305/1‑151 (MQ=255)
                             tgttgtCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGc                                                                                                                      >  2:650349/1‑151 (MQ=255)
                                tgtCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCtt                                                                                                                   >  7:391484/1‑151 (MQ=255)
                                     cgAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCg                                                                                                              >  5:1571098/1‑151 (MQ=255)
                                     cgAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCACCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGAGAGAAACAAGGTAACCTCCCCAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATccc                                                                                                              >  1:1585695/1‑150 (MQ=255)
                                      gAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCgg                                                                                                             >  5:285337/1‑151 (MQ=255)
                                        aTACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGcc                                                                                                           >  2:1190884/1‑151 (MQ=255)
                                        aTACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGcc                                                                                                           >  5:1699270/1‑151 (MQ=255)
                                          aCCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCt                                                                                                          <  8:1424064/150‑1 (MQ=255)
                                            cGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTATTGCCGGATGCGTCTACGCCTTATCCGGCCTACg                                                                                                       >  4:187283/1‑151 (MQ=255)
                                             gTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACgg                                                                                                      >  5:813861/1‑151 (MQ=255)
                                                gATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGt                                                                                                    >  4:315078/1‑150 (MQ=255)
                                                   aaaGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGGAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGt                                                                                                    >  6:38364/1‑147 (MQ=255)
                                                   aaaGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGAt                                                                                                <  4:1594059/151‑1 (MQ=255)
                                                    aaGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATa                                                                                               <  3:413601/151‑1 (MQ=255)
                                                     aGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATAc                                                                                              <  2:381403/151‑1 (MQ=255)
                                                     aGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATAc                                                                                              >  3:680107/1‑151 (MQ=255)
                                                      gAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACg                                                                                             <  7:1298504/151‑1 (MQ=255)
                                                             aTTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGt                                                                                       >  5:1190658/1‑150 (MQ=255)
                                                              tttCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGt                                                                                       >  4:1727824/1‑149 (MQ=255)
                                                               ttcttcATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAgg                                                                                    >  7:778995/1‑151 (MQ=255)
                                                                tcttcATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGc                                                                                   >  3:742958/1‑151 (MQ=255)
                                                                           cAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAagc                                                                          <  5:1077924/149‑1 (MQ=255)
                                                                              ggcATGAATATCTAGTTTTATAAGAATAATTTCTTTTGCGATAAACAATGTAACTTGCTCAGTACAAGAGAATGAGGCGGCTTTTGCCGGATGCGTATACGCCTTATAAGGCCTACTGGTAGATACGTCTCTTAGGCCTTATAagcgaagc                                                                     <  6:1832718/150‑1 (MQ=255)
                                                                              cgcATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGGCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAagcgcagc                                                                     <  5:626675/151‑1 (MQ=255)
                                                                               gcATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAagcgcagcg                                                                    >  6:1024909/1‑151 (MQ=255)
                                                                                   gACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATc                                                                >  1:1472700/1‑151 (MQ=255)
                                                                                     ccATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAg                                                              <  3:1428121/151‑1 (MQ=255)
                                                                                      cATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAgg                                                             <  3:1027902/151‑1 (MQ=255)
                                                                                       aTCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                            <  3:161627/151‑1 (MQ=255)
                                                                                         cACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGt                                                          >  4:1374717/1‑151 (MQ=255)
                                                                                         cACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGt                                                          >  4:1536432/1‑151 (MQ=255)
                                                                                         cACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGt                                                          >  7:863969/1‑151 (MQ=255)
                                                                                          aCGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGATTGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTa                                                         <  5:1788442/151‑1 (MQ=255)
                                                                                           cGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCTAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTaa                                                        <  2:1364658/151‑1 (MQ=255)
                                                                                               ttttACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGc                                                     >  3:1222696/1‑150 (MQ=255)
                                                                                               ttttACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGc                                                     >  2:1113504/1‑150 (MQ=255)
                                                                                                tttACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCa                                                    <  5:1847111/150‑1 (MQ=255)
                                                                                                  tACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATct                                                 >  5:330763/1‑151 (MQ=255)
                                                                                                   aCGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATctt                                                <  6:850853/151‑1 (MQ=255)
                                                                                                     gAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATcttac                                              >  4:628708/1‑151 (MQ=255)
                                                                                                      aataatTTCTTGTGCGAGAAACAAGGTAACCTACCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCAACGCCTTATCCGGCCTACGTGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATcttact                                             <  7:1782760/151‑1 (MQ=255)
                                                                                                         aatTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATcttactta                                           <  3:1163928/150‑1 (MQ=255)
                                                                                                         aatTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtag                                          <  5:904778/151‑1 (MQ=255)
                                                                                                         aatTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGATAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtag                                          <  3:1522911/151‑1 (MQ=255)
                                                                                                          atTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATcttactt                                            <  6:23941/148‑1 (MQ=255)
                                                                                                          atTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagt                                         <  7:1644640/151‑1 (MQ=255)
                                                                                                           tttCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagt                                         >  5:386796/1‑150 (MQ=255)
                                                                                                           tttCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagt                                         >  1:1081243/1‑150 (MQ=255)
                                                                                                           tttCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagt                                         >  2:907336/1‑150 (MQ=255)
                                                                                                           tttCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagt                                         >  8:1561248/1‑150 (MQ=255)
                                                                                                            ttCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagta                                        <  8:1026074/150‑1 (MQ=255)
                                                                                                            ttCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagta                                        <  1:512304/150‑1 (MQ=255)
                                                                                                             tCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTtagtag                                       <  1:1725476/150‑1 (MQ=255)
                                                                                                               ttGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGCGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAgctg                                    >  8:1621523/1‑151 (MQ=255)
                                                                                                                tgtgCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAgct                                     <  8:615565/149‑1 (MQ=255)
                                                                                                                tgtgCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAgc                                      <  1:203962/148‑1 (MQ=255)
                                                                                                                  tgCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGctgctg                                 >  2:1011600/1‑151 (MQ=255)
                                                                                                                   gagagaAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTgg                                <  6:1303701/149‑1 (MQ=255)
                                                                                                                   gCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTgg                                <  6:62916/151‑1 (MQ=255)
                                                                                                                     gagaAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGgcg                              >  3:878448/1‑151 (MQ=255)
                                                                                                                     gagaAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCACGCGTAAAGCATCTTACTTAGTAGCTGCTGgcg                              >  3:52030/1‑151 (MQ=255)
                                                                                                                      agaAACAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGgcgc                             >  2:1621062/1‑151 (MQ=255)
                                                                                                                        aaaCAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCTGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGctc                           >  3:1804687/1‑151 (MQ=255)
                                                                                                                        aaaCAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGctc                           <  2:1542944/151‑1 (MQ=255)
                                                                                                                         aaCAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGctct                          >  2:1304799/1‑151 (MQ=255)
                                                                                                                         aaCAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGctct                          <  8:300397/151‑1 (MQ=255)
                                                                                                                           cAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCtt                         >  6:982727/1‑150 (MQ=255)
                                                                                                                           cAAGGTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCtt                         >  8:452970/1‑150 (MQ=255)
                                                                                                                              ggTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCg                     <  6:1417547/151‑1 (MQ=255)
                                                                                                                              ggTAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCg                     <  2:509331/151‑1 (MQ=255)
                                                                                                                                tAACCTCCCGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTAACGTc                   >  7:456457/1‑151 (MQ=255)
                                                                                                                                      cccGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGt              <  1:1200137/150‑1 (MQ=255)
                                                                                                                                      cccGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTg             >  5:458901/1‑151 (MQ=255)
                                                                                                                                      cccGAGGAAAAGAGAATGAGGCGGCGTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGTTACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTg             <  1:848707/151‑1 (MQ=255)
                                                                                                                                       ccGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTg             >  6:1588822/1‑150 (MQ=255)
                                                                                                                                       ccGAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGa            <  6:492046/151‑1 (MQ=255)
                                                                                                                                         gAGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGAcc          >  4:881902/1‑151 (MQ=255)
                                                                                                                                          aGGAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGAccc         >  5:1402440/1‑151 (MQ=255)
                                                                                                                                            gAAAAGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGAccc         >  8:115855/1‑149 (MQ=255)
                                                                                                                                             aatagagAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGACCCAg       <  6:554787/147‑1 (MQ=255)
                                                                                                                                             aaaaGAGAATGAGGCGGCTTTTTCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTATGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTAACTTAGTAGCTGCTGGCGCTCTTACCGTCGCAGTGACCCAgc      >  6:1630040/1‑151 (MQ=255)
                                                                                                                                             aaaaGAGAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGACCCAgc      <  4:1209083/151‑1 (MQ=255)
                                                                                                                                                agagAATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTACTAGCTGCTGGCGCTCTTACCGTCGCCGTGACCCAGCGCt   >  8:82314/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGGGCGGCAATCTGCCCTTCGGAGATAGTGTTGTCGCGAATACCGTTGATAAAGAAACGAATTTCTTCATCGCTCAGCGCATGACCATCACGTTTTTTACGAATAATTTCTTGTGCGAGAAACAAGGTAACCTCCCGAGGAAAAGAGTATGAGGCGGCTTTTGCCGGATGCGTCTACGCCTTATCCGGCCTACGGGTAGATACGTCTCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAAGCATCTTACTTAGTAGCTGCTGGCGCTCTTACCGTCGCCGTGACCCAGCGCTT  >  DH1/3894675‑3894970

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