BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	613,597	0	G	C	26.3%	410.3 / 60.1	280	intergenic (+458/+565)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base G (77/102); new base C (4/60); total (86/194)
Fisher's exact test for biased strand distribution p-value = 9.91e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGT  >  DH1/613447‑613726

gCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAg                                                                                                                                   <  8:744693/151‑1 (MQ=255)

gCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAg                                                                                                                                   <  7:549601/151‑1 (MQ=255)

 cccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                                                                                                                  <  7:1098657/151‑1 (MQ=255)

     tACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACGTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                 >  3:1066238/1‑148 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat                                                                                                                             <  7:1112315/151‑1 (MQ=255)

        tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAATAAAC‑C‑G‑ATTCGGCGGCTTTGAGCTTTTGGAGTGCCAGTAAGaac                                                                                                                            >  8:1194952/1‑148 (MQ=37)

        tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGa                                                                                                                              >  6:187387/1‑148 (MQ=255)

         gAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga                                                                                                                           <  5:386751/150‑1 (MQ=255)

         gAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga                                                                                                                           <  4:1405195/150‑1 (MQ=255)

           aCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg                                                                                                                        <  7:891070/151‑1 (MQ=255)

           aCCCGGGTAGGCGGCTTGAGCCAGTGAGCGATTGCGGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg                                                                                                                        <  8:619811/151‑1 (MQ=255)

            ggcGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatga                                                                                                                       <  8:1407805/149‑1 (MQ=255)

            cccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatga                                                                                                                       <  7:604586/151‑1 (MQ=255)

             ccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCCGCTTATCAGTCAGTTTAACCTGAATTACGTAAAAAAC‑C‑G‑CATCAGCGGGATATACATTTTGGACGGGCAGAAAGatgaaaga                                                                                                                       >  2:1511083/1‑150 (MQ=25)

              gggttttggtgtttGAGCGTGGGATCGATTGATGGACTAGAGCAGACACTATACACCAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAATGACt                                                                                                                     <  1:1808549/139‑1 (MQ=11)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                    >  7:1269106/1‑151 (MQ=255)

                  tAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGATATATGAATGGCTGTcc                                                                                                                 >  3:623332/1‑151 (MQ=255)

                   aGGCTGCTTGAGCCAGTGTGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGCTCTGTGCTTTTCGAGCGTCAGAAAGatgcatgcatgg                                                                                                                   >  3:886890/1‑143 (MQ=25)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAGACGGCTTATCGGTCAGTTTCACCGGATTTACGTAAAAACC‑C‑G‑CTTCGGCCGGTTTTTGCTTGTGGAGGGGCAGATAGATGAATGACTGTccc                                                                                                                >  5:1345080/1‑150 (MQ=21)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                <  4:1333526/151‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                <  4:1903502/151‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAACCCGCTTATCGGTCAGTTTAACCTAAATAAGCAAAAAAAC‑C‑T‑CTACGGAAAGATTTTGATTTTGGTTGTAGAGAAAGATGAATGACTgaac                                                                                                                 >  2:1127958/1‑147 (MQ=255)

                    ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                               <  7:246764/151‑1 (MQ=255)

                     tgttgttGAGCCGGAGACCGATTGCTGGCCTCGACAAGTGACTGTCCCCTACAGAACGCGGCTTATCGGTCAGTTTCACCTGCTTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTTCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                               <  1:553724/145‑1 (MQ=14)

                     gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacg                                                                                                              <  1:1678337/151‑1 (MQ=255)

                      ctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATATAGGAAAAAACA‑C‑G‑ATTAGGCAGATTTTTTCTTATCGAGTGTAAGAAACATGAATGACAGtaaa                                                                                                                >  6:1662155/1‑145 (MQ=25)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                           <  7:518909/151‑1 (MQ=255)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                            <  5:1470247/150‑1 (MQ=255)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                            <  3:812409/150‑1 (MQ=255)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGCGAGTAAGagcaatgcgatgacaccccg                                                                                                           >  7:1758482/1‑138 (MQ=25)

                         ctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAACCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTAGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGtaaaagaaa                                                                                                           >  2:792542/1‑142 (MQ=255)

                            gAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATATACGTAAAAACC‑C‑T‑CATAGTCTGGAAATTGCTTTTGGACGGTCAGAAACATGAATGACTATcaaaaactaaat                                                                                                       >  6:572587/1‑140 (MQ=21)

                            gAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGCCTTATCGGTAAGTTTAACCTGAATTACGTAAAAAAC‑C‑T‑CTTCGGCGGGTTATTTATTTTTTATAGGCAGAAAGATGAATGACTGTCCacgaccgaat                                                                                                       >  6:155226/1‑146 (MQ=21)

                              gCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATAc                                                                                                     <  5:1128675/151‑1 (MQ=255)

                              gCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAAACCCATTATCAAACAGATTCACCTCATTTACGTAAAAAAC‑C‑G‑CTTCGTCGGCTTTTTGCTTTTGGATGGGCAGAAAGATGAATTACTATCCACGACACTATAc                                                                                                     >  8:1746754/1‑151 (MQ=18)

                                    gAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCaaaa                                                                                               <  2:158449/151‑1 (MQ=255)

                                     aGCGATTGCGGGCCTAGATGAATCCCTGTCCACGACAGAACCCGGCTTAGCGGGCATTTTCCCCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGCGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaag                                                                                              <  6:454486/151‑1 (MQ=18)

                                      gCGATTGCTGCCCTAGATGAATGCCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaaga                                                                                             <  4:250458/151‑1 (MQ=255)

                                       ggATTGCTGCCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                            <  6:1483549/150‑1 (MQ=255)

                                       cGATTGCGGGCCTAGATGAAGGACTGCCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTTATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGCGGGGCACAAAGATCAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                            <  8:848175/151‑1 (MQ=18)

                                        gttnGCTGGCCTNTATGAATTACTGGCCACGCCAGAACCCGGGTTATCGGTCAGTTTCACGTGATTTACGTAAAAACC‑C‑G‑GTTCGGCGGTTTTTTGCTTTTGCAGGGGCAGAAATATGATTGACTGTCCACGACGCTATACCCAaaagaaa                                                                                           <  5:1774007/149‑1 (MQ=255)

                                        gATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGAAAAAAGATGAAAGACTGTACACGAATCTATAACAAaaagaa                                                                                            >  8:529351/1‑150 (MQ=18)

                                        gATTGCTGCCCTAGATGAATGACTGCCCACGACAGAACCCGGCTTATCGGTCTGTTTCACCTGATTTACTTAAAAACC‑G‑G‑CTTCGGGGGTTTTTTGATTTTGGTGGGGCATAAATATGAATTACTGTCCACGACGCTATACCCAaaagaaa                                                                                           <  4:590400/151‑1 (MQ=255)

                                        gATGGCTGGCCTAGATGAATGACTGTCCACGACACAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAACCC‑C‑G‑CGTCGGCGGGTTTTTGCTTTTGGAGGGCCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaaa                                                                                           <  3:268321/151‑1 (MQ=255)

                                         aTTGCTGGCCTAGATCAATGACTCTCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTTCTTTTGGAGGGGCAGAAAGatgaatg                                                                                                                        <  3:1509192/121‑1 (MQ=255)

                                           tGCTGGCCTAGATGAATGACTGGCCACCACAGAACCCGGCTTATGGGTCAGTTTCA‑CTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCg                                                                                        <  1:1047407/150‑1 (MQ=255)

                                            gCTGGCCTAGTTGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGGGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCgg                                                                                       <  5:865488/151‑1 (MQ=255)

                                            gCTGGCCTAGATGAATGACTGTCCACGACATAACCCGCCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGTAGGGGCAGAAATATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCg                                                                                        <  4:441098/150‑1 (MQ=255)

                                               ggCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACGTGATTTACGTAAAAACA‑C‑G‑CGTCGGCGGGGTTTTCCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGACAGCGGCtt                                                                                    >  7:422313/1‑151 (MQ=18)

                                                gCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTAAAGAAAAAAAC‑C‑G‑AATAGTAGGGAATTTGATTATGGACGTGAAGAAAGATCAGTGAAAGTCCAATAATAAATAAAAAAAAGAAAGCGGCtt                                                                                    >  6:737857/1‑150 (MQ=255)

                                                 ccTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTAGTTTTGTAGAGCCAGATAGAGCACTGCATGGCCACGACGAAATACCCAAAAGACAGCGGCtttt                                                                                  >  7:658685/1‑149 (MQ=17)

                                                  ccggccccgcaaCTCTCCCCAACCGAACCCGGCTTATGGTTCATTTTCACCTGATTTACGTAAAACCC‑C‑G‑GTTCGGCGGGTTTTTCCTTTTGGAGGGGCCGAAAGATGAATTACTGTCCACGACACTATACCCAAAAGAAAGCGGCTTATc                                                                                 <  7:1611157/140‑1 (MQ=11)

                                                  cTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCATAAAGATGAATGACTGTCCACGAAGCTATACGAAAAAGAAAGAGGAATATc                                                                                 >  4:751828/1‑151 (MQ=18)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGGTTATCGGTCAGTTTCACCTGATTTAGGTAAAAACC‑C‑G‑CTTCGGCGGTTTTTTGATTTTGGAGGGCCATAAAGATCAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCg                                                                                <  7:1605737/150‑1 (MQ=17)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCg                                                                                <  2:401097/150‑1 (MQ=255)

                                                    aGATGAATGACTGGCCACGACAGAACCCGGATTATCGGGCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                 <  5:1018953/149‑1 (MQ=255)

                                                      ttgaatgaCTGTCCACGACAGAACCCGTCTTATCGTTAATTTTCACCTGTTTTACGTAAAACCG‑G‑G‑CTTCGGCGGGTTTTTGTTTTTGGCGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGt                                                                              <  6:644552/149‑1 (MQ=17)

                                                      atgaatgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTGAACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                             <  5:1537742/151‑1 (MQ=255)

                                                          atgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTCTCCACGACCCTATACCCAAAAGAAAGCGGCTTATCGGTcggat                                                                         >  7:872220/1‑147 (MQ=255)

                                                          atgaCTCTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGtt                                                                         <  6:1736256/151‑1 (MQ=255)

                                                             ccTGACCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGCGGGGCCGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCa                                                                      <  1:1138654/150‑1 (MQ=255)

                                                             aCTTTCCACGACATAACCCTTCTTATGGGTAATTTTCACGTGATTTACTTAAAAACC‑G‑G‑CTTCGGCGGGTTTTTGTTTTTGTTGGGGCTTAAAGATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTc                                                                       <  2:818387/150‑1 (MQ=14)

                                                                gcccaccaCAGAACCGGCCTACGCGGTAATTTGCCCCTGATTTACTTAAAAACC‑G‑T‑GTTCGGCGGTTTTTTGTTTTTGTCGGGGCTTAATCATGAATGAGTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCt                                                                   <  8:242705/149‑1 (MQ=2)

                                                                gcccaccaCAGAACCCGGATTATGGGGAAGTTGCACGTCATTTACTTAAAACCC‑C‑C‑CTGCGGCGGTTTTTTGTTTTTGTTGGGGCACAAAGATGAATGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCAc                                                                     <  1:950653/147‑1 (MQ=9)

                                                                gTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCTCAAGGATGAAAGTCTGTGAACGACACTATACCCAAAACAAAGGGGCTGATCGGTCAGTTGCaagg                                                                   >  4:639145/1‑148 (MQ=14)

                                                                 tCCACTACAGAACCCGGCTTATCGGTAAGTTTCACCTGATTTACTTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCACAATTATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCAccg                                                                   <  4:875426/150‑2 (MQ=14)

                                                                 cccACGACACAACCCGGCGTATCGGGCAGTTTCACCTGATTTCCGTAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAACATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTg                                                                  <  5:967915/150‑1 (MQ=14)

                                                                         aGAACCGGGATTATGGTTAATTTTCACCTGATTT‑TGTTAATACC‑C‑GTCTGCGGCGGGTTTTTGGTTTTGTTGGTGCTTTATTTTGATTGAGTCGACACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGt                                                          <  2:1273548/151‑1 (MQ=9)

                                                                           accccGGATCATCGGCCACTTGCCCCTGATTTCCGTACCCACC‑G‑G‑CTGCGGCGGGTTTTTGCTTTTGGCGGGGCCCACCCATGAAGGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaa                                                        <  2:239849/149‑1 (MQ=2)

                                                                              ccgcattATCGGGCAGTGGCGCCTGATTT‑CGGAAAACCCGC‑G‑CTTCGGGGGGTTTTTGATTTTGTGGGGGACGAAGCATTAAGGACTGTCCACGACGCTATACCCAAAAGAAAGCCGCGTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                     <  3:497406/146‑1 (MQ=255)

                                                                                  cTTATCGGTCAGTTTCACCTGATTTACTTAAAAACC‑C‑G‑CAGCTGAGGGTTTTTGGGTTTGGAGGGGCAGAAATATGAATGACTGTCCACGACAATATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTTTGCGtttaagcgcg                                                 >  1:1588559/1‑142 (MQ=11)

                                                                                    tATCGGGAATTTTCCGCTGATTT‑CCTAAAAACC‑CGG‑CTGCGGCGGGTTTTTGCTTTTGCTGGGGCTTAATCATGAATGACTGTCCCCGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCg                                                 <  5:1942674/149‑1 (MQ=11)

                                                                                     aTCGGGAAGTTTCACCTGATTTCCGTAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGCAGGGGCACAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCt                                               <  1:779034/150‑1 (MQ=11)

                                                                                       gggTCATTTGCACTTGATTTCCGTAAAACCG‑G‑G‑CTTCGGGGGTTTTTTGTTTTTGTCGGGGATAACCTTTGATTGACTCTCCACTACGCTATACCCTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                            <  6:1374560/150‑1 (MQ=2)

                                                                                       cGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGAGTTACgaagagcacacgtctg                                            >  1:1861589/1‑136 (MQ=11)

                                                                                           cAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTTTACGTAATAACCCGCTGcgccgg                                        >  7:1469508/1‑151 (MQ=255)

                                                                                           aaGTTTCCCGTGATTAACGTAACACCG‑T‑T‑CTTCTGCGGTTTTTTGTTTTTGTAGGAGCATAATTATGATTGACTTTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggcg                                         <  8:1006439/149‑1 (MQ=2)

                                                                                              tttCCCCTGATTTACTTAATAACG‑T‑T‑CTTCTGGGGTTTTTTGTTTTTGTTGGTGTTTAATTTTGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGt                                      <  6:717270/150‑1 (MQ=2)

                                                                                                tCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTCCTTTTGGAGCCGCACAAAAAAGAATCACTGTCCACGACCCTATACCCAAAAGAAATCATCTAATCGGTCACTTTCACCAGATAAACAAAAAAACCCGCTTCGGCGGCtttt                                   >  8:314227/1‑151 (MQ=2)

                                                                                                ggcggtGATTTACGAACAAACG‑C‑T‑CTGCTGCGGGTTTTTGTTGTGGTAGGTGTTCAACTATGCAGGCCTTTCCACGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttt                                    <  4:1093279/145‑1 (MQ=2)

                                                                                                 gtctgtttttagtaaaaaaCC‑C‑T‑CTTCGGCGGGTTTTTGGTTTTGTCGGTGCTTTAATCTGATTGACTTTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtttt                                   <  2:20117/136‑1 (MQ=11)

                                                                                                  cgcTGATTTCGTTAAACACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCTGAAAGATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTg                                 <  2:432517/149‑1 (MQ=11)

                                                                                                   ccTCATTTACGAAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGGTTTTGCAGGCGCCCAAACATGAATGACTCCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGc                                <  7:1677592/151‑1 (MQ=11)

                                                                                                      gATTAACGTAAAAACG‑C‑G‑CTGCGGCGGGTTTTTGCTTTTGGAGGCGCCGCAAGATGAAGGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                               <  3:435031/149‑1 (MQ=11)

                                                                                                      caTTTACGTAACACCC‑C‑G‑CTGCGGCGGGTTTTTGCTTTTGGCGGCGCACAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                               <  5:590733/148‑1 (MQ=11)

                                                                                                      caTTTACGAAAAACCC‑C‑C‑CGGCGGCGGGTTTTTGCTTTTGCCGGGGCACAACCATGAAGGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtt                              <  7:1817929/149‑1 (MQ=11)

                                                                                                         ttagttataaccc‑c‑T‑CTGCGGCGGGTTTTTGGTTTTGGCGGGGCTCAATGATGATTGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTg                           <  4:1280034/139‑1 (MQ=11)

                                                                                                            cGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGAGCAATGACTGTACACGAAGCTATACCAAAAAGAAAGAGGCTTATCGGTCAGTGGCAACAGTTTAACGTAAAAACCCGCTTCGGCGGGTGTTTGCTTTTGGAgg                       >  4:294232/1‑151 (MQ=11)

                                                                                                             gTAAAACCC‑C‑G‑CTTCGGAGGGTTTTTGGTTTTGGTGGGGCAGATAGACGAATCACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                      <  3:752650/151‑1 (MQ=11)

                                                                                                             cgttaatcccggg‑CTGCGGCGGGTTTTTGGTGTTGGTGGGGCATAATCTTGATGGCCTGTCGACGACGCTATACCCATAAGTAAGCGGCTTATCGGTCAGTTTCAGCTGGTTTAGGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                         <  2:528708/138‑1 (MQ=2)

                                                                                                              tAAAACCG‑C‑G‑CTTCTGCGGGTTTTTGCTTTTGCAGGGGCTCAATTATGATTGACTCTACACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                      <  6:1163090/150‑1 (MQ=11)

                                                                                                                 aaaaaacggccTTCGGCGGGTTTTTGCTTTTGTAGGTGCTCAATCATGAATGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                      <  8:1619555/140‑1 (MQ=11)

                                                                                                                  accc‑c‑G‑CTTCGCCGGGTTTTGGGTTTTGCGGGGGCCCAAACTTGAAGGCCTGGCCACGACGCTATACCCAAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                 <  7:803473/149‑1 (MQ=11)

                                                                                                                          ccTTCGGCGGGGTTTGGGTTTTGCGGGAGCGATCAGATGATGGCCTTGCCTCCACGCTCTACCCCAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat            <  5:1816707/150‑1 (MQ=2)

                                                                                                                           cTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCGCCGAAGATATACCGAAAGCAAACCGGGGTAAGTGAGAGTTTCACGTTGTTTACGTAAAAAACATATTACGAGGGTTTTTGCTTTTGGCGGGGCAGAAAGatg           >  1:511229/1‑151 (MQ=2)

                                                                                                                           cTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAATGTCCACGAACAGATAGCCAAAAGAAAGCGGCATATAGGTAAGGTTCACCTGGGTTAAAAAAAAAACCGCTTCGGCGGGTGTTTCCTTTTGGAGGGGCAGAATGatc           >  5:468270/1‑150 (MQ=11)

                                                                                                                              cgccggGTTTTTGGTTTTGCTGGCGCAAAAATATGATTGCCTCGCCAAGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat        <  1:130669/148‑1 (MQ=11)

                                                                                                                              ccgcggGTTTTTGTTTTTGCCGGTGGACCACCATGAAGGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat        <  1:21801/149‑1 (MQ=11)

                                                                                                                               cgcggGTTTTTGGTGTTGGCGGCGGCCAACTATGATGGACTCTCCACGACGCTCTACCCCAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg       <  3:1239203/150‑1 (MQ=11)

                                                                                                                                gcggGTTTTTGTTTTTGGCGGGGCACAAACATTCTGGACTGGCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatga      <  7:1426144/151‑1 (MQ=11)

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGT  >  DH1/613447‑613726

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