BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	613,613	0	T	G	28.3%	193.9 / 39.5	201	intergenic (+474/+549)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base T (42/95); new base G (4/50); total (50/151)
Fisher's exact test for biased strand distribution p-value = 5.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

GGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATAC  >  DH1/613463‑613740

ggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGt                                                                                                                                 <  2:487983/151‑1 (MQ=255)

   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                              <  1:212875/151‑1 (MQ=255)

    ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                                             <  3:1659666/151‑1 (MQ=255)

    ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                                             <  2:1140906/151‑1 (MQ=255)

     gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacg                                                                                                                            <  2:511630/151‑1 (MQ=255)

       tgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                          >  7:832649/1‑151 (MQ=255)

       ggctTGAGCCAGTGAGCGATTGCTGGGCTAGATGAATGACTGTCCACGACAGAACCAGGCTTATCGGTCAGTTTCACCTGATTTACGTAACAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                          <  4:1214785/150‑1 (MQ=255)

       ggctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                          <  2:679896/150‑1 (MQ=255)

        gctTGAGCCTGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACATAAGCCTGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                                         <  2:53207/151‑1 (MQ=255)

        gctTGAGCCAGTGAGCGATTGCTGGGCTAGATGAATGACTTTCCACTACAGAACCAGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                          <  8:617372/150‑1 (MQ=255)

        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacga                                                                                                                           <  7:869881/149‑1 (MQ=255)

          ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAAACCCCATAATAAATCAGGTTATCCAGATATAAATAAAAACCA‑GATTCG‑AAAGATTTTTCCTATAG‑GAGTAGCAGAAAGATGAATGACTATCAactaaaa                                                                                                                        >  2:285296/1‑145 (MQ=255)

            gAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCtat                                                                                                                     <  1:426518/151‑1 (MQ=255)

             aGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCATTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCtata                                                                                                                    <  3:1308921/151‑1 (MQ=255)

              gCCAGTGAGCGATTGCTGGCCTAGATTAATGACTGTCCACGACAGAACCCGTCTTATCGGTAATTTTCACCTGATTTACTTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTT‑TAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATAc                                                                                                                   <  3:1165014/151‑1 (MQ=21)

              gCCAGTGAGCGATTGCTGGCCTAGATTAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATAc                                                                                                                   <  1:841454/151‑1 (MQ=255)

              gCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGCTTATCGGTCAGTTTCACCTGATTTACCTAAAAACCC‑CCATCG‑GAGGGGATGACCTATGA‑GAGCGGCTCAAAGAACAAAGACTGTGCACCACGCTAGAc                                                                                                                   >  1:1825248/1‑151 (MQ=19)

                  gTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGTTTTTTGCTTTTG‑GAGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCaa                                                                                                               <  3:955903/151‑1 (MQ=255)

                   tGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGCTAATAGCTCACTTTCCCCAGATTTAACTAAAAACCA‑GCTTAC‑GCAGGATTTTGGTTTTG‑GAGGGGCAGAAAGATGCATCACTGTCAACGCCGCTATACCCaa                                                                                                               >  7:1485446/1‑150 (MQ=255)

                   ggAGCGATTGCTGGCCTAGATTAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCaaa                                                                                                              <  1:1069445/150‑1 (MQ=255)

                    ttgcttttGCTGGCCTAGATTAATGACTGTCCCGCACAGAACCCGGCTTATCGGTCTGTTGCACGTGATTTACGTAAAACCCC‑GGTTCG‑GCGTGTTTTTGCTTTTG‑GCGGGGCATAAAGATGATTGAGTGTCCACGACGCTATACCCaaa                                                                                                              <  3:85815/144‑1 (MQ=255)

                    gAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATGTTAATAAAAAACG‑GCTGAGCGCCGGTCTTAG‑TGATG‑GAGGGGGTGAGAGATGAATTACTGGAAAATACACTATACCCaaa                                                                                                              >  5:819763/1‑150 (MQ=255)

                    gAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCCGCTTATCGGTCAGTTTAACCGGAATTACATAAAAAAAC‑GCATAG‑GCGGATTTTATCATTTT‑AACTGTCAGAATTATGAATTACTGTAAACGAAACTATACACaaa                                                                                                              >  6:54548/1‑150 (MQ=255)

                     cgCGATTGCTGGCGTAGCTGAATGACTGGCCCCCACAGAACCCGGCTTAGCGGTCAGTTTCACGTGATTTACGTAAAACCCC‑GGTTCG‑GGGGGTTTTTGCTTTTG‑CAGGGGAACAAAGATTAATGCCTGTCCACGACGCTATACCCAaaag                                                                                                            <  4:1635752/150‑1 (MQ=255)

                       cGATTGCGGGCCTAGATTAATGACTGGCCACCACAGAACCCGGCTTATCGGTCAGTTGCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                                          <  1:1408986/151‑1 (MQ=255)

                        gATGGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTTTACCCAaaagaaa                                                                                                         <  7:1762294/151‑1 (MQ=255)

                        gATGGCGGGCCTAGATCATCCCCTGCCCCCCACACACCCGGGGTTATCGGTAAGTTCCACCTTATTTACCTAAAACCCC‑GGTGCG‑GGGGTTTTTTTTTTTTG‑GCGGGCCAGAAAGATGAAGGACTGCCCACGACGCTATACCCAaaagaaa                                                                                                         <  7:90992/151‑1 (MQ=255)

                         aTTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                                         >  6:1698161/1‑134 (MQ=255)

                          ttgggggcgtattttaatgaCTGCCCACCCCACAACCCGGGTTTTCGGTCTTTTTCACGTTATTTACTTAAAAAGGC‑GCTGGG‑GGGGTTTTTTGTTTTTG‑TTGGGGCACAAAGTTTAATGAGTGTCCACGACGCTATACCCAaaagaaag                                                                                                        <  5:93645/135‑1 (MQ=255)

                          ttGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCATAAAAATCAATGACTGTCCACGAAGCTATACCCAaaataaatc                                                                                                       >  2:743048/1‑149 (MQ=18)

                           tGCTGGCCTAGATGAATGACGGGCCACGACACAACCCGGCTTATCGGTCAGTTTCACCTTTTTTACGTAAAACCCC‑GCTGCG‑GGGGGTTTTTGCTTTTG‑GAGGGGCACAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGc                                                                                                       <  5:1296358/150‑1 (MQ=17)

                           ggCTGGCCTAGATGAATGACTGCCCACGACAGAACCCGGCTTATGGGTCTGTTTCACCTGATTTACTTAAAAACCC‑GCTTCG‑GGGGTTTTTTGTTTTTG‑GCGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCg                                                                                                      <  8:1127364/150‑1 (MQ=18)

                              gggCCTATATGAATGACTGTCCACGACATAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                                   <  4:417845/150‑1 (MQ=255)

                              gggCCTAGATGAATCACTTTCCACGACAGAACCCGGCTTAGCGGTCTGTTTCACCTGATTTACGTAAAAACCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑GTGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                                   <  6:1317738/150‑1 (MQ=17)

                                gCCTCGATGAATCACAGCCCACGACATAACCCGGGTTAGCGGTAATTTTCACGTGATTTCCGTAAAAACCC‑GGTTGG‑GCGGTTTTTTGTTTTTG‑GTGGGGCAGAAAGATGAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTa                                                                                                 <  3:1801935/151‑1 (MQ=14)

                                  cgggcggcgagaCTATCCACAACAGGACCCGGGTTATCGGTCTGTTTCACGGGGGTTACTTAAAACCCG‑GCTGCG‑GCGTGTTTTTGGTTTTG‑TTGGGGCTGAAAGATGAATGACTGTCCACGACACTATACCCAAAAGAAAGCGGCTTATc                                                                                               <  3:272273/141‑1 (MQ=2)

                                   tAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAATTAAAAAAGA‑GAATAG‑TAGGGTTTTTGATTTTG‑TATTTGAATAAATATGAATGACTGTCCACTACTATATAAAAAAAAGAAAGCGGCTTATCg                                                                                              >  8:408114/1‑151 (MQ=12)

                                    aGATGAATGACTGTCCACGACAGAACACGGCTTATCGGTCAGTTTCGCCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAGAGCGGCTGAGCgg                                                                                             >  1:561082/1‑151 (MQ=255)

                                      attaatgaCTGTCCACTACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAGTTAAAAACCC‑GCTTCG‑GCGGTTTTTTGCTTTTG‑TAGGGTCATAAAGATTAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                           <  4:17517/148‑1 (MQ=17)

                                      agtaatgaCTGACCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGGTTTTG‑GCGGGGCAGAAAGATGAAGGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                           <  5:1068972/148‑1 (MQ=14)

                                       tgaatgaCTGTCCACGACAGAAACCGGCTAATCGGTCAGTTTCACCTGAATTACGTAAAAACCA‑GCTTCG‑TCGGTTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTCACAACGAATCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                           >  2:1292586/1‑150 (MQ=14)

                                        gaatgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                          <  1:133847/150‑1 (MQ=255)

                                        gaatgaCTGCCCACCACACAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAAACCC‑TCTTCG‑GCGGGTTTTTGCTTTTG‑GTGGGGCCCAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                          <  2:1637596/150‑1 (MQ=14)

                                        gaatcaggggccaccaCACAACCGGGCTTAGGGGCCATTTGCACGTG‑TTTAACTAAAAACCCGGGTTGG‑GGGGTTTTTTGTTTTTG‑GCGGGCCAAAAACTTCAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                          <  7:1639612/140‑1 (MQ=2)

                                        aaattcgtttaccaaaCATAACCCGTCTTATCGGTAATTTTAACCTGTTTTAGTTAAAAACCC‑GATTTT‑GCGGGTTTTTGTTTTTG‑TTGGGGCTTAAATATGATTGACTTTCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                         <  8:841009/136‑1 (MQ=2)

                                         aatcaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                         <  6:913444/150‑1 (MQ=255)

                                          atgaCTCGCGACGACAGAACGAGGTTTATGGGAAAGTTTCACTTCTTTTCCGTAAAAACCC‑GCTTCG‑GGGTGTTTTTGCTTTTG‑TTTGGGCATTAATATTAATCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGtt                                                                                       <  4:728136/151‑1 (MQ=2)

                                            gtgtatcgcaaaCAGAACCAGTCTTTTGGTGCAATTTCACGTGCTTAAGGTAAAACGGC‑TCTGCG‑GCGTGTTTTTGGTTTTG‑ATGGTGGCCAAATTTGAATGACTTGCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                      <  1:1547683/139‑1 (MQ=2)

                                            gaCTGCCCACTACAGAACCCGGCTTAGCGGGCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCATAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTc                                                                                     <  5:1199301/151‑1 (MQ=255)

                                             aCTGTCCACGACTTAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAACCCG‑GCTTCG‑GCGGGTTTTTGTTTTTG‑TAGGGGCTCATATATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCa                                                                                    <  6:1787363/151‑1 (MQ=14)

                                              cTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTACA‑TCGGAATATAGAATTTT‑AAGGGAAAGAAATATCAATAACTGTACAACACGCTATACCCAAAATAAATAAAAAAATCGGTCAGTTTCAc                                                                                   >  2:927860/1‑151 (MQ=9)

                                               cctcctttaCAGAACCAGTCTTATCTGTCAGTTTCACTGTCTTTAGTTAAACCCGG‑GCTGCG‑GCGGGTTTTTGGTTTTG‑TTGGGTCCGAAATATGATTTACTGCCCACGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCaac                                                                                  <  8:1279177/143‑3 (MQ=2)

                                                tcccaccaCACAACCCGGCTTATCGCTCAGTTTCACGTTATTTACTTAAAAACCC‑TCTTCG‑GGGGGTTTTTGTTTTTG‑GCGGGGCTTAAATATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCAcc                                                                                  <  4:1832011/148‑1 (MQ=14)

                                                   cACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGGTTTTG‑GAGCGGCAGAAACATCAATCAATGGACACGAAGCTATACCCAAAAGAAAGCGGCTTATCGGTAAGTTTCACCGGGt                                                                              >  5:1634088/1‑151 (MQ=14)

                                                    ccGACACAACCCGGCTTATCGGTCATTTTCCCCTGATTTAGTTAAAAACCG‑GCTTCG‑GCGGGTTTTTGTTTTTG‑GTGGGGCTCAACTATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGtt                                                                             <  2:1490641/150‑1 (MQ=14)

                                                     caacacAACCCGGCTTATGGGCAATTTGCCCGTCATTTACTTAAAACCGG‑TCTGCG‑GCGGGTTTTTGTTTTTG‑TTGGGGCACAATTATGAATGACTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                            <  2:603763/149‑1 (MQ=9)

                                                       aCAGAACCCGGCTTATCGGTCATTTTCACGTGATTTAGGTAAAAACCC‑GCTGCG‑GCGGGTTTTTGGTTTTG‑GCGGGGCATAAAGATGAATGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                            <  7:765919/149‑1 (MQ=14)

                                                       aCAGAACCCGGCTTATCGGTCATTTGCACCTGATTTAGTTAAAAACCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑GCGGGGCCCAAACATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                            <  4:1663201/149‑1 (MQ=14)

                                                       aCAGAACCCGGATTAACGTGCATTTTCACGTG‑TTTTCGTTATAACCCGGCTTCG‑GCGGGTTTTTGGTTTTG‑GAGGGGGATAAAGATGAATGACTGGCCACGTCGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTAc                                                                          <  5:1188099/151‑1 (MQ=14)

                                                          gAACCCGGCTTATGGGGCATTTTCACGTCATTTAGGAAAAACCGG‑GCTGCG‑GCGGTTTTTTGTTTTTG‑GTGGGGCCCAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTa                                                                       <  1:707200/151‑1 (MQ=14)

                                                              ccGGATTATGGGTCAGTTTCCCGTGATTTACTAAAAAACCC‑CCTTCG‑GGGGGTTTTTGGTTTTG‑CTGGCGCAGAAAGATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                                   <  7:587333/151‑1 (MQ=11)

                                                                ggCTTATCGGTCAGTTTCACCTGATTACCGTAAAACCCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑CAGGCGCAGAAAGATGAAGGACTGGCCACGACGCTATACCCAAACGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAcc                                                                 <  7:1845602/151‑1 (MQ=11)

                                                                 gcgtATGGGACAGTTGCCCGTGATTTACGAAACACCCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑CCGGGGCCCAAACATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAccc                                                                <  1:1501531/148‑1 (MQ=11)

                                                                 gCTTATCGGTAATTTTAACGTTATTTAGGTAAAAACCG‑GGTTGG‑GGGGGTTTTTGCTTTTG‑CTGGGGCATAATTTTTAATGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAccc                                                                <  1:1365091/151‑1 (MQ=11)

                                                                 gCTTATCGGGCATTTTCACCTGATTTCCTTAAAACCCG‑CCTTCG‑CCGGGTTTTTGTTTTTG‑TTGGCGCCAAACCTTGAAGGCCTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAccc                                                                <  2:1174737/151‑1 (MQ=2)

                                                                    taccGGCCAGTTTCCCCTGATTTACTTATTAATCC‑TCTTCG‑GGGGTTTTTTGATTTTT‑TATGTTATTATATATCTATTACTGTCCACGACGCTATTCCGAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGc                                                              <  4:878118/147‑1 (MQ=2)

                                                                     atgggaatgTTTCCCCTCATTACCTTAAAAACCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑CCGGGGCCTAAAGTTGATTGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCtt                                                            <  6:666360/143‑1 (MQ=11)

                                                                     aaggggaaGTTGCACCTGATTTAGGTAACCCCCG‑GCTGCG‑GCGGGTTTTTGGTTTTG‑GCGGGGCACAACCATGAATGCCTGTCCACGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCtt                                                            <  3:101325/144‑1 (MQ=11)

                                                                       ggggaattttgcccTGATTTACGTAAAAACCG‑CGTGCG‑GGGGGTTTTTGGTTTTG‑GTGGGGCAAAAATTTGAATGCCTGGCCACGACCCTATACCCAAAAGAAAGCGGCGTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                                          <  7:1591403/139‑1 (MQ=11)

                                                                       cGGACATTTTCCCCTGATTTAGTTAAAACCCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GTGGGGCTTAATGATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTc                                                           <  8:486569/150‑1 (MQ=11)

                                                                          caattttCACGTGATTTCCGAACCAACGG‑TCTTAT‑GCGGGTTTTTGTTTTTG‑TTGGTGGTAAATCTTGAATGACTTTACTAGTCTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggcg                                                       <  4:1428124/147‑1 (MQ=2)

                                                                            actttcggttGATTTACTTAAAAACCG‑GCTGCG‑GCGGGTTTTAGGTTTTG‑CAGGGGGACAACGATGAATGACTTGCCACGACGCTATACCCAACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCggg                                                     <  1:96640/142‑1 (MQ=11)

                                                                               ttCACCTTATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTGTACGTAATAACCCGGTGCGGCGGGttt                                                  >  7:1081806/1‑151 (MQ=11)

                                                                               ttCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGCTGTACGTAAAAACCCGGTGCGGCGGGttt                                                  >  7:374853/1‑151 (MQ=255)

                                                                               tgccccTGATTTACGAAAAAACCC‑GCTTCG‑GCGGGTTTTTGGTTTTG‑CAGGGCCAGAAAGATGAAGCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGG‑TTACGTAAAAACCCGCTTCGGCGGGttt                                                  <  3:416740/146‑1 (MQ=11)

                                                                                tCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTTACCTGATTTACGTAATAACCCGGTGCGGCGGGtttg                                                 >  3:575940/1‑150 (MQ=255)

                                                                                 cccctcagtcaCGAAAAACCCC‑CCTGCG‑GGGGGTTTTTGCTTTTG‑CCGGGGCCCTAACATCACGCACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGGCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                                <  8:1692392/141‑1 (MQ=2)

                                                                                 ccccTGATTTACGTAAAAACCC‑GCTGCG‑GCGGGTTTTTGGTTTTG‑GTGGCGCATAATCATGATTGACTTTCCACTACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtttt                                                 <  2:1772116/148‑1 (MQ=11)

                                                                                  aCCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAAACCGGTTCGGCGGGtttggg                                               >  5:210672/1‑148 (MQ=9)

                                                                                      gATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCCCGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACATGGTGTACGTAATAAACCGCTTCGGCGGGTTTGTGCttg                                           >  7:1138773/1‑150 (MQ=11)

                                                                                      gATTAACGTAAAACCCC‑GCTGCG‑GCGGGTTTTTGCTGTTG‑GCGGGGCACAAACATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtt                                            <  3:1522859/150‑1 (MQ=11)

                                                                                      gATTAACGTAAAAACCC‑CCTTCG‑GGGTGTTATTGATTTTG‑CAGTGCATTAACTATGAAGGACTTTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtt                                            <  6:920784/150‑1 (MQ=11)

                                                                                       attaacttattaaCCG‑GCTTAC‑GGGGGTTTTTGTTTTTG‑GCGCTGTATATATATCATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGCTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                                          <  2:1443139/140‑1 (MQ=2)

                                                                                       aTTTCCGTAAAACCCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑GCGGCGCCGAAACATGAAGGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                                          <  5:80562/151‑1 (MQ=11)

                                                                                          tccgaaaaacCCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑GCGGGGCCCAAACATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTgg                                        <  3:460829/145‑1 (MQ=11)

                                                                                           aCGTAAAACCCC‑GCTGCG‑GCGGGTTTTTGCTTTGG‑TCGGGGCAGAACGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                                       <  5:1517228/150‑1 (MQ=11)

                                                                                             gTAAAACCCC‑CCTTGC‑GGGGGTTTTTGATTTGG‑GCGCGGCAGAACTATGATGGCCTCGCCACCACGCTATACCCAACAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                                    <  1:740746/151‑1 (MQ=11)

                                                                                              tAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGCTGTACGTAAAAAACCGTTGCGGCGGGTTTCTGCTTGTGGAgggg                                   >  3:570002/1‑151 (MQ=11)

                                                                                               aaccccgc‑GCTGCG‑GCGGGTTTTTGCTTTTG‑CCGGGGCACAACACTGAAGGCCTCGCCACGACGCTATCCCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGc                                  <  7:901462/144‑1 (MQ=11)

                                                                                               aaacaCCC‑GCTGCG‑GCGGGTTTTTGCTTTTG‑ACGGGGCCCACACATGAATGACTGGCCACGACGCTATACCCAACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGc                                  <  1:315614/151‑1 (MQ=11)

                                                                                                aacccgc‑GCTGCC‑GCGGGTTTTTGGTCTTG‑CCGGGGCCCAACCATGATTGCCTGGCCCCGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCa                                 <  1:342220/145‑1 (MQ=11)

                                                                                                aaaaaccgttctgctgcggGTTTTTGTTTTTG‑TTGGTGTTTAATTTTGAAGGACTCTCCTAGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAgggg                                   <  4:263076/136‑1 (MQ=2)

                                                                                                aaaaCCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGAAGGATCACTGGGCACGACGCGATACCCAAAACAAACGGGCTTATCGGTGAGTTTCCCCGGGTTTACGTAAATCCCGGCTTCGGGGGGTGTTGGCGTGGGGAgggggc                                 >  5:580922/1‑149 (MQ=2)

                                                                                                  ccccgttctgag‑gcggGTTTTTGCTTTTGTGGGTGGCA‑TATCTTGAATGCCTCGCCACGACGCTATACCCTACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGTGGCa                                 <  6:1534093/139‑1 (MQ=2)

                                                                                                    ccc‑cCTGCG‑CCGGGTTTTTGCTTTTG‑CAGGGGCACAACGATGAAGGACTCTCCACCACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                             <  3:749749/151‑1 (MQ=11)

                                                                                                     ccgcctgcg‑gcggGTTTTTGGGTTTG‑CCGGGGCCCCAACTTGCTGGACTGGCCACCACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                             <  5:1779263/144‑1 (MQ=11)

                                                                                                          ggtcgggttttttttGGTTTTG‑GTGGGGCTTAAATATGAATGCGTTCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg                         <  7:101438/141‑1 (MQ=11)

                                                                                                                  ggTTTTTGCTTTTG‑CAGGGGCAGAAACATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                 <  5:236931/151‑1 (MQ=25)

                                                                                                                         ggttttG‑CGGGGGGGGAAAGATTAATGCGTTGCGAAGACTCTATACCCAAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACg           <  1:1592782/148‑1 (MQ=11)

                                                                                                                                  aGGGGCAGAAAGATGAATGAATGTCCAAGACGCTATAACAAAAAGAAAGCGGATTATCGATCAGTTTAACATGGTTTACGTAAAAACAATATTATAAGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACGGTcggcagcgcgtggg  >  6:466763/1‑138 (MQ=2)

GGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCC‑GCTTCG‑GCGGGTTTTTGCTTTTG‑GAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATAC  >  DH1/613463‑613740

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