BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	DH1	1,971,067	A→C	K9Q (AAG→CAG)	cspC →	cold‑shock protein CspC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	1,971,067	0	A	C	98.9%	1759.0 / ‑1.5	451	K9Q (AAG→CAG)	cspC	cold‑shock protein CspC
Reads supporting (aligned to +/- strand): ref base A (3/2); new base C (236/210); total (239/212)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 451 total aligned reads displayed.

AAACACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTAAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGCC  >  DH1/1970920‑1971214

 aaCACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTg                                                                                                                                                 >  7:478660/1‑151 (MQ=255)

   cacGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGt                                                                                                                                               >  6:1637423/1‑151 (MQ=255)

      gttgacgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCa                                                                                                                                            <  5:470321/151‑1 (MQ=255)

          acgttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACg                                                                                                                                         <  1:643743/150‑1 (MQ=255)

            gttgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGt                                                                                                                                      >  5:64196/1‑151 (MQ=255)

              tgaATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCt                                                                                                                                    >  5:1274544/1‑151 (MQ=255)

                aaTACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTa                                                                                                                                   >  7:389289/1‑150 (MQ=255)

                   aCCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCTTGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAgg                                                                                                                               <  7:1676271/151‑1 (MQ=255)

                     cGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGtt                                                                                                                             <  7:1524958/151‑1 (MQ=255)

                      gCTTCTCTGCCGAGCCTTATATTGGTGCATCATGCAGTAATGTGTCAGTTTTATCTATCTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTAAAATGGCAAAGATTAAAGGTCAGGTTCAGAGGTTCAACGAGTCTAAAGGttt                                                                                                                            >  2:1208343/1‑151 (MQ=255)

                        ttCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTg                                                                                                                          >  8:147920/1‑151 (MQ=255)

                          ctctGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGc                                                                                                                        >  2:1074825/1‑151 (MQ=255)

                            ctGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGc                                                                                                                        <  6:14196/149‑1 (MQ=255)

                              gCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCCGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCt                                                                                                                       <  8:154959/148‑1 (MQ=255)

                              gCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCt                                                                                                                       <  8:24129/148‑1 (MQ=255)

                              gCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTc                                                                                                                     <  4:883496/150‑1 (MQ=255)

                              gCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCa                                                                                                                    <  2:1590723/151‑1 (MQ=255)

                              gCCGAGCCTTATATTGGGGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCtt                                                                                                                      <  6:852542/149‑1 (MQ=255)

                               ccGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCAt                                                                                                                   >  6:211043/1‑151 (MQ=255)

                               ccGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCAt                                                                                                                   >  7:1361817/1‑151 (MQ=255)

                                 gAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTTTGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCAtt                                                                                                                  >  7:1588244/1‑150 (MQ=255)

                                 gAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCAtt                                                                                                                  >  1:1409527/1‑150 (MQ=255)

                                 gAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTa                                                                                                                 <  3:1309996/151‑1 (MQ=255)

                                    ccTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTc                                                                                                              >  3:1395382/1‑151 (MQ=255)

                                      ttATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCg                                                                                                            >  4:1029001/1‑151 (MQ=255)

                                       tttatTGGTGCCTCATGCAGTAATGTGTCAGTTTTATGTATGTTATGCCTGCGGGCGAATAATAGAATCTAAGTAATTTTTCAAATTGGAAATATTAATGGTCATTTTCAGTTGTTAAACGAGTCTAAAGGTTTTTTCTTCATTAATCCg                                                                                                            <  2:1358018/148‑1 (MQ=255)

                                       tataTTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTcc                                                                                                             <  2:1660278/149‑1 (MQ=255)

                                       tataTTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCgg                                                                                                           >  3:404900/1‑151 (MQ=255)

                                          atTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCt                                                                                                         <  8:1627007/150‑1 (MQ=255)

                                          atTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCt                                                                                                         <  1:1032462/150‑1 (MQ=255)

                                           ttGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCGACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTgg                                                                                                       >  6:1097039/1‑150 (MQ=255)

                                            tgttGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGa                                                                                                       <  4:1357781/147‑1 (MQ=255)

                                             ggTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATg                                                                                                     <  4:1223143/151‑1 (MQ=255)

                                                gCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGgc                                                                                                   <  5:1676415/150‑1 (MQ=255)

                                                 ccTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATATAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGgc                                                                                                   >  4:597969/1‑149 (MQ=255)

                                                       gCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCaa                                                                                              >  7:203654/1‑148 (MQ=255)

                                                       gCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGa                                                                                           <  3:1070121/151‑1 (MQ=255)

                                                             aTGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGtt                                                                                     <  3:477558/151‑1 (MQ=255)

                                                                 gtCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGt                                                                                  >  1:1255038/1‑150 (MQ=255)

                                                                  tCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTac                                                                                >  7:908520/1‑151 (MQ=255)

                                                                  tCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTac                                                                                >  6:782083/1‑151 (MQ=255)

                                                                   cAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTac                                                                                <  5:1466989/150‑1 (MQ=255)

                                                                   cAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTac                                                                                >  2:953898/1‑150 (MQ=255)

                                                                      ttttATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACtt                                                                            >  5:787570/1‑151 (MQ=255)

                                                                         tATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTtctc                                                                         >  8:704959/1‑151 (MQ=255)

                                                                          atctatctTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATCTCTTCGTACACTTCTcc                                                                        >  6:1363178/1‑151 (MQ=255)

                                                                             tATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCt                                                                     >  7:635067/1‑151 (MQ=255)

                                                                             tATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCt                                                                     >  8:1247196/1‑151 (MQ=255)

                                                                              aTGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCt                                                                     >  2:304651/1‑150 (MQ=255)

                                                                              aTGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCt                                                                     >  8:405238/1‑150 (MQ=255)

                                                                              aTGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTa                                                                    <  1:970162/151‑1 (MQ=255)

                                                                                gTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATc                                                                  <  7:939342/151‑1 (MQ=255)

                                                                                     gCCTTCGGGCGAAGAAAACAATCTAAGCAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTACGCTATCCAgcg                                                             >  4:164907/1‑149 (MQ=255)

                                                                                            ggCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGc                                                      >  5:905284/1‑151 (MQ=255)

                                                                                             gCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCt                                                     >  4:291580/1‑151 (MQ=255)

                                                                                             gCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCt                                                     <  5:1569597/151‑1 (MQ=255)

                                                                                                aagaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTc                                                   >  8:1084122/1‑150 (MQ=255)

                                                                                                aagaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCa                                                  <  6:341217/151‑1 (MQ=255)

                                                                                                 agaaAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCa                                                  >  3:923409/1‑150 (MQ=255)

                                                                                                  gaaAACAATCTAAGGAATTTTTCAAATGCCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCaaa                                                <  4:1507399/151‑1 (MQ=255)

                                                                                                     aaCAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAAct                                             >  6:550791/1‑151 (MQ=255)

                                                                                                     aaCAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAAct                                             >  1:118139/1‑151 (MQ=255)

                                                                                                     aaCAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTACAGTGGTTCAACGAGTCAAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCTAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGGTTCAAAAc                                              >  4:953892/1‑150 (MQ=255)

                                                                                                        aaTCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTg                                          <  7:1152245/151‑1 (MQ=255)

                                                                                                        aaTCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTg                                          <  7:702499/151‑1 (MQ=255)

                                                                                                          tCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGc                                        >  3:561674/1‑151 (MQ=255)

                                                                                                          tCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGc                                        >  4:802557/1‑151 (MQ=255)

                                                                                                            tAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTg                                      >  3:182906/1‑151 (MQ=255)

                                                                                                            tAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTg                                      >  3:828083/1‑151 (MQ=255)

                                                                                                             aaGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTGCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTg                                      >  3:1492854/1‑150 (MQ=255)

                                                                                                             aaGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTgg                                         <  2:666576/147‑1 (MQ=255)

                                                                                                             aaGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTg                                      <  5:976478/150‑1 (MQ=255)

                                                                                                             aaGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTg                                      >  6:1179327/1‑150 (MQ=255)

                                                                                                               ggAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCATAACTCTGGCTGa                                     >  2:855328/1‑149 (MQ=255)

                                                                                                                 aaTTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGt                                 >  2:754291/1‑151 (MQ=255)

                                                                                                                 aaTTTTTCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGt                                 >  3:783505/1‑151 (MQ=255)

                                                                                                                     tttCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTc                                >  5:319545/1‑148 (MQ=255)

                                                                                                                     tttCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTc                                >  2:572751/1‑148 (MQ=255)

                                                                                                                       tCAAATGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAAc                           >  7:552574/1‑151 (MQ=255)

                                                                                                                         aaaTGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGt                         <  3:647667/151‑1 (MQ=255)

                                                                                                                           aTGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTg                       <  7:807603/151‑1 (MQ=255)

                                                                                                                           aTGGCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTg                       >  4:400718/1‑151 (MQ=255)

                                                                                                                            tGGCAAAGATTAAAGGTCAGGTTCAGTGCTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGCTGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTg                       >  6:951399/1‑150 (MQ=255)

                                                                                                                              gCAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGt                    >  1:1548642/1‑151 (MQ=255)

                                                                                                                               cAAAGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGtt                   <  2:1495634/151‑1 (MQ=255)

                                                                                                                                 aacaTTAAAGGTAAGGTTCACTCGTTCACCGACTCTAAAGGTATAGGCTTCATTACTCACGCTCATGGCATCAAAGATCTGTTCGTACACTTCTCCGCTATCCAGCGTAAAGCATTCAAAACTCTCGCTGAAGCACAGAACGTAGAGATcc                 >  6:1602526/4‑150 (MQ=255)

                                                                                                                                 aaGATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCg                 >  6:1281863/1‑151 (MQ=255)

                                                                                                                                   gATTAAAGGTCAGGTTCTGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGGAATGATGTGTTCTTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAAGGTTGAGTTCGa                <  8:875883/150‑1 (MQ=255)

                                                                                                                                   gATTAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGa                >  4:531869/1‑150 (MQ=255)

                                                                                                                                     ttAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAAt             >  7:286883/1‑151 (MQ=255)

                                                                                                                                     ttAAAGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAAt             >  1:242697/1‑151 (MQ=255)

                                                                                                                                        aaGGTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCa          <  6:611849/151‑1 (MQ=255)

                                                                                                                                          ggTCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAgg        <  7:1207465/151‑1 (MQ=255)

                                                                                                                                            tCAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGAc      >  1:71658/1‑151 (MQ=255)

                                                                                                                                             cAGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACg     <  8:74528/151‑1 (MQ=255)

                                                                                                                                              aGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGAc      >  7:817763/1‑149 (MQ=255)

                                                                                                                                              aGGTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACgg    >  1:899277/1‑151 (MQ=255)

                                                                                                                                               ggTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATTTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGc   >  3:1161366/1‑151 (MQ=255)

                                                                                                                                               ggTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGc   >  8:516662/1‑151 (MQ=255)

                                                                                                                                                gTTCAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACCCTTCTCCGCTATCCAAGGTAAGGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGcc  >  4:314801/1‑151 (MQ=255)

AAACACGTTGACGTTGAATACCGCTTCTCTGCCGAGCCTTATATTGGTGCCTCATGCAGTAATGTGTCAGTTTTATCTATGTTATGCCTGCGGGCGAAGAAAACAATCTAAGGAATTTTTCAAATGGCAAAGATTAAAGGTCAGGTTAAGTGGTTCAACGAGTCTAAAGGTTTTGGCTTCATTACTCCGGCTGATGGCAGCAAAGATGTGTTCGTACACTTCTCCGCTATCCAGGGTAATGGCTTCAAAACTCTGGCTGAAGGTCAGAACGTTGAGTTCGAAATTCAGGACGGCC  >  DH1/1970920‑1971214

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