BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	96,897	A→G	6.1%	intergenic (+201/+466)	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,897	0	A	G	6.1%	385.8 / 10.3	115	intergenic (+201/+466)	USA300TCH1516_00088/dus_2	hypothetical protein/putative tRNA‑dihydrouridine synthase
Reads supporting (aligned to +/- strand): ref base A (53/55); new base G (3/4); total (56/59)
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.90e-01

Only 100 of 117 total aligned reads displayed.

AATAGTATAAGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTA  >  USA300TCH1516_ALE/96807‑96985

aaTAGTATAAGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTa                                                                                            >  1:1099689/1‑91 (MQ=255)

 aTAGTATAAGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACGGCAGTaa                                                                                           >  1:2065302/1‑91 (MQ=255)

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

        aaGATTGAAAAAGTTGTTACTGTTTCAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑Attt                                                                                   <  1:3868520/91‑1 (MQ=255)

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

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

               aaaaaGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAAtt                                                                            <  1:2668467/91‑1 (MQ=255)

                aaaaGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑g                                                                          >  1:2511582/1‑91 (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)

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

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

                     ttgttACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTCCTAATT‑GAAAat                                                                     <  1:1620908/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:2794567/1‑91 (MQ=255)

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

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

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

                         tACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAatcatc                                                                 >  1:758230/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:1471305/91‑1 (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)

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

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

                                     tGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGttt                                                     >  1:2427302/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:2488025/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:2813559/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                 aGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTta                                         >  1:231098/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)

                                                  gTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTtat                                        <  1:557873/91‑1 (MQ=255)

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

                                                  gTCATTCAATAAGAACGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTtat                                        >  1:3528990/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:3751424/91‑1 (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)

                                                            aaGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCct                              <  1:3959230/91‑1 (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:3940224/91‑1 (MQ=255)

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

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

                                                                  gATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAat                        >  1:909743/1‑91 (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:1666123/1‑91 (MQ=255)

                                                                       tGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTAttt                   >  1:3887639/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)

                                                                             tAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAaat             >  1:3248070/1‑91 (MQ=255)

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

                                                                                aaaCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAaataat          <  1:1704779/91‑1 (MQ=255)

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

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

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

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

                                                                                   cagcagTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTTTAAAtaattta       <  1:2042120/91‑3 (MQ=25)

                                                                                   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)

                                                                                      cagTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATgg    <  1:1318909/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:2196955/91‑1 (MQ=255)

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

AATAGTATAAGATTGAAAAAGTTGTTACTGTTTTAAATGATCACGATGAAGTCATTCAATAAGAATGATTATGAAAATAGAAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATCATCTCACTGCTGTTTTTTAAAGGTTTATACCTCATCCTCTAAATTATTTAAAAATAATTAATGGTA  >  USA300TCH1516_ALE/96807‑96985

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