BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	613,597	0	G	C	23.9%	224.4 / 26.1	152	intergenic (+458/+565)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base G (56/49); new base C (2/31); total (60/92)
Fisher's exact test for biased strand distribution p-value = 4.56e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

CCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTAT  >  DH1/613448‑613738

cccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                                                                                                                              <  2:156825/151‑1 (MQ=255)

  cGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTACATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                             >  2:390270/1‑150 (MQ=255)

   gTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAATATCGGG‑CTTCGGCGGGTGTTTGCTGTTGGAGGGGCTGaa                                                                                                                                             >  1:866021/1‑149 (MQ=37)

     cctcaccccGGGTAGGCTGCTTCAGCCAGGGAGCGTTGGCTGGCCTATATGAACGACTGGCCACCACAGACCCCGGCTTAGCGGTCAGTTTCCCCTGATTTACTAAAAAACCCG‑CTTCGGCGGTTTTTTGCTTTTGGAGGGGCAGaaaca                                                                                                                                          <  8:941862/144‑3 (MQ=18)

     aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat                                                                                                                                         <  3:809396/151‑1 (MQ=255)

     aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa                                                                                                                                          <  7:521795/150‑1 (MQ=255)

     aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa                                                                                                                                          <  6:595849/150‑1 (MQ=255)

      cTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg                                                                                                                                        >  1:383577/1‑151 (MQ=255)

       tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGTAGGGGCAGAAAGatg                                                                                                                                        >  2:144487/1‑150 (MQ=255)

         aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                                                                                                                                     <  1:438495/151‑1 (MQ=255)

           cccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg                                                                                                                                    >  3:44864/1‑150 (MQ=255)

               ggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGt                                                                                                                               >  7:946354/1‑151 (MQ=255)

               ggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGt                                                                                                                               <  1:575409/151‑1 (MQ=255)

                gTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTc                                                                                                                              <  5:873772/151‑1 (MQ=255)

                 tAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCGACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATAACtctc                                                                                                                              >  5:256486/1‑150 (MQ=255)

                  aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTcc                                                                                                                             <  3:825693/150‑1 (MQ=255)

                  aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAACACGACTTATAGGTCAGTTTAAACTGATTTACGTAAAAACCAG‑ATTAGGAGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTcgc                                                                                                                            >  3:329620/1‑149 (MQ=21)

                      tgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGATGGTTATTGCTGTGGGAGGCGCAGAACGCTGAAGGACTGTCCccgac                                                                                                                        >  7:26688/1‑151 (MQ=25)

                      ggctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                        <  6:452259/150‑1 (MQ=255)

                       gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGGCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                                       <  5:266125/151‑1 (MQ=255)

                         ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACGGAAGCCGGCTTATCGGTCAGTTTCACCTTATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAATGATGAATGACTGTCGacgtcggt                                                                                                                     >  3:2786/1‑149 (MQ=21)

                         ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACGTGATGTAGCTAAAACCGCT‑CATCGGCGGTGTGTTCCGATTGGAGGGGCAGCACGATGAATGACTGTCCACGACGCt                                                                                                                     >  3:899024/1‑151 (MQ=25)

                          tGAGCCTGTGCCCTATTGCTGGCCTAGATGAATGACTGTCCACGAAATAACCCGGCTTATCGGTCATTTTCACGTGATTTACTTAAAAACCCG‑CTGGGGGGGGTTTTTGCTTTTGGTGGGGCACAAAGATGAATGACTGTCCACGACGCt                                                                                                                     <  2:294937/150‑1 (MQ=18)

                           gcgCCTGTGAGCTATGGCTGGCATAGATGAATCACTGTCCACGACAGAACCCGGCTTATCGGTCATTTTCACCTTATTTACTTAAAAACCCG‑GTTCGGCGGTTTTTTGTTTTTGGAGGGGCATAAATATTAATTACTGTCCACGACGCtat                                                                                                                   <  8:494084/149‑1 (MQ=255)

                            cgCCAGTGAGCGATTGCTGGCCTAGATTAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCtata                                                                                                                  <  4:101468/150‑1 (MQ=255)

                            aGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCAG‑CTTATGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTAAacgaacatat                                                                                                                   >  6:398660/1‑144 (MQ=21)

                             gCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACATAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATAc                                                                                                                 <  4:663725/151‑1 (MQ=255)

                               cAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCT‑CTTCGGCGGGTTTTTGCGTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACCCTATAccc                                                                                                               >  5:801809/1‑151 (MQ=255)

                               cAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAAACCG‑ATACGGCGGGTTTAAGATTTTTGAGTGGCATATATATAAAATACTATGCACGACGTTATAcg                                                                                                                >  4:236103/1‑149 (MQ=21)

                               cAGTGAGCGATTGCTGGCCTAGAGGAATGACTGTCCCCGACAGAACCCGGCTTATCGGTCAGTTTCACCTTATTTACGTAAAAACCCG‑CTGCGGGGGTTTTTTGCTTTTGGAGGGGCAGAAAGATCAATGACTGGCCACGACGCTATAccc                                                                                                               <  7:821967/151‑1 (MQ=21)

                                aGTGAGCGTTTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCa                                                                                                              <  1:147141/151‑1 (MQ=255)

                                aGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGCTCAGATTAACCTTAATAAAGTAAAAAACCG‑ATTATGAGGGATATTGATTTTTGAGCGGAAGAAAGATGAATGAAAGTACAATACACTATAcac                                                                                                               >  6:464591/1‑148 (MQ=255)

                                 gTGCGCGATTGCGGGCCTAGTTGAATTACTGTCCACGACCGAACCGGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCATAAACATGAATGACTGTCCACGACGCTATAccc                                                                                                               <  2:510554/149‑1 (MQ=17)

                                    aGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaag                                                                                                          <  4:814951/151‑1 (MQ=255)

                                     gctaTTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCGGGGTTATCGGTCAGTTTCACCTGATTTACGTAAAACCCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGGCCACGACGCTATACCCAaaaga                                                                                                         <  1:666675/148‑1 (MQ=18)

                                      cGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAGAACCGCTAAGCTGTCAGGTTCACGTGCTGTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGATGGGCAGAAAGATGAATGACTGTCGACGgcgcggtcgcgcacaggg                                                                                                        >  5:628162/1‑137 (MQ=17)

                                      cGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATGGGTCAGTTTCACCTGATTTACTTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                                        <  8:873196/151‑1 (MQ=255)

                                      cGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATAAATCTCAACAACGAAATACCCaaaaaa                                                                                                         >  2:440330/1‑148 (MQ=21)

                                       gtttGCTGGCCTATATGAATGACTGCCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTTATTTACGAAAAACCCCG‑CTTCGGGGGTTTTTTGCTTTTGGCGGGGCAGAAAGATTAATGACTGCCCACGACGCTATACCCAaaagaaa                                                                                                       <  7:165784/149‑1 (MQ=17)

                                       gtttGCTGGCCTAGATGAATGACTGCCCACGACAGAACCCGGCTTAGCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATTAATGACTGGCCACGACGCTATACCCAaaagaaa                                                                                                       <  5:1071989/149‑1 (MQ=255)

                                       gATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCCGTCAGGTTCACCTGGATTACGTAAAAACCCG‑CTTCGGCGCCTTTTGGGGTTTGGAGGGGCAGATAGATGAATGCCTGTCCACGACGCTATACCCAaaagca                                                                                                        >  5:287458/1‑148 (MQ=18)

                                       gATTGCGGGCCTAGATGAATGACTGGCCACCCCAGAACCGGGCTTATCGGGCTTTTTCACCTGATTTACGTAAAAACCGG‑GTTGGGCGGTTTTTTGGTTTTGGAGGGGCCGAAAGATGAATGACTGGCCACGACGCTATACCCAaaagaaa                                                                                                       <  7:191609/151‑1 (MQ=255)

                                        aTTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGGGGGTTTTTGCTTTTGGAGGGGCACAAAGATCAATGACTGTCCACGACGCTATACCCAaaagaaag                                                                                                      <  1:134753/151‑1 (MQ=255)

                                         ttGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGTGAAAGAAAGATCAAAGACTATCCAATACGCTATAAACAaaagaa                                                                                                        >  8:806526/1‑148 (MQ=18)

                                         ttGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGCTCAGTTTCACCAGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGAAACTATACCCAAAAGAAAGc                                                                                                     >  4:361316/1‑151 (MQ=255)

                                          ggCTGGCCTAGTTTAATGAGTGTCCACGACAGATCCCTGATTCTCGGCCAGTGTGACGTTATTTATTTAAAATTCCT‑GTTCGGGGGTTTTTTGCTTTTTTAGGGGGCGAATTATTATTGACTGTCCACGACGCTATACCCAAAAGAAAGc                                                                                                     <  4:1043786/149‑1 (MQ=255)

                                           gCTGGCCTAGATGAATGACTGCCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAAACCCG‑CTTCGGCGGTTTTTTGCTTTTGGAGGGGCAGAAAGATGAATCACTGTCCACGACGCTATACCCAAAAGAAAGCgg                                                                                                   <  7:274123/151‑1 (MQ=255)

                                              ggCCTAGATGAATGACTTTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTATTTAATTTTATG‑CTTCGGGGGTTTTTTGATTTTTTATGTGCAGAATTATGTATTACTGTACACGACGCTATACCCAAAAGAAAGCGGCt                                                                                                 <  2:910102/150‑1 (MQ=16)

                                              ggCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAATAAACCA‑CTTCGGCCGCTGTATTCTTATCGAGGGGCAGAAAGAAGAATCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCtt                                                                                                >  7:615133/1‑151 (MQ=18)

                                              ggCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGTAGGGGCAGAAATATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCtt                                                                                                <  6:307206/151‑1 (MQ=255)

                                               gCCTAGATGAATGACTGTCCACGACACAATACGCAAAATATCGCAGGATCACCGGAATTAACTAAAGACACT‑CGTCAGAGGCTTGTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCGCCGACGCTCTACCCacacgatcggcggagc                                                                                               >  3:607884/1‑136 (MQ=9)

                                                ccTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGATCAGTTTCACCAGATATAAATAAAAACCAG‑CTTCGAAAGCTTTTAGATTTTGGAGGTGCACAAAGATGAATGACTGTCAAAGAAAATATACCCAAAAGAAAGCGGCTTAt                                                                                              >  6:878996/1‑151 (MQ=17)

                                                  tAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAAACCG‑CGTCGCCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTACACGACGCTATACCCAAAAGAAAGCGGcgttggg                                                                                            >  3:297938/1‑145 (MQ=17)

                                                   aGATGACTGACTGTCCACGACAGAACCGGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                             <  4:608290/149‑1 (MQ=255)

                                                   aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGGGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                             <  1:233605/149‑1 (MQ=255)

                                                       gaatgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACACCGCTATACACAAAAGAAAACGGCTAATCTGTCat                                                                                       >  8:661165/1‑150 (MQ=17)

                                                       gaatgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAAAGAAAGTAAAAGCCTCTATAACCGaaagaaagacttataaagtta                                                                                         >  2:80202/1‑136 (MQ=17)

                                                       gaaggaggggccACGACAGAACCCGGCTTAGGGGGCAGTTTCACCTGATTTACGTAAAACCCCG‑CTTCGGCGGGTTTTTGGTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                        <  3:753307/140‑1 (MQ=14)

                                                       aacagtctttCCACCACATAACCCTTCTTATCTGTAATTTTCACGTGTTTTACGTAAAAACGGT‑CTTCTGCGGGTTTTTGGTTTTGTTGGGGCTCAAATATGAATGACTGTCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                       <  2:303251/142‑1 (MQ=2)

                                                        acagcctcaccACAACACAACCCGGCTTAGGGGGCATTTTCCCGTCTTTTACGTAAAACCCGG‑CTGCGGCGGGTTTTTGCTTTTGTCGGCGCCCAACGATGATTGACTGGCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGt                                                                                      <  1:345889/142‑1 (MQ=2)

                                                         atgaCGGCCCACGACAGAACCCGGCTTATCGGTCAGTTGCACCTCATTTCCGAAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGtt                                                                                     <  1:588155/151‑1 (MQ=14)

                                                           gaCTGTCCACGACAGAACCCGGCTTATCGGTCATTTTCAGCTGATTTACGTAAAACCCCG‑CTTCGGCGGGTTTTTGGTTTTGTAGGGGCTTAAAGATGAATGAGTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                    <  6:440425/150‑1 (MQ=17)

                                                            agtggccACGACAGAACCCGGCTTATCGGTCAGTTTCACGTTATTTACTTAAAACCCCG‑CTGCGGGGGGTTTTTGCTTTTGTTGGGGCAAAACTATGATTGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTc                                                                                   <  5:399298/145‑1 (MQ=14)

                                                            aCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAATAAGGCGGATGATCCGTCTGttgtc                                                                                  >  7:690331/1‑148 (MQ=17)

                                                                     caCAGAACCCGGCTTATGGGTCAGTTTCACGTGATTTCCTTAAAACCCCG‑CTGCGGCGGGTTTTTGGTTTTGCCGGGGCCCAAACATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                                         <  1:118467/150‑1 (MQ=14)

                                                                      ccAGAACCCGGATTATGGGTAATTTGCCCGTCATTTACGTAAAACCGCG‑CTTCGGGGGGTTTTTGCTTTTGGCGGGGCTCAAAGATGAATGACGGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTAc                                                                        <  1:905272/150‑1 (MQ=14)

                                                                       cAGAACCCGGCTTTTGGGTCAGTTTCACCTGATTTACGTAAAACCCGT‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACg                                                                       <  4:811589/151‑1 (MQ=14)

                                                                       cAGAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacg                                                                                                                          <  1:914929/100‑1 (MQ=2)

                                                                       cAGAACCCGGCTTATCGGTCAGTTTCACCTCATTTACGTAAAAACCCG‑CTGCGGCGGGTTTTTGCTTTTGGCGGGGCACAAAGATCAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACg                                                                       <  1:803659/151‑1 (MQ=14)

                                                                         caccccGGTTTCTGGGCCAGTTGCCCCTGATTTCGGTAATACGGGG‑CTGCGGGGGGTTTTTGTTTTTGGTGGCGCCTAAACATGATTGAGTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTa                                                                     <  2:1035051/148‑1 (MQ=2)

                                                                          aaCCCGGATTATCGGTCAGTTTCCCGTGATTTAGGTAAACCCCCC‑CTTCGGCGGGTTTTTGTTTTTGCTGGTGCTCAATTATGAATGCCTGTCCACGACCCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaa                                                                    <  3:1018200/151‑1 (MQ=11)

                                                                            cccGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGAGCGGATCCGCAAAAGAAAGCGCGACATCCGGAAGGGACACCTGGTTTACGTaaa                                                                   >  7:737689/1‑150 (MQ=11)

                                                                            cccGGATTATCGGTAATTTTCACCTGATTTCCTTAAAACCCGT‑CTTCTGCGGGTTTTTGCTTTTGTTGGCGCTCAATTATTATTGACTTTCAACGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaa                                                                  <  2:54791/151‑1 (MQ=2)

                                                                              cGGATTATGGGTCATTTTCACCTGATTTACTTAAAACCCCG‑CTTCCGCGGGTTTTTGCTTTTGCCGGTGCCCAAAGATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAc                                                                <  2:443470/151‑1 (MQ=11)

                                                                                    agggttaatttgcACGTGATTTTCTTAAACACCGT‑AGTCGGGGGGTTTTTGTTTTTTTTGGGGGAGAAGTATGATTGACTTGCGTCGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCt                                                           <  4:362915/138‑1 (MQ=2)

                                                                                     gcGGTCAGTTGCACCTGATTTACGTAAACACCCG‑CTTCGGCGGGTTTTTGCTTTTGCCGGGGCACAATTTTGAATGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTc                                                         <  7:862905/150‑1 (MQ=11)

                                                                                      ggccaattttCACCTGATTAACTTAAAAACGCT‑CTTCGGCGGGTTTTTGCTTTTGAAGCTGCTTAACCATGAATCACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                                        <  6:493806/144‑1 (MQ=11)

                                                                                        gccacttgcacgttaTTTCCGTAAAACCCCC‑CTTCGGCGGTTTTTTGCTTTTGGAGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggc                                                      <  7:148100/137‑1 (MQ=11)

                                                                                          cAGTTTCACCTGATTTACGTAAAAACCCG‑CTTAGAACGGTATTTGCTATTGGAGGGGAAGAAAGATTAATGAATTTCAAATACCCTATAACCAAAAAAAATATAAATAAAGGTCAGTTTAACCTGGTTTACGTAAAAACCCGGTTcggtgg                                                    >  2:319554/1‑151 (MQ=255)

                                                                                             tttCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTGTACGTAATAACCCGGTGCGGCggggt                                                 >  5:1006871/1‑149 (MQ=11)

                                                                                              tacACCTGATTAACGTAAAACCCCC‑CTTCGGCGGGTTTTTGCTTTTGCTGGGGCACAACCCTGAATGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGt                                                  <  3:541689/147‑1 (MQ=11)

                                                                                               gcgggtcatttcgttAAAACCCGG‑CTTCGGCGGGTTTTTGCTTTTGTCGCCGCCCACGCATGAAGGCCTCTCCAAGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtttt                                               <  1:884324/137‑1 (MQ=2)

                                                                                                   cTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGGAAGATGAATGACTGTCGACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGGGGGGGTTTggg                                           >  5:769027/1‑149 (MQ=11)

                                                                                                   cTGATTTACGTAAAAACCCG‑CTTAGAGCGGAATATCCTATTGTAGGGGAAGAAAGATTAATGAATGTAAACTAACCTATAACCAAAAAAAATATAAAAATAGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTAGGGGGGttgagggg                                           >  6:1005001/1‑145 (MQ=255)

                                                                                                    tGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAACGCGGCTGATCGGTCAGTTTCGCCTGGTTTACGTAAAAACCCGCTTCGGGGGGTTGTTGCtt                                          >  7:383175/1‑151 (MQ=255)

                                                                                                    tGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAAATACTAACAACGACGCTAAACCCAAAAGAAACCGGCTAATAGGTCAGATTCAACTTATTAACATAAAAAAAAAAATCGGCGGGTTTTTGGtt                                          >  8:286929/1‑151 (MQ=2)

                                                                                                     gATTTACGTAAAACCCGG‑CTTCGGCGGTTTTTTGCTTTTGGCGGGGCACAACCATGAATGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCttt                                         <  7:392767/151‑1 (MQ=11)

                                                                                                       tttACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACCCAAAGCTCGGAGGCTGAGCGGTCTGTTGGCCCGGGTGTCCGTAGAATCGCGGTTGGGCGggttttgggttttg                                       >  7:967867/1‑151 (MQ=0)

                                                                                                          tcGTAATAACCCGTCTGCGGCGGGTTTTTGCTTTTGTCGGCGCACAATGATGAATGACTCTACACTACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTgg                                      <  4:520131/149‑1 (MQ=11)

                                                                                                            gTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAAACCGATGCGGCGGGTTTATGATTTTTGAgtg                                  >  2:824960/1‑149 (MQ=11)

                                                                                                                cacccgg‑CTTCGGCGGTTTTTTGCGTTTGCCGGCGCCGAAAGTTGAAGGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAg                              <  5:857921/145‑1 (MQ=11)

                                                                                                                    ccgccTTCCGCGGGTTTTTGGTTTTGCCGCGGCGCCACCATGATTGACTCGCCACCACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                           <  1:642212/147‑1 (MQ=11)

                                                                                                                       tcTTCTGCGGGTTATTGTTTTTGTCGGTGTATAATTTTGAATGACTTTAGTAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCCGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAg                          <  4:173777/148‑1 (MQ=11)

                                                                                                                         ttCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAACCCGATACGGCGGGTTTAAGATTTTAGATCGGaagaaaga                         >  8:364779/1‑148 (MQ=11)

                                                                                                                          gcggcggGTTTTTGTTCTTGTCGGTGCACAACCATGAATGACTCTCCAAGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat                        <  6:533689/147‑1 (MQ=11)

                                                                                                                           cggcggGTTTTTGGTTTTGGCCGCGCACAAACTTGAAGGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                    <  7:755013/151‑1 (MQ=11)

                                                                                                                             gcggGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGGCTCTGAACACAAAACAAAGCGGCTGATCGCACTCGGGAACCTGGGGTACTTAAAAAGCCGGTTCTGATGGTTATTTATTTTGGAgggggggatagcgggaggc                  >  5:295405/1‑136 (MQ=2)

                                                                                                                              cggGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGCTGTACGTAATAACCCGTTACGGCGGGTTTAAGCTTTTAGAGGTGGAGaaaactcaataa                  >  1:826246/1‑141 (MQ=11)

                                                                                                                                  tttttGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGTTGGGGAGGCAAGATGTATGACttgc             >  3:780155/1‑148 (MQ=11)

                                                                                                                                              gAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGCTGCTGAAAGCTGAATGACtctcgccccgccgat  >  1:435696/1‑139 (MQ=11)

CCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTAT  >  DH1/613448‑613738

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