BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	96,901:1	+G	5.9%	intergenic (+205/+462)	USA300TCH1516_00088 → / ← dus_2	hypothetical protein/putative tRNA‑dihydrouridine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	96,901	1	.	G	5.9%	523.5 / 25.0	118	intergenic (+205/+462)	USA300TCH1516_00088/dus_2	hypothetical protein/putative tRNA‑dihydrouridine synthase
Reads supporting (aligned to +/- strand): ref base . (53/58); new base G (3/4); total (56/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 121 total aligned reads displayed.

GTATAAGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTATTTGAG  >  USA300TCH1516_ALE/96811‑96991

  ataAGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGat‑at                                                                                           <  1:3193547/91‑1 (MQ=255)

     aGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑Atttt                                                                                        >  1:2689060/1‑91 (MQ=255)

      gATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAATATAGAAACAGCAGTAAGAT‑ATTTTc                                                                                       <  1:1240012/91‑1 (MQ=255)

             aaaGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑Ga                                                                               >  1:2689086/1‑91 (MQ=255)

               aGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑Gaaa                                                                             <  1:2865935/91‑1 (MQ=255)

               aGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑Gaaa                                                                             <  1:2865796/91‑1 (MQ=255)

               aGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑Gaaa                                                                             >  1:136949/1‑91 (MQ=255)

               aGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAAAAACAGCAGTAAGAT‑ATTTTCTAATT‑Gaaa                                                                             >  1:2106046/1‑91 (MQ=255)

                 ttgttACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAat                                                                           <  1:4050062/91‑1 (MQ=255)

                 ttgttACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTGCTAATT‑GAAAat                                                                           <  1:3175075/91‑1 (MQ=255)

                  tgttACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAatc                                                                          >  1:2219133/1‑91 (MQ=255)

                  tgtgACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTGAGATGATTTTCAATTC‑GAAAat                                                                           <  1:3220543/91‑1 (MQ=25)

                    ttACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAatcat                                                                        >  1:3941554/1‑91 (MQ=255)

                    ttACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAatcat                                                                        >  1:2516930/1‑91 (MQ=255)

                    ttACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAatcat                                                                        >  1:2530785/1‑91 (MQ=255)

                        tGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCa                                                                    <  1:2981473/91‑1 (MQ=255)

                        tGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCa                                                                    >  1:2499819/1‑91 (MQ=255)

                          ttttAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCAct                                                                  <  1:1471452/91‑1 (MQ=255)

                          ttttAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCAct                                                                  >  1:801476/1‑91 (MQ=255)

                             tAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTGAGATGATTTTC‑AATTAGAAAAT‑ATCTTActgct                                                               <  1:3156519/91‑1 (MQ=17)

                              aaaTGATCACGATGAAGTCATTCAAGAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCActgctg                                                              <  1:3847022/91‑1 (MQ=255)

                               aaTGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGt                                                             >  1:3792674/1‑91 (MQ=255)

                                aTGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTGAGATGATTTTC‑AATTAGAAAAT‑ATCTTACTGCTGtt                                                            <  1:4037644/91‑1 (MQ=17)

                                aTGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGtt                                                            >  1:1595533/1‑91 (MQ=255)

                                 tGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGttt                                                           <  1:3958776/91‑1 (MQ=255)

                                  gATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTGAGATGATTTTC‑AATTAGAAAAT‑ATCTTACTGCTGtttc                                                          <  1:1504208/91‑2 (MQ=14)

                                  gATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGtttt                                                          <  1:3215805/91‑1 (MQ=255)

                                  gATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGttac                                                          <  1:1951703/91‑3 (MQ=255)

                                  gATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGttac                                                          <  1:1951777/91‑3 (MQ=255)

                                  gATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTGCTAATT‑GAAAATCATCTCACTGCTGtttt                                                          >  1:2305876/1‑91 (MQ=255)

                                   aTCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGttttt                                                         >  1:3364676/1‑91 (MQ=255)

                                    tCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGtttttt                                                        >  1:3361472/1‑91 (MQ=255)

                                    tCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGtttttt                                                        >  1:2488025/1‑91 (MQ=255)

                                    tCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGtttttt                                                        >  1:2488038/1‑91 (MQ=255)

                                    tCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGtttttt                                                        >  1:2488566/1‑91 (MQ=255)

                                      aCGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTaa                                                      >  1:1226860/1‑91 (MQ=255)

                                       cGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTaaa                                                     <  1:797588/91‑1 (MQ=255)

                                       cGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTaaa                                                     >  1:2532799/1‑91 (MQ=255)

                                       cGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTaaa                                                     >  1:2813559/1‑91 (MQ=255)

                                        gATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAg                                                    <  1:3081840/91‑1 (MQ=255)

                                        gATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAg                                                    >  1:602620/1‑91 (MQ=255)

                                        gATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAg                                                    <  1:1529623/91‑1 (MQ=255)

                                         aTGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAgg                                                   <  1:1945377/91‑1 (MQ=255)

                                         aTGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAgg                                                   >  1:1063001/1‑91 (MQ=255)

                                          tGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGt                                                  <  1:359040/91‑1 (MQ=255)

                                          tGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGt                                                  <  1:4006818/91‑1 (MQ=255)

                                           gAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGtt                                                 >  1:46572/1‑91 (MQ=255)

                                           gAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGtt                                                 <  1:719129/91‑1 (MQ=255)

                                             aGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTta                                               >  1:231098/1‑91 (MQ=255)

                                             aGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTta                                               <  1:1749050/91‑1 (MQ=255)

                                             aGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTta                                               >  1:1369828/1‑91 (MQ=255)

                                             aGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTta                                               <  1:2976718/91‑1 (MQ=255)

                                              gTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GACAATCATCTCACTGCTGTTTTTTAAAGGTTtat                                              >  1:3529788/1‑91 (MQ=255)

                                              gTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTtat                                              >  1:3797915/1‑91 (MQ=255)

                                              gTCATTCAATAAGAACGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GACAATCATCTCACTGCTGTTTTTTAAAGGTTtat                                              >  1:3529787/1‑91 (MQ=255)

                                                cATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATAc                                            >  1:2597192/1‑91 (MQ=255)

                                                 aTTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATAcc                                           >  1:1631796/1‑91 (MQ=255)

                                                   tCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTc                                         <  1:697061/91‑1 (MQ=255)

                                                       tAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATcc                                     >  1:3046976/1‑91 (MQ=255)

                                                       tAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATcc                                     >  1:2506789/1‑91 (MQ=255)

                                                          gAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCctct                                  >  1:1655877/1‑91 (MQ=255)

                                                           aaTGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTa                                 >  1:1216305/1‑91 (MQ=255)

                                                             tGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTaaa                               <  1:2694504/91‑1 (MQ=255)

                                                              gATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAat                              >  1:909743/1‑91 (MQ=255)

                                                              gATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAat                              <  1:3940224/91‑1 (MQ=255)

                                                              gATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAat                              <  1:4094614/91‑1 (MQ=255)

                                                              gATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAat                              <  1:1095986/91‑1 (MQ=255)

                                                                ttATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAatta                            <  1:1199603/91‑1 (MQ=255)

                                                                  aTGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAattatt                          <  1:1025461/91‑1 (MQ=255)

                                                                  aTGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAattatt                          <  1:1307729/91‑1 (MQ=255)

                                                                   tGAAAATAGAAACTGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCAGCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTAttt                         <  1:4149768/91‑1 (MQ=255)

                                                                   tGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAACTTAttt                         >  1:1665975/1‑91 (MQ=255)

                                                                   tGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTAttt                         >  1:3887639/1‑91 (MQ=255)

                                                                   tGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTAttt                         >  1:1666123/1‑91 (MQ=255)

                                                                      aaaTAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTaaa                      <  1:4144379/91‑1 (MQ=255)

                                                                      aaaTAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTaaa                      >  1:1541083/1‑91 (MQ=255)

                                                                       aaTAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTaaaa                     >  1:111025/1‑91 (MQ=255)

                                                                        aTAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTgaaaa                    <  1:3274771/91‑1 (MQ=255)

                                                                        aTAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTaaaaa                    <  1:765862/91‑1 (MQ=255)

                                                                          aGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAaata                  <  1:1013843/91‑1 (MQ=255)

                                                                             aaCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAAtaatt               >  1:3391066/1‑91 (MQ=255)

                                                                             aaCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAAtaatt               <  1:2670649/91‑1 (MQ=255)

                                                                             aaCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAAtaatt               <  1:415364/91‑1 (MQ=255)

                                                                             aaCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAAtaatt               <  1:2726724/91‑1 (MQ=255)

                                                                              aCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAAtaatta              <  1:3697885/91‑1 (MQ=255)

                                                                              aCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAAtaatta              >  1:593223/1‑91 (MQ=255)

                                                                               cagcagTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAAtaattaa             >  1:48245/1‑91 (MQ=255)

                                                                                agcagTGAGATGATTTTCAATTA‑GAAAAT‑ATCTTACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAGAAAtaattaat            >  1:3768473/1‑91 (MQ=17)

                                                                                agcagTGAGATGATTTTCAATTA‑GAAAAT‑ATCTTACTGCTGTTTTTTAAAGGTTTATACCTCATCATCTAAATTATTTAAAAAtaattaat            >  1:1809542/1‑91 (MQ=17)

                                                                                agcagTGAGATGATTTTC‑AATTAGAAAAT‑ATCTTACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTATAAAtaattaat            >  1:3768397/1‑91 (MQ=17)

                                                                                 gcagTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCGCTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATg           <  1:3680317/91‑1 (MQ=255)

                                                                                  cagTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATgg          <  1:3259617/91‑1 (MQ=255)

                                                                                   agTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGt         >  1:4118379/1‑91 (MQ=255)

                                                                                    gTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTa        >  1:2820169/1‑91 (MQ=255)

                                                                                    gTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTa        <  1:2196955/91‑1 (MQ=255)

                                                                                       aGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTAttt     <  1:2886475/91‑1 (MQ=255)

                                                                                        gAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTATTTg    <  1:1448953/91‑1 (MQ=255)

                                                                                        gAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTATTTg    <  1:1449186/91‑1 (MQ=255)

                                                                                        gAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTATTTg    <  1:1449234/91‑1 (MQ=255)

                                                                                          t‑atTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTATTTGAg  <  1:623841/91‑1 (MQ=255)

GTATAAGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTATTTGAG  >  USA300TCH1516_ALE/96811‑96991

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