BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,885,369	T→C	5.3%	intergenic (‑151/+175)	ebh_2 ← / ← moeZ_2	Extracellular matrix‑binding protein ebh/putative adenylyltransferase/sulfurtransferase MoeZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,885,369	0	T	C	5.3%	712.5 / 20.6	224	intergenic (‑151/+175)	ebh_2/moeZ_2	Extracellular matrix‑binding protein ebh/putative adenylyltransferase/sulfurtransferase MoeZ
Reads supporting (aligned to +/- strand): ref base T (100/112); new base C (5/7); total (105/119)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 256 total aligned reads displayed.

ACGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGT  >  USA300TCH1516_ALE/1885279‑1885459

aCGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGt                                                                                            <  2:2882137/91‑1 (MQ=255)

aCGCTCCCTCATATTATATTTATAGCTATATTATTTTATACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGt                                                                                            >  1:3266549/1‑91 (MQ=255)

 cGCTCCCTCATATTATGTTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGtt                                                                                           <  2:3705539/91‑1 (MQ=255)

 cGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGtt                                                                                           <  1:1577084/91‑1 (MQ=255)

  gCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttt                                                                                          >  2:2492226/1‑91 (MQ=255)

  ccTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttt                                                                                          >  2:86430/2‑91 (MQ=255)

    tCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACACAATTTTATCATCATCATTTTATAAAATAAAGtttta                                                                                        <  1:3286316/91‑1 (MQ=255)

       cTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttttattt                                                                                     >  1:2332062/1‑91 (MQ=255)

        tCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAActttatttt                                                                                    <  1:1897943/91‑1 (MQ=255)

         cATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAActttatttta                                                                                   <  2:3773237/91‑1 (MQ=255)

         cATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAActttatttta                                                                                   >  2:2618381/1‑91 (MQ=255)

           tattatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTttattttatt                                                                                 >  1:543888/1‑91 (MQ=255)

           tattatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAACtttattttata                                                                                 <  2:3028576/91‑2 (MQ=255)

           tattatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAACtttattttata                                                                                 <  1:981451/91‑2 (MQ=255)

           tattatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAACtttattttata                                                                                 <  1:1437292/91‑2 (MQ=255)

             ttatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTAc                                                                               <  2:1407350/91‑1 (MQ=255)

             ttatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTAc                                                                               >  2:914154/1‑91 (MQ=255)

             ttatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTAc                                                                               <  1:1018973/91‑1 (MQ=255)

              tatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACa                                                                              <  2:2269940/91‑1 (MQ=255)

              tatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAACtttattttataaaa                                                                              <  1:3626948/91‑5 (MQ=38)

               atattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACaa                                                                             >  1:2740943/1‑91 (MQ=255)

                tattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAAt                                                                            <  1:1278512/91‑1 (MQ=255)

                 attTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAAtt                                                                           >  2:2788949/1‑91 (MQ=255)

                   ttATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTa                                                                         <  1:2559804/91‑1 (MQ=255)

                     ataGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGt                                                                       <  1:1307348/91‑1 (MQ=255)

                     ataGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGt                                                                       >  1:170933/1‑91 (MQ=255)

                          tataTTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttttatttta                                                                                   <  1:3527895/74‑1 (MQ=255)

                          tataTTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCg                                                                  <  2:1137969/91‑1 (MQ=255)

                           atatTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGt                                                                 >  1:2519920/1‑91 (MQ=255)

                           atatTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGt                                                                 >  1:2359305/1‑91 (MQ=255)

                           atatTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGt                                                                 >  2:770542/1‑91 (MQ=255)

                            tattatTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTa                                                                         >  2:1717658/1‑82 (MQ=255)

                             attattTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGttt                                                               <  2:2386718/91‑1 (MQ=255)

                              ttattTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGtttt                                                              >  2:3328004/1‑91 (MQ=255)

                               tattTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTg                                                             >  2:1049704/1‑91 (MQ=255)

                                attTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGt                                                            <  1:1495756/91‑1 (MQ=255)

                                  tttGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTaa                                                          >  2:3525306/1‑91 (MQ=255)

                                   ttGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAAt                                                         <  2:1882482/91‑1 (MQ=255)

                                     gTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTa                                                       <  2:2359305/91‑1 (MQ=255)

                                       aCAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTaaa                                                     >  2:2265266/1‑91 (MQ=255)

                                       aCAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTaaa                                                     >  2:3811582/1‑91 (MQ=255)

                                        cAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTaaaa                                                    <  2:2888316/91‑1 (MQ=255)

                                         aaTATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAAt                                                   >  2:3285968/1‑91 (MQ=255)

                                          atatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATc                                                  <  1:1759365/91‑1 (MQ=255)

                                           tatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCt                                                 >  2:3574957/1‑91 (MQ=255)

                                           tatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCt                                                 >  1:167197/1‑91 (MQ=255)

                                            atTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCtt                                                >  2:543495/1‑91 (MQ=255)

                                               ttAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTtaat                                             >  2:2180192/1‑91 (MQ=255)

                                               ttAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTtaat                                             <  2:540438/91‑1 (MQ=255)

                                                tAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTtaata                                            <  1:320382/91‑1 (MQ=255)

                                                tAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTtaata                                            <  2:2312426/91‑1 (MQ=255)

                                                  aGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATaaa                                          >  2:1573622/1‑91 (MQ=255)

                                                  aGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATaaa                                          >  1:763100/1‑91 (MQ=255)

                                                    ttCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAAc                                        >  2:1330027/1‑91 (MQ=255)

                                                    ttCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAAc                                        >  2:2241683/1‑91 (MQ=255)

                                                     tCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACa                                       <  1:3404181/91‑1 (MQ=255)

                                                     tCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACa                                       <  2:1205196/91‑1 (MQ=255)

                                                     tCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACa                                       <  1:2654690/91‑1 (MQ=255)

                                                       aCTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATa                                     <  2:2257974/91‑1 (MQ=255)

                                                          cAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACa                                  <  2:3105629/91‑1 (MQ=255)

                                                          cAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACa                                  <  2:2978366/91‑1 (MQ=255)

                                                           aaTTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACaa                                 <  1:700946/91‑1 (MQ=255)

                                                           aaTTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACaa                                 <  1:1770024/91‑1 (MQ=255)

                                                            aTTTTATCATCATCATTTTATAAACTAAAGTTTTATTTTATTACAATTTATTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAg                                <  2:3597104/91‑1 (MQ=255)

                                                            aTTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAg                                >  2:3689888/1‑91 (MQ=255)

                                                            aTTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAg                                <  2:2613543/91‑1 (MQ=255)

                                                            aTTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAg                                <  2:3722173/91‑1 (MQ=255)

                                                              tttATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAATACATAACAAGaa                              <  1:2968802/91‑1 (MQ=255)

                                                                tATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACg                            >  2:3461934/1‑91 (MQ=255)

                                                                 atcatcatcATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGc                           >  2:2954929/1‑91 (MQ=255)

                                                                  tcatcatcaTTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCt                          >  1:1082750/1‑91 (MQ=255)

                                                                    atcatcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCttt                        >  2:2470351/1‑91 (MQ=255)

                                                                    atcatcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCttt                        >  1:2811618/1‑91 (MQ=255)

                                                                     tcatcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTa                       <  1:3751030/91‑1 (MQ=255)

                                                                      catcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAAGGCTTTAg                      >  1:1945509/1‑91 (MQ=255)

                                                                      catcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAg                      <  1:2459977/91‑1 (MQ=255)

                                                                       atcttTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGa                     <  2:3626964/91‑1 (MQ=255)

                                                                        tcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGAt                    <  1:508584/91‑1 (MQ=255)

                                                                         catTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGAtt                   <  1:1211958/91‑1 (MQ=255)

                                                                           ttttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGt                 <  2:2044645/91‑1 (MQ=255)

                                                                           ttttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGt                 <  2:593857/91‑1 (MQ=255)

                                                                           ttttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGt                 <  2:485588/91‑1 (MQ=255)

                                                                            gtaataaaataaaACTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGtt                >  1:1622928/4‑91 (MQ=39)

                                                                             ttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt               >  2:252422/1‑91 (MQ=255)

                                                                             ttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt               <  1:3811364/91‑1 (MQ=255)

                                                                               ataaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTg             <  1:2901371/91‑1 (MQ=255)

                                                                               ataaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTg             >  1:852680/1‑91 (MQ=255)

                                                                               ataaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTg             >  2:1288898/1‑91 (MQ=255)

                                                                                taaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGt            >  2:1756066/1‑91 (MQ=255)

                                                                                 aaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTa           >  2:2711803/1‑91 (MQ=255)

                                                                                 aaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTa           <  1:1139575/91‑1 (MQ=255)

                                                                                 aaaataaaGTTTTATTTTATTACAATTTAATTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTa           <  2:1466560/91‑1 (MQ=255)

                                                                                    ataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACa        <  1:1330027/91‑1 (MQ=255)

                                                                                     taaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAg       <  2:1926594/91‑1 (MQ=255)

                                                                                      aaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATGACAAGAACGCTTTAGATTGTTTTGTAACAGa      <  2:41074/91‑1 (MQ=255)

                                                                                       aaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGaa     >  2:3520648/1‑91 (MQ=255)

                                                                                       aaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGaa     >  1:2490442/1‑91 (MQ=255)

                                                                                       aaGTATTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACGGaa     >  1:3064075/1‑91 (MQ=255)

                                                                                          ttttattttaTTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGt  <  2:3800440/91‑1 (MQ=255)

                                                                                          ttttattttaTTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGt  <  1:2064525/91‑1 (MQ=255)

ACGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGT  >  USA300TCH1516_ALE/1885279‑1885459

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: `‑`