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	99.4%	2029.2 / ‑3.8	539	S51F (TCT→TTT)	yjtD	RNA methyltransferase
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (304/231); total (306/233)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 539 total aligned reads displayed.

TTCCTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAGATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTTATAC  >  DH1/3871788‑3872083

  ccTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATcc                                                                                                                                                 >  8:1631373/1‑151 (MQ=255)

  ccTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATcc                                                                                                                                                 <  6:1466141/151‑1 (MQ=255)

  ccTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATc                                                                                                                                                  <  8:430143/150‑1 (MQ=255)

     ctAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCa                                                                                                                                                <  7:431642/149‑1 (MQ=255)

     ctAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCa                                                                                                                                                <  6:1304969/149‑1 (MQ=255)

     ctAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCCtg                                                                                                                                              >  1:1814243/1‑151 (MQ=255)

       aaCAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCAtg                                                                                                                                              <  4:1066768/149‑1 (MQ=255)

         cAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCg                                                                                                                                          >  3:1216969/1‑151 (MQ=255)

           gCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGAc                                                                                                                                        >  6:854335/1‑151 (MQ=255)

             ggCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGAccc                                                                                                                                      >  6:1054126/1‑151 (MQ=255)

                 ccAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAATGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCg                                                                                                                                  <  5:1019465/151‑1 (MQ=255)

                  cAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCgg                                                                                                                                 <  2:1915014/151‑1 (MQ=255)

                  cAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATAAATAATATCACCAAATCCATGTGCGACCCAGCgg                                                                                                                                 >  1:874399/1‑151 (MQ=255)

                    gTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGt                                                                                                                               <  4:1185974/151‑1 (MQ=255)

                     ttCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtg                                                                                                                              >  8:1141533/1‑151 (MQ=255)

                     ttCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtg                                                                                                                              >  7:1673685/1‑151 (MQ=255)

                      tCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtgg                                                                                                                             >  1:324396/1‑151 (MQ=255)

                       cAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggc                                                                                                                            <  6:1950311/151‑1 (MQ=255)

                        aaCTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggct                                                                                                                           <  2:654876/151‑1 (MQ=255)

                         aCTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggctg                                                                                                                          <  2:1606415/151‑1 (MQ=255)

                            ggcgtggcgtAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGtggctggc                                                                                                                        <  2:413224/150‑1 (MQ=255)

                               gtggcgtAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCa                                                                                                                    <  8:258720/151‑1 (MQ=255)

                                tggcgtAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTAAAAATAATATCACCAAATCCATGTGCGCCCCAGCGGGTGGCTGGCTcc                                                                                                                     >  3:1416662/1‑149 (MQ=255)

                                   cgtAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGgtg                                                                                                                >  6:1633329/1‑151 (MQ=255)

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

                                    gtagtaATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGgtgt                                                                                                               <  1:577815/151‑1 (MQ=255)

                                      agtaATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGtgtg                                                                                                              <  4:1575453/150‑1 (MQ=255)

                                      agtaATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGc                                                                                                             >  7:694500/1‑151 (MQ=255)

                                         aaTGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCCTGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTg                                                                                                          <  1:1382369/151‑1 (MQ=255)

                                             atatTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACt                                                                                                       >  5:866700/1‑150 (MQ=255)

                                             atatTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTa                                                                                                      <  7:183348/151‑1 (MQ=255)

                                                  tCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCTGGTGTGCCTGACTATCGAc                                                                                                 >  3:1106849/1‑151 (MQ=255)

                                                  tCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGAc                                                                                                 >  2:121004/1‑151 (MQ=255)

                                                     cccGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAAt                                                                                              <  6:318803/151‑1 (MQ=255)

                                                     cccGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAAt                                                                                              >  7:1031492/1‑151 (MQ=255)

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

                                                        gACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCg                                                                                           >  5:1489531/1‑151 (MQ=255)

                                                         aCTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCTGGGGTGCCTGACTATCGACAATCCGc                                                                                          >  4:830166/1‑151 (MQ=255)

                                                           tGCGCGCAGTGGTGGCGACAGTGAAATCGACATCCTGTAACGATTAAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCAACAAATCCATGTCCGACCCAGCGCGAGGCTCGCTCCAGCAGTCCCTGACTATCGACAATCAGCat                                                                                        >  6:781035/1‑150 (MQ=255)

                                                            gcgcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCa                                                                                         <  7:361369/149‑1 (MQ=255)

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

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

                                                              gcgcAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATc                                                                                     <  2:1024727/151‑1 (MQ=255)

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

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

                                                                     ggtggCGACTGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTaaa                                                                               <  5:517904/150‑1 (MQ=255)

                                                                          cGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCAt                                                                         <  4:35910/151‑1 (MQ=255)

                                                                             cAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGt                                                                      >  6:396436/1‑151 (MQ=255)

                                                                              aGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGtt                                                                     <  2:296357/151‑1 (MQ=255)

                                                                                 gAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGtttt                                                                   <  2:763957/150‑1 (MQ=255)

                                                                                       gACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGcc                                                            >  4:488/1‑151 (MQ=255)

                                                                                       gACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGcc                                                            >  3:896246/1‑151 (MQ=255)

                                                                                         cATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCcg                                                          <  4:1831280/151‑1 (MQ=255)

                                                                                         cATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCACATCGGTAAACCCCATCGTTTTCATTGCCcg                                                          >  8:890522/1‑151 (MQ=255)

                                                                                         cATCGTGTAACGATTCAGCCAATCACGGGAAAACTTTAATATTATCAATAATATCACCAAATCAATGTGCGACCAAGCCGCTGGCTCCATCCAGCTGTCACTGAATATAGACAATACGCAGATAGATAAACCCAATCGTTTTCAtttcaac                                                          >  4:651043/1‑146 (MQ=255)

                                                                                          aTCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCcgc                                                         >  8:1315452/1‑151 (MQ=255)

                                                                                           tCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCcgcg                                                        >  7:527784/1‑151 (MQ=255)

                                                                                             gtgtAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGcgcc                                                      <  4:1364279/151‑1 (MQ=255)

                                                                                              tgtAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGcgccg                                                     >  3:1001348/1‑151 (MQ=255)

                                                                                               gtAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATGTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGcgccg                                                     <  5:1569045/150‑1 (MQ=255)

                                                                                               gtAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGcgccgc                                                    >  5:194412/1‑151 (MQ=255)

                                                                                               gtAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGcgccgc                                                    >  5:688371/1‑151 (MQ=255)

                                                                                                      ttCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTCGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTCCCCGCGCCGCTGcccc                                              >  4:884831/1‑150 (MQ=255)

                                                                                                         aGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAat                                           >  7:1263103/1‑150 (MQ=255)

                                                                                                           ccAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTCTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATAtt                                        <  1:1197791/151‑1 (MQ=255)

                                                                                                            cAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATAttt                                       >  3:261837/1‑151 (MQ=255)

                                                                                                            cAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATAttt                                       >  6:949237/1‑151 (MQ=255)

                                                                                                                gTCGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTg                                   >  3:460088/1‑151 (MQ=255)

                                                                                                                  cGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGct                                 <  5:1189276/151‑1 (MQ=255)

                                                                                                                  cGGGAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGc                                  <  5:173617/150‑1 (MQ=255)

                                                                                                                     gAAAACTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTg                              <  4:195385/151‑1 (MQ=255)

                                                                                                                      aaaaCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTgg                             >  7:1392675/1‑151 (MQ=255)

                                                                                                                        aaCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCg                           >  2:1306570/1‑151 (MQ=255)

                                                                                                                        aaCATTAATATTATCAATAAAATCACCAAATCCATGTCCGACCCAGCGCGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAAAATTTTCTGGTCTGGCg                           >  8:1449079/1‑151 (MQ=255)

                                                                                                                         aCTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCgg                          >  8:1215286/1‑151 (MQ=255)

                                                                                                                          cTTTAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCggg                         >  6:1776572/1‑151 (MQ=255)

                                                                                                                          atttAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCggg                         <  4:297830/150‑1 (MQ=255)

                                                                                                                           tttAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGt                        >  6:573337/1‑151 (MQ=255)

                                                                                                                           tttAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGt                        >  4:1660582/1‑151 (MQ=255)

                                                                                                                            ttAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCggg                         <  2:527824/149‑1 (MQ=255)

                                                                                                                             tAATATTATCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTg                       <  4:1153192/150‑1 (MQ=255)

                                                                                                                                 attatCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAACCCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGAcc                  >  3:1530289/1‑151 (MQ=255)

                                                                                                                                 attatCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGccc                  >  7:902499/1‑151 (MQ=255)

                                                                                                                                 attatCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGccc                  >  3:693307/1‑151 (MQ=255)

                                                                                                                                 attatCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGAcc                  <  5:1144021/151‑1 (MQ=255)

                                                                                                                                    atCAATAATATCCCCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGa               <  5:38591/151‑1 (MQ=255)

                                                                                                                                    atCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGGGAcc                  >  1:42555/1‑148 (MQ=255)

                                                                                                                                    atCAATAATATCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATAGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGAcc                  >  3:395006/1‑148 (MQ=255)

                                                                                                                                        ataataTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGaata            >  4:1494772/1‑150 (MQ=255)

                                                                                                                                         taataTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGaataat          >  6:1524407/1‑151 (MQ=255)

                                                                                                                                         taataTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGaataat          >  6:137970/1‑151 (MQ=255)

                                                                                                                                          aataTCACCAGATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGGAGATCGCTAAACCCCTTCGTTTTCATTGCCCGCGCCGCTGCCGCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATc         <  2:1528346/151‑1 (MQ=255)

                                                                                                                                          aataTCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATc         >  3:1620564/1‑151 (MQ=255)

                                                                                                                                           atatCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCg        >  3:363129/1‑151 (MQ=255)

                                                                                                                                              tCACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGCCAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGtt      >  7:146033/1‑150 (MQ=255)

                                                                                                                                               cACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTtat    <  3:1826738/151‑1 (MQ=255)

                                                                                                                                               cACCAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTtat    >  7:1574183/1‑151 (MQ=255)

                                                                                                                                                 ccAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTTATAc  <  1:1594330/151‑1 (MQ=255)

                                                                                                                                                 ccAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTTATAc  >  5:812663/1‑151 (MQ=255)

                                                                                                                                                 ccAAATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCACCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTTATAc  >  5:1867443/1‑151 (MQ=255)

TTCCTCTAACAGCGGCACCAGTTCAACTGGCGTGGCGTAGTAATGATATTTCGCCCGACTGCGCGCAGTGGTGGCGACAGTGAAATCGACATCGTGTAACGATTCAGCCAATGTCGGGAAAACTTTAATATTATCAATAATATCACCAGATCCATGTGCGACCCAGCGGGTGGCTGGCTCCAGGTGTGCCTGACTATCGACAATCCGCAGATCGCTAAACCCCATCGTTTTCATTGCCCGCGCCGCTGCCCCAATATTTTCTGCTCTGGCGGGTGCGACCAGAATAATCGTTATAC  >  DH1/3871788‑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: `‑`