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.4%	333.6 / 54.0	258	intergenic (+458/+565)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base G (62/91); new base C (1/54); total (74/184)
Fisher's exact test for biased strand distribution p-value = 3.09e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCC  >  DH1/613447‑613728

gCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAg                                                                                                                                     <  8:1585069/151‑1 (MQ=255)

 cccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGa                                                                                                                                    <  2:1449881/151‑1 (MQ=255)

    gTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAg                                                                                                                                 <  5:147989/151‑1 (MQ=255)

    gTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAg                                                                                                                                 <  1:148932/151‑1 (MQ=255)

    gTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAg                                                                                                                                 <  6:1679547/151‑1 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTCCTTTTGGAGG‑GGCAGaa                                                                                                                                   >  7:79158/1‑147 (MQ=255)

          aagccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGAGAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGatgaat                                                                                                                           <  6:641847/148‑1 (MQ=255)

           gtcttgctatctttctTGAGCGAGTGATCTATTGATGGCCTAGATTATTGACTCTACACCAAAGAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGG‑GGCAGAAAGatgaatg                                                                                                                          <  4:2031667/137‑1 (MQ=14)

           aCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGatgaatg                                                                                                                          <  1:1583551/151‑1 (MQ=255)

            cccGGGTAGGCTGCTTGAGCCAGTGAGTGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGatgaatga                                                                                                                         <  6:15949/151‑1 (MQ=255)

            cccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGatgaatga                                                                                                                         <  5:441237/151‑1 (MQ=255)

            cccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGatgaatga                                                                                                                         <  8:215390/151‑1 (MQ=255)

               gggTATGCTGGTTGCGCCAGGGAGCGATTGCTGGCCTAGATCACACACTCTCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGG‑GGCAGAAAGATGAATGACTg                                                                                                                      <  3:537776/151‑1 (MQ=17)

                  tattctttttgatgcAGGGATCGATTGATGGCCTAGATCACTGACTCTCCACCACAGAACGCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTTCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTc                                                                                                                    <  7:843126/136‑1 (MQ=17)

                  tAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACACGCTTCGGCGGCTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTcc                                                                                                                   >  5:106381/1‑151 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCTCGGGATTTACGTAAAAACACGCGGCGTCTGCTTTTTGCTTTTGGTGT‑GGCAGAAAGATGAATGACTGTccc                                                                                                                  >  5:18791/1‑150 (MQ=25)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTcc                                                                                                                   <  2:119436/150‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTcc                                                                                                                   <  7:321564/150‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCa                                                                                                                  <  1:1916189/151‑1 (MQ=255)

                     gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCacg                                                                                                                <  7:1437662/151‑1 (MQ=255)

                       tgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAAGTAAAAACCCGCTTCGGCGGGGTTTTGCTGTTGGAGG‑GGCAGAAAGATGAATGACTGTCCacgac                                                                                                              >  3:1228056/1‑151 (MQ=255)

                       tgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACATAAAAACCCACTTCGGAAGGTAATTCGTTTTGGAAT‑GAGAGAAACAAGAATGACTGTCCaacac                                                                                                              >  4:92778/1‑147 (MQ=255)

                       ggctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCacgac                                                                                                              <  7:1001859/150‑1 (MQ=255)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCacgac                                                                                                              <  4:1413988/150‑1 (MQ=255)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCacgac                                                                                                              <  1:1896815/150‑1 (MQ=255)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCacgac                                                                                                              <  1:415675/150‑1 (MQ=255)

                          ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATCACTGACCacgacccg                                                                                                           >  7:190629/1‑148 (MQ=255)

                          ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATAAATGAcaatacactactct                                                                                                           >  2:585273/1‑138 (MQ=255)

                            gAGCCAGTTTGGTATTCCTGGCCTATATTAATCACTGTTCTCCCCTTTACCGGGATTATCGGTTAGTTTAACCTTATTTTCGTAAAACCCCGCTTCGGCGGTTTTTTTTTTTTGTAGG‑GGCACATCGATTAATGACTGTCCACGACGCta                                                                                                          <  8:1248067/150‑1 (MQ=255)

                            gAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCtat                                                                                                         <  2:1610488/151‑1 (MQ=255)

                            gAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGATTATAGGTCAGTTTCACCTGATTTAAGTAAAAACCAGCTTCGTCGGCTTTTTGCTTTTGGAGT‑GGCAGAAAGATGAATGACTGTCCACGACGCtat                                                                                                         >  4:617026/1‑151 (MQ=21)

                             aGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCacgaggggatc                                                                                                        >  5:285502/1‑144 (MQ=255)

                               ccAGTGAGCGATTGGTGGCCTAGATGAATGACTGTCCACTACAGAACCCGGCTTATCGGTCATTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGTAGG‑GGCTGAAAGATTAATTACTGTCCACGACGCTATAcc                                                                                                      <  8:1698482/151‑1 (MQ=21)

                                 aGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTATCACCTCATTTTCGTAAAAAACAGCTTAGGCGGGTTTTTCCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCGACGCCGCTATAcccc                                                                                                    >  1:1372658/1‑150 (MQ=21)

                                  gtgaacctttGCTGGCCTAGGCTACACACTATCCCCAACAGAACCCGGCTTAGCGGGCAGTTTCACGTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGTAGG‑GGCAGAAAGATGAATGACTGTCCACGACACTATACCCaa                                                                                                   <  1:751922/143‑1 (MQ=11)

                                    gAGCGATGGCGGGCCTACAGGAATGACTGGCCACGCCAGAACCCGGGTTAGCGGGCCGTTCCACCTGTTTTACGTAAAACCCCGCTGCGGGGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCaaaa                                                                                                 <  3:722171/151‑1 (MQ=17)

                                       cGATTGCTGGCCTAGATGAATGCCTGTCCACTACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTGCGGCGGGTTTTTGCTTTTGGCGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                              <  5:2065317/151‑1 (MQ=255)

                                         aTTGCTGGCCTATTTTAATGACTGTCCACGACAGAACCCGGTTTATGGGTCAGTTTCACCTTATTTAGTTAAAAACCCGCTTCGGGTGTTTTTTGTTTTTGGAGG‑GGCATAAATATGAATGACTGTCCACGACGCTATACCCAaaagaaa                                                                                             <  8:865063/150‑1 (MQ=255)

                                         aTTGCTGGCCTAGATGAATGACTGCCCACTACAGAACCCGGCTTTTGGGTCAGTTTCACGTTATTTACGTAAAACCCCGCTTCGGGGTGTTTTTGTTTTTGGAGG‑GGCAGAAAGATGAATCACTGCCCACGACGCTATACCCAaaagaaag                                                                                            <  5:764671/151‑1 (MQ=18)

                                            gCTGGCCTAGATGAATGACTGTCCACGACAGAAAACCCATAATCAGTCAGATTAAACTGAATAACGTAAAAAACAGCTTCGGAGTGAATATAATTTTAGAGC‑TCCATAAAGATGAATGAATTTACAAAAAACTATACCCAAAAGAAAGCgg                                                                                         >  4:1167690/1‑151 (MQ=255)

                                               ggCCTAGATGAATGACTGTCCCCCACAGAACCCGGCGTATCGGTAAGTTTCACCTGATTTACTTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAACGAAAGCGGCt                                                                                       <  2:705235/150‑1 (MQ=17)

                                                cccTATATGAATGACTGCCCACGACAGACCCCGGATTACCGGGCAGTTGCACCTGATTTACGTAAAACCCCGCTGCGGCGGGTTTTTGCTTTTGGCGG‑GGCACAAAGATGAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTa                                                                                     <  1:2002855/150‑1 (MQ=17)

                                                 cctcgaaaacacccTGTCCCCAACATAACGCGGCTTATCGGTCATTTTCACGTGTTTTACTTAAAAACCCGCTTCGGCGGTTTTTTGGTTTTTGAGG‑GGCACAAAGATTAATGACTGTCCACGACACTATACCCAAAAGAAAGCGGCTTAt                                                                                    <  1:1201363/138‑1 (MQ=1)

                                                 cattgATTAATGACTGTCCACGACATAACCCGGTTTATGGGTCAGTTGCACGTGATTTACTTAAAAACCCGCTTCGGGGGTTTTTTGGTTTTGTAGG‑GGCATAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTAt                                                                                    <  1:1415232/147‑1 (MQ=17)

                                                  atAGATGAATCACTGTCCACTACAGAACCCGGCTTATCGGTAAGTTTCACCTTATTTACGTAAAAACCCGCTTCGGGGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                   <  7:1220424/150‑1 (MQ=17)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTGTACCCAAAAGAAAGCgggatcgggg                                                                                 >  5:1672759/1‑143 (MQ=255)

                                                    aGATGAATGACTGGCCACGACAGAACCGGGCTTTTCTGTCAGTTTCCCCTGATTTACGTAAAACCGCGGTTCGGCGGGTTTTTGCTTTTGGTGG‑GGCAGAATGATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCg                                                                                  <  1:218534/150‑1 (MQ=17)

                                                    aGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGGGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCgg                                                                                 <  8:1171293/151‑1 (MQ=255)

                                                       tgaatgaCTGTCCACGACAGAACCCGCATTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAGATCGGCTGAGCgggc                                                                               >  5:899954/1‑148 (MQ=17)

                                                         aatgaCTGGCCACGACAGACCCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGTAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                             <  2:408950/150‑1 (MQ=255)

                                                          atgaCTGGCCACGACAGAACCCGGCTTATGGGGCAGTTTCACCTGATTTACGTAAAAACCCGCTGCGGCGGTTTTTTGGTTTTGGAGG‑GGCATAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGtt                                                                           <  3:917626/151‑1 (MQ=17)

                                                            gacgggccACGACAGAACCCGGCTTATCGGTCAGTTGCACCTGATTTACGTAAAAACCGGCTGCGGCGGGTTTTTGCTTTTGTAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                          <  1:1392550/144‑1 (MQ=14)

                                                            cccTGCCCACCACATAACCCGGCTTATCGGTCTGTTTCACCTGATTAAGGTAAAACCCCGCTGCGGCGGGTTTTTGCTTTTGTCGG‑GGCACAACTATCAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                          <  6:787732/148‑1 (MQ=14)

                                                             aCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTGACCTGATATAACTAAAAACCCTCTTCTACGGGATTATGCTTTAGCAGA‑GGCAAAAAGATGAATGACTGACCAAGAATATATAAAAAAAAGAAAGCGGCTTATCGGTCAGtttt                                                                         >  4:346502/1‑149 (MQ=9)

                                                             aCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACATAAAAACACCCTTCGGAAGGTATTTACTTTTTAAAT‑TACAGAAAAAATAATTACTCTCCAAGAAAATATAAAAAAAAGAAAGAGGCTTATCGGTCAGtttg                                                                         >  8:541873/1‑149 (MQ=255)

                                                                ggccaccaCACAACCCGGCTTATCGTTCAGTTGCACCTGATTACCGAAAAACCCCCCTGCGGCGGGTTTTTGTTTTTGTCGG‑GGCACAAACTTGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCt                                                                     <  3:356490/149‑1 (MQ=14)

                                                                gTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGAAAAACCAAGCGTCGGCGTCTGATGGCGATTCGAGG‑TGCTGAAAGATGAATGACTGTCAACGAGTCTATAACCAAAAGAAAGCGGCTGATCGGTCGGtttggcgg                                                                     >  7:449741/1‑146 (MQ=9)

                                                                 gcgaCGACTGAACCCGGATTAGGGGGCAGTTTCACCTGATTTACGTAAAAACCGGCTTCGGCGGGTTTTTGCTTTTGCTGG‑CGCACAAATTTGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTg                                                                    <  5:859069/148‑1 (MQ=14)

                                                                 gccACGACAGAACCCGGATTATCGGTCAGTTTCACCTGATTTACTTAAAAACCCTCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAATATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTg                                                                    <  5:1975198/150‑1 (MQ=255)

                                                                  ccACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAACCCGCTACGGCGGGTTTATGCTTTTTGAGT‑GCCATAAACATAAATGACTCTACACGCCACTAAAACCAAAAAAAAGATTATTATAGTTCAGTTTAACCtt                                                                    >  8:1879078/1‑149 (MQ=255)

                                                                   cACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACCCTATACCCaaa                                                                                                  >  3:618450/1‑119 (MQ=255)

                                                                      caCAGAACCCGGCTTATCGGGAATTTTCACGTGATTTACTTAAACACCCGCTGCGGCGGGTTTTTGCTTTTGGAGG‑GGCATAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                               <  1:1126254/150‑1 (MQ=14)

                                                                      aactgAACCCGGATTATGGGCAATTTTCAGCTGATTTCCGTACACACCGTCTGCGGCGGGTTTTTGCGTTTGTCGG‑CGCCCAACCTTGATGGACTCTACACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                               <  4:792900/147‑1 (MQ=2)

                                                                       ccAGAACCCGGCTTAGCGGGCAGTTGCACCTGATTTACGTAAAAACCCGCTGCGGCGGGTTTTTGCTTTTGGCGG‑GGCCGAACGATGAATGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTAc                                                              <  3:253761/150‑1 (MQ=14)

                                                                           aaCCCGGCTTATCGGGCAGTTTCCCCTGATTTCCGAAAACCCCCGCTGCGGCGGGTTTTTGGTTTTGTCGG‑CGCACAACGATGAATGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaa                                                          <  1:590017/151‑1 (MQ=14)

                                                                            aCCCTGCTTATCGGTGATTTTCAGCTTTTTTAATTAAAATGGCGCTTCGCCGTGTTTTTGGTTTTGGTGG‑GGGGCAACCTTGAATGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaa                                                          <  5:1332361/150‑1 (MQ=2)

                                                                             cccGGCTTATCGGGCATTTGCCCGTGATTTACGTAACACCCCGCTGCGGCGGGTTTTTGCTTTTGCAGG‑GGCACAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaa                                                        <  5:335017/151‑1 (MQ=11)

                                                                             cccGGCTTATCGGGCAGTTTCACCTGATTTACGTAAAACCCCGCTGAGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAAGGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaa                                                        <  1:85638/151‑1 (MQ=255)

                                                                               cGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTTCTTTTCGAGG‑GTCAGAAACATTATTAACTCAAGACGAATGGATAAACAAAAGAAGGCGGCTTATGGGTCAGTAAAACCTGGTTTACGTaaactc                                                      >  3:666835/1‑148 (MQ=2)

                                                                                   ttATCGGTCAGTTTCACCTGATTTACGTAAAAAAACGCTTCGGCGGGGTTTTGCTTTTGGAGG‑GGAAGAAAGATGAATGACTGTCCACGACGCTATCCAAAAAAGAAAGCGGCTTATCGGTCAGTGTGACCTGGTGTACGTaatagctcgg                                                  >  5:1190879/1‑143 (MQ=11)

                                                                                    aatggttaaGTTTCCCGTGATTTACGTAAAACCCGTCTTAGGCGGGTTTTTGCTTTTGTTGG‑TGCTAAATTTTGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGc                                                  <  4:1637051/142‑1 (MQ=11)

                                                                                     atgggTAAGTTTCCGGTGATTTACTTAAAAACCGGCTTCGGCGGGTTTTTGTTTTTGTTGG‑GGCACAAATATGAATGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCtt                                                <  8:1944580/148‑1 (MQ=11)

                                                                                     aTCGGTCAGTTTCACCTGATTTACGTAAAATCCCGCTTCGGAGCTCTTTTGCTGATCGCTC‑TGAAGAAAAATCACATACAGAACACGCCGATATACCCAAAAGAAAGCGGATTATCGGTCAGGTTCCCCTGGGTTTCGTaaaaggcggctg                                                >  7:1813504/1‑143 (MQ=255)

                                                                                           aatttgcACGTGATTTAGGTAAAAACCTGCTTCTGCGGGTTTTTGTTTTTGCTGG‑TGACCAAATATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggcgg                                          <  1:509307/145‑1 (MQ=11)

                                                                                            attttctggtGATTTTCTAAAAACCCCTGTTCGGCGGTTTTTTGTTTTTGAGGG‑TGCTAATTATTGATTGACTTGACTCGTCTCTATACCCAAAAGAACGCGGCTTATCGGTCAGTTTCACCTTGTTTACGTAAAAACCCGCTTCGGCggg                                         <  6:1822053/142‑1 (MQ=2)

                                                                                            attttCACGTGATTTACTAAATACCCCTCTTCGGGGGTTTTTTGCTTTTGTTGG‑CGTCCAAAGATGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCggg                                         <  7:1655389/149‑1 (MQ=11)

                                                                                              tttccgcTGATTTCCTTATAAACGTGCTGCTGCGGGTTTTTGCTTTTGTTGG‑CGCTTAATTATGATTGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGt                                        <  2:1814276/144‑1 (MQ=2)

                                                                                              tttCACCTGATTTACGTAAAAACCCGATTAGAACGGATTATCCTATTGGAGG‑GGCAGAGAGATAAAATAATTTCAACTAACCTATAACCAAAACAAAGATACTAATCAGTCAGTTTGACCTGATGTACGTAAATAACCggtcagaagggat                                       >  6:244406/1‑139 (MQ=255)

                                                                                                tCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCGTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGCTGTCCGTAATAACGCTCGGCGGCggggagg                                     >  5:1301857/1‑147 (MQ=11)

                                                                                                 gccatGATTTCCGTAACAACCGCCTTCGGCGGTTTTTTGTTTTTGTCGG‑GGCTTAATTATGATTGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtttt                                     <  6:1268240/146‑1 (MQ=11)

                                                                                                 cACCTGATTTACTTATAAACGCTCTTCGGCGGGTTTTTGTTTTTGCTGG‑CGCTTAATTTTGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                    <  4:472559/151‑1 (MQ=11)

                                                                                                 cACCTGATTACCGTAAAACCCCCCTGCGGCGGGTTTTTGGTTTTGTTGC‑TGCAAAATAATGACGTACTCGCGACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                    <  7:947833/151‑1 (MQ=11)

                                                                                                    cTGATTTACGTAAAAACCCGCTTCGGCGGCTATGTGATGTTCGAGG‑GGCTGAAAGAGGAATGCCTGGACACGACGCTTTACCCATAAGAACGCTGCTAATAGCTCAGTGGGACCGGGTGGACGTAAACTGCCGCGTGGGCgggggtggggg                                 >  7:517317/1‑143 (MQ=2)

                                                                                                      gATTTACGTAAAACCCCGCTGCGGCGGGTTTTTGCTTTTGTCCG‑GGCACAACCATGAATGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                                 <  3:1154182/149‑1 (MQ=11)

                                                                                                      caTTTACGTAAAACCCCGCTGCGGCGGGTTTTTGGTTTTGGAGG‑GGCCCAACGATCAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                                 <  3:1955677/148‑1 (MQ=11)

                                                                                                       atatACGTAAACCCCCGCTTCGCCGGGTTTTTGCTTTTGGAGG‑CGCAGAACGATGCATCCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                              <  3:1857072/148‑1 (MQ=11)

                                                                                                       aTTTCCGTAACAACCGGCTGCGGCGGGTTTTTGTTTTTGTTGG‑GGGAGAATGTTGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCttt                               <  8:479368/150‑1 (MQ=11)

                                                                                                       aTTTCCGTAAAACCCCGCTTCGGCGGGTTTTTGATTTTGCTGG‑TGCACAATTATGATTGACTTTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                              <  2:936014/151‑1 (MQ=11)

                                                                                                        tttcgttAAAACCCGGCTGCGGCGGGTTTTTGGTTTTGTTGG‑GGCTACATCATGATTGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                              <  2:2041417/144‑1 (MQ=11)

                                                                                                         ttACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTGCCCGACGGGAGACCCATAAGAACGAGGCTGATCGGTCAGTTTCCCCTGGGTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTGTgg                            >  7:1098001/1‑151 (MQ=11)

                                                                                                           ccGTAAACCCCGGCTTCGGCGGGTTTTTGGTTTTGGCGG‑GGCACAACCATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAg                          <  1:991977/150‑1 (MQ=11)

                                                                                                              tAAAAACCCTGTTCGGGGGGTTTTTGCTTTTGGTGG‑TGCAGAATGATGATTGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                        <  6:1922873/150‑1 (MQ=11)

                                                                                                                cagccggtctgctgcggGTTTTTGGTTTTGTCGG‑GGCCTCATTATGATGGCCTTGCCTCGACTCTATCCACTAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCa                     <  7:1987871/138‑1 (MQ=2)

                                                                                                                    ccgtcTTCGGCGGGTTTTTGGTTTTGGCCG‑GGCACCCAGATGAATGACTCGCCACGACGCTATACCCTACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                 <  1:28867/147‑1 (MQ=11)

                                                                                                                     gcacTTCTGCTGGTTTTTGTTTTTGTAGG‑CGTATTATCTTGATTGACTCTCGTCGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAg                <  2:392322/148‑1 (MQ=11)

                                                                                                                      gtcTTCGGGGGTTTTTTGCTTTTGGAGG‑GGCTTTCTTATGATTGAGTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa               <  4:236426/149‑1 (MQ=11)

                                                                                                                       gctgccgcggGTTATTGTTTTTGTTGG‑TGCTCAATTATGAATGACTTGACTAGACTCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat              <  2:51513/145‑1 (MQ=11)

                                                                                                                          gcggcggGTTTTTGTTTTTGGTGG‑AGGCTAATCATGAATGACTCGACACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg             <  3:1507329/148‑1 (MQ=11)

                                                                                                                          gatgcggGTTTTTGCTTTTGGGGGTGTCATAAA‑TTGAAGGCCTAGCCAAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg             <  5:943261/146‑1 (MQ=11)

                                                                                                                             gcggGTTTTTGGTTTTGATGC‑GGCGACAATTTGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatga        <  1:33992/151‑1 (MQ=11)

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGG‑GGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCC  >  DH1/613447‑613728

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: `‑`