BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	DH1	3,871,936	G→A	S51F (TCT→TTT)	yjtD ←	RNA methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	3,871,936	0	G	A	98.7%	1924.2 / 6.6	522	S51F (TCT→TTT)	yjtD	RNA methyltransferase
Reads supporting (aligned to +/- strand): ref base G (2/5); new base A (304/210); total (306/216)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 522 total aligned reads displayed.

TCCTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAGATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTTATAC  >  DH1/3871789‑3872083

tCCTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAAt                                                                                                                                                   <  5:1634237/150‑1 (MQ=255)

     tAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCAtgt                                                                                                                                             >  2:1232868/1‑151 (MQ=255)

      aaCAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCAtgtg                                                                                                                                            >  6:1120543/1‑151 (MQ=255)

           cGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAGACCCATGTGCGAcc                                                                                                                                       <  8:398575/151‑1 (MQ=255)

             gCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGAccc                                                                                                                                      >  8:1018297/1‑150 (MQ=255)

              cACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGAccc                                                                                                                                      >  7:533958/1‑149 (MQ=255)

                ccAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCg                                                                                                                                  <  7:1016321/151‑1 (MQ=255)

                ccAGTTCAACTGGCGTGGAGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCg                                                                                                                                  <  7:681889/151‑1 (MQ=255)

                  aGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCggg                                                                                                                                <  4:1406915/151‑1 (MQ=255)

                   gTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGt                                                                                                                               >  7:416572/1‑151 (MQ=255)

                     tCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtg                                                                                                                              <  5:1171718/150‑1 (MQ=255)

                       aaCTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggct                                                                                                                           <  5:787372/151‑1 (MQ=255)

                       aaCTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggct                                                                                                                           <  7:1244112/151‑1 (MQ=255)

                        aCTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACATTAATATTATCAATAATATCACCAAATACATGTGCGACCCAGCgggagact                                                                                                                           >  2:37448/1‑145 (MQ=255)

                         cTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggctgg                                                                                                                         >  7:749881/1‑151 (MQ=255)

                          tggcgtggcgTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCAATGTGCGACCAAGCGGCtggctggc                                                                                                                        >  3:1611376/1‑151 (MQ=255)

                           ggcgtgtcgtTGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTTAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggctgg                                                                                                                         <  6:6213/149‑1 (MQ=255)

                              ggggcgtAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGACGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCa                                                                                                                    <  2:1599292/149‑1 (MQ=255)

                                ggcgtAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAgg                                                                                                                  >  4:737196/1‑151 (MQ=255)

                                  cgtAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGgtg                                                                                                                <  1:213567/151‑1 (MQ=255)

                                   gtagtaATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGgtgt                                                                                                               >  8:239643/1‑151 (MQ=255)

                                   gtagtaATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGCGGCTGGCTCCAGgtgt                                                                                                               >  7:1816108/1‑151 (MQ=255)

                                       taATGATATATAGCCCCACAGCGCGCAGAGCAGGCGACAGTGAAATCGACATCGTCTAACGATTCAGCAAATCTCGGCAAAACATTAATATTATCAAAAATAACAACAAAGCAAAGTGCGTCCAAGCGAGAGGATGCCTCCAggatagact                                                                                                           >  2:972370/1‑144 (MQ=255)

                                             tatTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTAt                                                                                                     >  1:1656671/1‑151 (MQ=255)

                                                ttCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCg                                                                                                   <  2:820534/150‑1 (MQ=255)

                                                ttCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCg                                                                                                   >  2:177080/1‑150 (MQ=255)

                                                  cGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGAc                                                                                                 <  3:150195/150‑1 (MQ=255)

                                                  cGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACa                                                                                                <  5:1842877/151‑1 (MQ=255)

                                                    cccGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAAt                                                                                              >  4:413021/1‑151 (MQ=255)

                                                      cGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATcc                                                                                            <  7:1372966/151‑1 (MQ=255)

                                                      cGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATcc                                                                                            <  8:46806/151‑1 (MQ=255)

                                                       gACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCg                                                                                           >  3:227576/1‑151 (MQ=255)

                                                        aCTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGc                                                                                          >  6:454069/1‑151 (MQ=255)

                                                         cTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGc                                                                                          >  3:997726/1‑150 (MQ=255)

                                                           gcgcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGc                                                                                          >  5:72339/1‑148 (MQ=255)

                                                            cgcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGc                                                                                          >  1:153254/1‑147 (MQ=255)

                                                            cgcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGc                                                                                          >  3:886619/1‑147 (MQ=255)

                                                            cgcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGAt                                                                                      <  8:749118/151‑1 (MQ=255)

                                                             gcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGc                                                                                          >  6:375552/1‑146 (MQ=255)

                                                             gcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATc                                                                                     <  8:1337354/151‑1 (MQ=255)

                                                               gcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGc                                                                                   >  5:1037/1‑151 (MQ=255)

                                                                    ggtggCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAAc                                                                              <  5:941274/151‑1 (MQ=255)

                                                                    ggtggCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAAc                                                                              <  4:1687332/151‑1 (MQ=255)

                                                                     gtgcCGCCAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCCGATCGCTAAAcc                                                                             >  8:1673612/1‑151 (MQ=255)

                                                                      tggCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAAccc                                                                            >  8:148721/1‑151 (MQ=255)

                                                                       ggCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGCTCGCTAAAcccc                                                                           >  7:1518131/1‑151 (MQ=255)

                                                                       ggCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAAcccc                                                                           >  3:1028649/1‑151 (MQ=255)

                                                                         cGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCAt                                                                         >  3:1675160/1‑151 (MQ=255)

                                                                          gACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGCTGTCCCTGAGAATAGACAATCCGCAGATCGCTAAACCCAATc                                                                        >  6:744296/1‑151 (MQ=255)

                                                                           aCAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAACCCCCATCg                                                                       >  1:1360759/1‑151 (MQ=255)

                                                                             aGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGtt                                                                     <  4:1126890/151‑1 (MQ=255)

                                                                             aGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTCGCTCCAGGTGTGCCTGACTATCGACAATCCCCAGATCGCTAAACCCCATCGtt                                                                     >  6:869925/1‑151 (MQ=255)

                                                                               tGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGtttt                                                                   >  7:1189965/1‑151 (MQ=255)

                                                                                gAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTc                                                                  >  2:702630/1‑151 (MQ=255)

                                                                                 aaaTCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTc                                                                  >  4:197296/1‑150 (MQ=255)

                                                                                  aaTCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCAt                                                                <  7:1631672/151‑1 (MQ=255)

                                                                                  aaTCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCAt                                                                <  8:187887/151‑1 (MQ=255)

                                                                                   aTCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAAACCATGTGCGACCCAGCGGCTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCAtt                                                               >  5:334093/1‑151 (MQ=255)

                                                                                       aCATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGccc                                                           >  3:1115156/1‑151 (MQ=255)

                                                                                        cATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCcg                                                          <  8:1192618/151‑1 (MQ=255)

                                                                                        aaTCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCcg                                                          <  7:947734/150‑1 (MQ=255)

                                                                                           cGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCgcgc                                                       <  3:480404/151‑1 (MQ=255)

                                                                                            gtgtAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGcgcc                                                      <  3:409373/151‑1 (MQ=255)

                                                                                                aaCGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCt                                                   <  2:1521205/150‑1 (MQ=255)

                                                                                                aaCGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCt                                                   <  1:667678/150‑1 (MQ=255)

                                                                                                 aCGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGc                                                 >  1:1108301/1‑151 (MQ=255)

                                                                                                  cGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAGTATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGcc                                                >  2:1343815/1‑151 (MQ=255)

                                                                                                      tCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCa                                             >  6:231393/1‑150 (MQ=255)

                                                                                                       cAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAat                                           >  3:1234911/1‑151 (MQ=255)

                                                                                                         gCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAatat                                         >  5:1791404/1‑151 (MQ=255)

                                                                                                          ccAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATAtt                                        >  7:1104293/1‑151 (MQ=255)

                                                                                                           cAATGTCGGGTAAACTTTAATATTATCAATAATATCACAAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGAGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCAACCTTTTCATTGCCCGCGCCGCTGCCCCAATAttt                                       >  3:776277/1‑151 (MQ=255)

                                                                                                            aaTGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATGCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATAtttt                                      <  5:1526826/151‑1 (MQ=255)

                                                                                                            aaTGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCCTCGTTTTCATTGCCCGCGCCGCTGCCCCAATAtttt                                      >  1:866246/1‑151 (MQ=255)

                                                                                                            aaTGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATAtttt                                      >  8:1541472/1‑151 (MQ=255)

                                                                                                               gTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTg                                   <  8:1096604/151‑1 (MQ=255)

                                                                                                                 cGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGct                                 >  1:416896/1‑151 (MQ=255)

                                                                                                                 cGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGct                                 >  7:323676/1‑151 (MQ=255)

                                                                                                                    gAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCGGGCTCCAGGTGGGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTg                              <  3:292815/151‑1 (MQ=255)

                                                                                                                     aaaaCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTgg                             >  6:1363483/1‑151 (MQ=255)

                                                                                                                      aaaCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGc                            <  1:1000492/151‑1 (MQ=255)

                                                                                                                       aaCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGct                           >  2:472705/1‑150 (MQ=255)

                                                                                                                       aaCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCg                           <  5:57996/151‑1 (MQ=255)

                                                                                                                        aCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCgg                          >  6:1744464/1‑151 (MQ=255)

                                                                                                                        aCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCgg                          >  5:837175/1‑151 (MQ=255)

                                                                                                                             aaTATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCg                     <  5:1130040/151‑1 (MQ=255)

                                                                                                                              atatTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCgc                    >  1:558570/1‑150 (MQ=255)

                                                                                                                                  tatCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGAcc                  >  3:743062/1‑149 (MQ=255)

                                                                                                                                   atCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGAcc                  >  2:226817/1‑148 (MQ=255)

                                                                                                                                    tCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGGTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGa               >  6:1522183/1‑150 (MQ=255)

                                                                                                                                     cAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGaat             <  1:787354/151‑1 (MQ=255)

                                                                                                                                      aataatATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGaat             >  5:1488898/1‑150 (MQ=255)

                                                                                                                                       ataataTCACCAACTCCATGTGCGACCCAGCGGGAGGCTGGCTCCACGACAGCCTGACTATCGACAATCCGCAGATCGCTACACCCAATCGTTTTCATTCCCCCCTCCGATGCCCCAATATTTACTGCTCAGGAGTGTGCGACCAGaata            >  8:690793/1‑150 (MQ=255)

                                                                                                                                        taataTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTCCCCCAATATTTTCTGCTCTGGCGGCTGCGACCAGaaaaa           >  2:625100/1‑150 (MQ=255)

                                                                                                                                         aataTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATc         >  4:592953/1‑151 (MQ=255)

                                                                                                                                         aataTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATc         <  8:1397216/151‑1 (MQ=255)

                                                                                                                                          atatCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCg        <  4:267095/151‑1 (MQ=255)

                                                                                                                                           tatCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGt       >  5:509507/1‑151 (MQ=255)

                                                                                                                                            atCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGtt      <  3:65663/151‑1 (MQ=255)

                                                                                                                                                ccAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGCCGGGTGCGACCAGAATAATCGTTATAc  >  1:124458/1‑151 (MQ=255)

TCCTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAGATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTTATAC  >  DH1/3871789‑3872083

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