BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	613,597	0	G	C	22.9%	493.4 / 48.3	309	intergenic (+458/+565)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base G (79/110); new base C (3/53); total (90/219)
Fisher's exact test for biased strand distribution p-value = 4.03e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACC  >  DH1/613447‑613741

cccccgaCTGACCCGGGTTATGCGTCTTGAGCCAGCGAGAGATTGCTGCCCTAGACCACACACTCTACACAACAGAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAg                                                                                                                                                  <  1:1262225/145‑1 (MQ=2)

 cccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                                                                                                                                 <  3:65609/151‑1 (MQ=255)

  ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                                <  6:334025/151‑1 (MQ=255)

  ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                                <  3:1238846/151‑1 (MQ=255)

  ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                                <  2:911375/151‑1 (MQ=255)

  ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                                <  1:807160/151‑1 (MQ=255)

  ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTACTTTTGTAGGGACAGa                                                                                                                                                 >  8:863560/1‑150 (MQ=255)

   cGTACTGACCCCGGTTAGGCGGCTTGAGCCAGCGACCCCTTGCTGGCCTAGACGACACACTCTCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGaaa                                                                                                                                               <  1:1189411/151‑1 (MQ=17)

   cGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                                                                                                                                               <  1:526919/151‑1 (MQ=255)

    gTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                                >  2:815930/1‑149 (MQ=255)

    gTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAg                                                                                                                                              <  7:1477963/151‑1 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa                                                                                                                                             <  3:1095032/150‑1 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACCTAAAAAAC‑C‑G‑CTTCGGCGGATTTTTACTTTTGGATGGGCAGaa                                                                                                                                                >  8:828843/1‑147 (MQ=255)

        tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAAGTAAAAACA‑A‑G‑CTAAAGCGTGTTTTTCATATTGTAGGTGaagaaaga                                                                                                                                             >  8:152570/1‑148 (MQ=37)

          aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                                                                                                                                        <  6:1926582/151‑1 (MQ=255)

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

              cGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACt                                                                                                                                    <  7:817373/151‑1 (MQ=255)

              cGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACt                                                                                                                                    >  7:922843/1‑151 (MQ=255)

               gggTAGGGTGGTTGAGCGAGTGAGCGATTGCTGGGCTAGATGAATGACTGTCCACGACAGAACGCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                                   <  6:1329035/151‑1 (MQ=255)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGGATGCCtt                                                                                                                                   >  3:1944859/1‑150 (MQ=255)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                                   <  7:783511/151‑1 (MQ=255)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                                   <  1:478683/151‑1 (MQ=255)

                ggtttggttgttGAGGCAGTGAGCTATTGATGGACTAGATCACACACTCTACACCAAAGAAACCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTTCTTTTGGAGGGGCAGAAAGATGAATGACTGt                                                                                                                                  <  3:1644425/141‑1 (MQ=11)

                 gTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTc                                                                                                                                 <  2:1427642/151‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTcc                                                                                                                                <  6:1090500/150‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTcc                                                                                                                                <  7:762495/150‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                               <  3:900153/151‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑A‑G‑CTTAGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTcc                                                                                                                                >  5:798670/1‑150 (MQ=255)

                    ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                                              <  4:1563172/151‑1 (MQ=255)

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

                       tgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                                                                                                                                 <  7:940417/129‑1 (MQ=255)

                       tgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTCCTTTTGGAGGGGCTGAAAGATGAATGACTGTCCacatg                                                                                                                           >  5:1790981/1‑148 (MQ=255)

                       ggctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                           <  3:1728104/150‑1 (MQ=255)

                         ctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAACCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CGGGGGCGGGGTTTGGCTTTTGGAGGGGGAGAGAGATGGGTGGCTGTCGACGTCGc                                                                                                                         >  4:1214675/1‑151 (MQ=21)

                            gAGCCTGTGAGCGATTGCTGGCCTAGATGAATCACTGTCCACGACATAACCCGGCTTATCGGGCAGTTTCACCTGATTTACTTAAAAACC‑C‑G‑CTTGGGGGGGTTTTTGCTTTTGGAGGGGCACAAAGATCAATGACTGTCCACGACGCtat                                                                                                                      <  8:466293/151‑1 (MQ=21)

                             cgcgagtgacgtaTTGCTGGCCTAGATGAATGTCTGTCCCCGACAGAACCCGGCTTATCGTTCTGTTTCACCTTATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGTTTTTGGTGGGGCTTAAAGATGAATGACTGTCCACGACGCtata                                                                                                                     <  4:818543/139‑1 (MQ=18)

                             aGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑A‑G‑CTACGAAGAGATTATTGTTTTGCAAGCGAATAAAGATTAATGAATATCAACGAGGatat                                                                                                                      >  8:1321155/1‑150 (MQ=19)

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

                              gCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAAGAAAAAGAC‑A‑G‑ATCACGCGTGTGGTGTCTTTTGCACGGGCAGAAACATGAATGACTGAAAAAGAAGCTATAc                                                                                                                    >  3:332766/1‑151 (MQ=255)

                                     aGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGCTGTAGCTAAAAACG‑C‑T‑GATCTCCCGTCTGATGGCATTGGAGCGGCAGAAAGATGCTTGGCTGTCAACTACGGTATACCCata                                                                                                               >  1:1662347/1‑147 (MQ=255)

                                        gATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACGGTCCACGACGCTATACCCAaaagaaa                                                                                                          <  1:1574690/151‑1 (MQ=255)

                                           tGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCATTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGTTTTTGGAGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGc                                                                                                        <  2:546502/150‑1 (MQ=255)

                                             ggggCCTAGATGAATGACTGTCCACGACAGAACCCGGGTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCATAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGc                                                                                                     <  2:1831630/149‑1 (MQ=255)

                                                gCCTAGATGAATGCCGGCCCACGACAGAACCCGGCTTTTGGGTCAGTTTCACCTGATTTACGTAAAACCC‑C‑G‑GTGGGGGGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTa                                                                                                  <  7:1376774/151‑1 (MQ=17)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                                <  3:1001471/149‑1 (MQ=255)

                                                      attaatgaatgGCCACGACTGAACCCGGATTTTCGTTCTGTTTCACGTGTTTTACTTAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGTTTTTGGCGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                            <  1:960881/148‑1 (MQ=14)

                                                      atcactgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGatgaa                                                                                                                                         <  5:1654918/101‑1 (MQ=255)

                                                        taatgaCTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                           <  1:1073301/149‑1 (MQ=255)

                                                        taatgaCTGCCCCCCACATAACCCGGATTATGGGTCAGTTTCACCTGATTTACTTAAAAACC‑C‑G‑CTTCGGCGGTTTTTTGTTTTTGTAGGGGCCCAAATATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                           <  2:1623110/149‑1 (MQ=14)

                                                            gacttgccACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGAAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                       <  7:1279991/144‑1 (MQ=255)

                                                             aCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGATATTAAAAGCGAAGAAAAATGAAAGACTGTAAACAAAACTATACAAAAAAGAAATCGCCGTATcattaaaaaaa                                                                                      >  4:1384408/1‑140 (MQ=9)

                                                               tGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg                                                                                                                                       >  5:1890188/1‑99 (MQ=255)

                                                               tGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAACAAGAATGACTGTCCACTACGCTATACCGAATAGAAATAGTCTTATAGTTCATTTTCAcc                                                                                   >  4:1951334/1‑151 (MQ=14)

                                                                ggccACGACATAACCCGGCTTATGGGTCAGTTTCACCTGATTTACTAAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCCGAAACATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCt                                                                                  <  1:226156/149‑1 (MQ=14)

                                                                 tCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAAAC‑C‑G‑GTACCGCGGCTTTTAACTATTCTGACGTCACCAAGAAGAATGATTGTGCCGGACCCTATAGCCAAATGAAAGAGGCTTATCGGTCAGtgtggccgg                                                                                 >  7:29258/1‑143 (MQ=255)

                                                                 gccaccaCAGAACCCGGCTTATCGGTCAGTTGCACGTGATTTCCTTAAAACCC‑G‑G‑CTGCGGCGGGTTTTTGTTTTTGGCGGGGAATAATCATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTg                                                                                 <  6:873136/150‑1 (MQ=14)

                                                                     acacacAACCCGGCTGATCGGGAAGTTTCACCTGATTTAGGTAAAACCC‑G‑G‑CTTCGGCGGGTTTTTGGTTTTGCCGGGGCCAAAAGTTGAATGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                             <  1:195697/149‑1 (MQ=14)

                                                                      tactgAACCCGGCTTATCGGTCATTTGCTCCTGATTACCTTACAAACG‑C‑T‑CTTCTCGGGGTTTTTGCTATTGTTGGCGTCCTAAGTTGAATGACTCTCCTCTTCCCTATACCCAATAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                             <  2:346450/146‑1 (MQ=2)

                                                                      gACAGAACCCGGCTTATCGGTCAGTTTCACCTGAATAACGTAAAAACC‑A‑C‑CTTATGCGGAATTTTGGTTATATATGGGGAGAAATATGAATGACTCTCAACAAAAAAAAAAAAAAAAGAAAGCGGCTTATCGGTCAGTTTGCCCTggagtt                                                                            >  6:1381840/1‑147 (MQ=2)

                                                                      caCAGAACCCGGATTAGGGGGAAGTTGCCCCTCA‑TTACCGAAAAACC‑C‑GCCTGCGGCGGGTTTTTGGTTTTGCTGGCGCCCAACCATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                                            <  1:570544/150‑1 (MQ=2)

                                                                        cAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTGCGAAAAAACC‑C‑C‑CTGCTGGGCGGTTTTGCATAAGGACCGGCACTAAGATGAATGACTGTGCTCGAAAATAAACCCAAAAGAAAGCGGCGTAGCGGGCCGGTTCACGTGGTTTGCg                                                                          >  1:1746684/1‑151 (MQ=2)

                                                                          gaccgggggttAGCGGTCTGTTGCACCTGA‑TTACGTAAACACCGT‑C‑CTGCGCCGGGTTTTTGTTTTTGCCGGCGGCTCAGCATGATGGCCTCGACACGACGCTCTACCCAAAAGAAACCGGCTTATCGGTCAGTTTCACCTGGTTTACg                                                                          <  3:1138658/140‑1 (MQ=2)

                                                                           aaCCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑A‑CATCGGAAGATATTAGATTATGGAAAGACAGAAAAAATAATTAAAGAAAACGAAGGTATACCAAAAATAAAAAAAAATATCGGTCAGTTTCACCTGGTTTACGTa                                                                        >  4:588677/1‑150 (MQ=9)

                                                                            aCCCGGATTATGGGGAATTTTCCCGTGATTTCCTTCACACCG‑C‑T‑CTTCGGCGGGTTTTTGGTTTTGACGCCGCCCAAGTTTGATTGACTTGCCAAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaa                                                                      <  1:1366433/151‑1 (MQ=2)

                                                                             ccgggattATCGGAAATTTTCCCGTGATTTAGTTAAAACCG‑C‑T‑CTTCGGGGGGTTTTTGTTTTTGTAGGGGCATAAATTTGATTGACTGTACAAGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaa                                                                     <  4:1015969/145‑1 (MQ=2)

                                                                                    tATGGGTAATTTTCCCCTGATTT‑CGTTAACCCC‑GGG‑CTTCTGCGGGTTTTTGCTTTTGTTGGGGCCCAAATATGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGc                                                               <  6:1665039/150‑1 (MQ=2)

                                                                                     aTCGGTCAGTTGCCCCTGAATTAGTTAAAACCC‑C‑G‑GTTCGGCGGGTTTTTGCTTTTGTAGGGGCACAAATATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCt                                                              <  1:367713/150‑1 (MQ=11)

                                                                                       gggTAAGTTTCTCGTGATTTACTTAAAACCC‑C‑T‑CTGCGGCGGGTTTTTGCTTTTGTTGGGGCCTAACGATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                                           <  5:274117/150‑1 (MQ=11)

                                                                                       cgggaaGTTTCCCCTCATTTTC‑TAAACACCGC‑T‑CTTCTGCGGGTTTTTGTTTTTGATGGGGTATAATTATGAATGACTTTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTc                                                            <  6:1509722/145‑1 (MQ=2)

                                                                                       cGGTAAGTTTCCCCTGATTTACGTAAAACCC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGCGGGGCATAAAGATGAATGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTc                                                            <  5:41920/150‑1 (MQ=11)

                                                                                        ggTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTGTACGTAATAACCCGGTGcgg                                                          >  1:1859217/1‑151 (MQ=255)

                                                                                        ggTCAGTTTCACCTGATTT‑CGTTAAAACCGG‑G‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCAGAAATATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcgg                                                          <  6:573157/151‑1 (MQ=11)

                                                                                         gtaattttCACGTTATTTAGTTAAAAACC‑G‑T‑CTTCGGCGGTTTTTTGTTTTTGTAGGTGTTTAATTATGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggc                                                         <  4:642165/146‑1 (MQ=11)

                                                                                            aGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAAACCGCTTCGGCggg                                                      >  7:1112466/1‑151 (MQ=14)

                                                                                                ccACGTGATTTACTAAAAACCC‑C‑C‑CTGCGGCGGGTTTTTGGTTTTGCCGGGGCAGAATTATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttt                                                   <  3:849980/149‑1 (MQ=11)

                                                                                                 cACGTGATTTACTTAAAAACG‑G‑T‑CTGCGGCGGGTTTTTGTTTTTGTTGGGGCTTAAATTTGAATGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                                 <  2:1519422/151‑1 (MQ=11)

                                                                                                   ccTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTTTACGTAATAACCCGGTCCGGCGGGtttg                                                  >  3:1478440/1‑147 (MQ=255)

                                                                                                   ccTCATTTACGTAACCCCC‑C‑C‑CTGCGGCGGGTTTTTGCTTTTGCCGGGGCCCAACGATGAAGGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGc                                               <  5:1005184/151‑1 (MQ=11)

                                                                                                    cTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTATATCCTAATAGAGTCGAAGAAAGATAAATTAATTTCCAAAACGCTAAACACAAAAGAAAGCGTATAATATATCAAATTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttggtgcg                                              >  4:393540/1‑145 (MQ=2)

                                                                                                     tGATTTAGTTAAAAACG‑G‑G‑CTTCGGCGGGTTTTTGTTTTTGTTGGTGCATAATTATGATTGACTGGCCACGTCTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                                              <  8:1892752/150‑1 (MQ=11)

                                                                                                     tGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTGTACGTAAAAACCCGCTCCGGCGGGTTTCTGCtt                                             >  5:541862/1‑151 (MQ=255)

                                                                                                       aTTTCCGTAAAACCC‑C‑G‑CTGCGGCGGGTTTTTGCTTTGGCCGGGGCACAAACATCACTGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCGCCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                                           <  1:1452075/151‑1 (MQ=11)

                                                                                                       aTTTCCGAAAAACCC‑C‑C‑CTGCGCCGGGTTTTTGGTTTTGCCGGCGCACAACAATGAAGGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                                           <  5:1471165/151‑1 (MQ=11)

                                                                                                        tttcggTATAAACG‑G‑G‑CTGCGGCGGGTTTTTGCTTTTGCTGGTGCTTTCTTTTGATTGACTGTCGACTTCGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                                           <  2:1772656/145‑1 (MQ=2)

                                                                                                         aaccGTAAAACCC‑G‑G‑CTTCGGGGGGTTTTTGCTTTTGGAGGGCCACAAAGATCAATGCCTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTgg                                         <  7:636127/148‑1 (MQ=11)

                                                                                                           tccgtaaccccgc‑G‑CTGCGGCGGGTTTTTGGTTTTGCCGGCGCCCACCCATGATGGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTgg                                         <  3:1569448/138‑1 (MQ=11)

                                                                                                           aCGTAAAACCC‑C‑T‑CTTCTGCGGGTTTTTGCTTTTGGTGGTGTATAAATATGATTGACTGTCCACGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                                        <  8:683998/150‑1 (MQ=11)

                                                                                                           aCGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTGTTGGAGGGGCAGAAAGCTGAATGACTGTCCCCGCCGCTATACGCAAAAGATAGCGGCGTATCGGTCGGTTTCACCTGGTTTACGTAAAAACCCGCGTCGGCGGGGTGTTGCtggggggg                                       >  3:1292206/1‑144 (MQ=11)

                                                                                                               aaacaCC‑C‑G‑CTGCGGCGGGTTTTTGGTTTTGCCGCGGCCCAACCATGAATGACTCGCCACGACGCTATACCCAACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGc                                   <  3:1836096/151‑1 (MQ=11)

                                                                                                               aaaacCC‑C‑G‑CTGCGGCGGGTTTTTGGTTTTGGCGGGGCACAAACATGAATGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                                     <  8:710/149‑1 (MQ=11)

                                                                                                                  accc‑c‑G‑CTGCGGCGGGTTTTTGCTTTTGTAGGGCCACACACTTCAATGACGGACCACGACTCTATACCCAAAAGAAACCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                <  7:702211/149‑1 (MQ=11)

                                                                                                                       ccttctgatgcggGTTTTTGTTTTTGTTCGTGCACAACTATGATTGACTTGCCACGACGCTATACCCTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                              <  4:1130122/142‑1 (MQ=11)

                                                                                                                       c‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGATGTACGTAAAAACCCGATACGGCGGGTTTCTGCTTGTTGAg                                       >  3:49189/1‑140 (MQ=11)

                                                                                                                         ccctgcggcggGTTTTTGGTTTTGGCGGGGCCCAAACATGCAGGCCTCCCCCCGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa                            <  7:478395/146‑1 (MQ=11)

                                                                                                                          cctgcggcggGTTTTTGTTTTTGGAGGGGAAGAAAGTTGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat                           <  3:720247/147‑1 (MQ=255)

                                                                                                                             gcggcggGTTTTTGCTTTTGAGGGGGCTTAACTATGATTGACTGGCCACGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga                         <  5:221079/149‑1 (MQ=11)

                                                                                                                             gagccggGTTTTTGCTTTTGGCGCCGCACAACCATGAATGCCTCGCCACCACGCTATACCCAACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg                          <  5:1093821/145‑1 (MQ=11)

                                                                                                                              cggcggGTTTTTGTTTTTGTCGGTGCATAATTTTGATGGACTCTACACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                       <  8:1209308/151‑1 (MQ=11)

                                                                                                                                     gttttggttttggTGGGGCAGTATGTTGAGTGACTGGCGTCGTCGCTATACCCCACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTc                <  6:693678/150‑1 (MQ=11)

                                                                                                                                      ttttGTTTTTGCCGGGGCACAAACATGAAGGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTc                <  7:181320/150‑1 (MQ=11)

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACC‑C‑G‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACC  >  DH1/613447‑613741

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