BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	950,098	T→A	18.1%	intergenic (+92/‑270)	yugI_2 → / → USA300HOU_RS04740	General stress protein 13/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	950,098	0	T	A	18.1%	403.5 / 80.0	149	intergenic (+92/‑270)	yugI_2/USA300HOU_RS04740	General stress protein 13/putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (48/74); new base A (14/13); total (62/87)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 159 total aligned reads displayed.

AACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCT  >  USA300TCH1516_ALE/950010‑950188

aaCAAATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGAC‑A‑AAAAGAGATTAGTATCTAATAAAATTTAATAGa                                                                                        <  1:3632959/91‑1 (MQ=18)

 aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                                         <  1:234822/91‑1 (MQ=38)

 aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                                         <  1:1327351/91‑1 (MQ=38)

     aTAAATCGTACCGAAAATCATACAACGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATAc                                                                                     <  1:3859224/91‑1 (MQ=25)

     aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                     >  1:1456394/1‑91 (MQ=255)

     aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                     >  1:1456360/1‑91 (MQ=255)

     aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                     >  1:1123979/1‑91 (MQ=255)

       aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTa                                                                                   <  1:2810610/91‑1 (MQ=255)

        aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                                  <  1:1430157/91‑1 (MQ=255)

        aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                                  <  1:132517/91‑1 (MQ=255)

        aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                                  >  1:3674747/1‑91 (MQ=255)

         aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAt                                                                                 >  1:1039448/1‑91 (MQ=255)

          tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                                <  1:1546167/91‑1 (MQ=255)

           cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAAtct                                                                               <  1:1409552/91‑1 (MQ=25)

             tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                             <  1:2320094/91‑1 (MQ=255)

             tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                             >  1:3451376/1‑91 (MQ=255)

             tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                             <  1:3903179/91‑1 (MQ=255)

             tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAATctct                                                                             <  1:1925366/91‑1 (MQ=25)

              aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtt                                                                            >  1:1263977/1‑91 (MQ=255)

               ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                           >  1:3387485/1‑91 (MQ=255)

               ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                           >  1:1823617/1‑91 (MQ=255)

               ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAATCTCttt                                                                           <  1:1836118/91‑1 (MQ=21)

                cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                          <  1:3713297/91‑1 (MQ=255)

                cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                          <  1:1894636/91‑1 (MQ=255)

                cGAAAATCATACAAAGGGTCTGAAATGAAAGGTTCTTAGACTATGAAAGAGATTAGTGTCTATTAAATTTTATTAGATACTAATCTCtttt                                                                          <  1:3310122/91‑1 (MQ=255)

                 gAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGAC‑A‑AAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAATCTCttttata                                                                       <  1:3039284/91‑4 (MQ=14)

                     aTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCt                                                                     >  1:1088536/1‑91 (MQ=255)

                      tCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTa                                                                    >  1:3570507/1‑91 (MQ=255)

                       cATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTAc                                                                   <  1:479753/91‑1 (MQ=255)

                         tACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGa                                                                 >  1:1350559/1‑91 (MQ=255)

                          aCAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGAt                                                                >  1:3940286/1‑91 (MQ=255)

                          aCAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGAt                                                                <  1:2415672/91‑1 (MQ=255)

                             aaGGGTCTGAAGTGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAc                                                             <  1:3752489/91‑1 (MQ=255)

                             aaGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAc                                                             >  1:459778/1‑91 (MQ=255)

                               gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                                           <  1:1892388/91‑1 (MQ=255)

                                 gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                                         <  1:3131815/91‑1 (MQ=255)

                                  tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGCAat                                                        >  1:3341865/1‑91 (MQ=255)

                                   cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                       >  1:2845841/1‑91 (MQ=255)

                                   cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                       >  1:3871142/1‑91 (MQ=255)

                                   cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                       >  1:2656221/1‑91 (MQ=255)

                                    tgaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAatat                                                      <  1:2653472/91‑1 (MQ=255)

                                      aaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtga                                                    >  1:1020030/1‑91 (MQ=255)

                                       aatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                                   <  1:2816959/91‑1 (MQ=255)

                                       aatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                                   >  1:2090366/1‑91 (MQ=255)

                                        atgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgatt                                                  <  1:555638/91‑1 (MQ=255)

                                         tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                                 <  1:3942952/91‑1 (MQ=255)

                                         tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                                 >  1:2522517/1‑91 (MQ=255)

                                         tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                                 <  1:621606/91‑1 (MQ=255)

                                         tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTGTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                                 <  1:3302629/91‑1 (MQ=255)

                                           aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                               <  1:1256949/91‑1 (MQ=255)

                                           aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                               <  1:307382/91‑1 (MQ=255)

                                           aaaGTTTCTTAGACTATAAAAGAGATTAGTAGCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAGGTAATAtgatagat                                               >  1:3741939/1‑91 (MQ=255)

                                            aaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATgattgatt                                              <  1:590255/91‑1 (MQ=255)

                                             aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                             >  1:2716303/1‑91 (MQ=255)

                                             aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                             <  1:3026594/91‑1 (MQ=255)

                                                    tAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTac                                      <  1:1115607/91‑1 (MQ=255)

                                                       aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                                   <  1:535744/91‑1 (MQ=255)

                                                       aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                                   >  1:1965389/1‑91 (MQ=255)

                                                       aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                                   <  1:2857340/91‑1 (MQ=255)

                                                       aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                                   <  1:2518975/91‑1 (MQ=255)

                                                        cTATAAAAGTGATTAGTATCTATTAAATTTTATTAGATACTATTCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                                  <  1:2513383/91‑1 (MQ=255)

                                                          ataAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTAc                                <  1:1507494/91‑1 (MQ=255)

                                                            aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                              <  1:3430016/91‑1 (MQ=255)

                                                            aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                              <  1:723440/91‑1 (MQ=255)

                                                             aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                             <  1:1079546/91‑1 (MQ=255)

                                                             aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                             <  1:2378682/91‑1 (MQ=255)

                                                             aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                             <  1:2378418/91‑1 (MQ=255)

                                                                gagaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATgg                          <  1:1497065/91‑1 (MQ=255)

                                                                  gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt                        <  1:1508598/91‑1 (MQ=255)

                                                                    ttAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa                      <  1:1843535/91‑1 (MQ=255)

                                                                    ttAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa                      >  1:480276/1‑91 (MQ=25)

                                                                      aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg                    >  1:2621664/1‑91 (MQ=37)

                                                                      aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg                    >  1:462284/1‑91 (MQ=37)

                                                                      aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg                    >  1:3845382/1‑91 (MQ=37)

                                                                       gTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg                   >  1:2309723/1‑91 (MQ=37)

                                                                        tATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGt                  <  1:2393935/91‑1 (MQ=255)

                                                                         aTCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg                 <  1:4044651/91‑1 (MQ=255)

                                                                          tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa                <  1:598892/91‑1 (MQ=255)

                                                                          tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa                <  1:329629/91‑1 (MQ=255)

                                                                          tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa                <  1:1455320/91‑1 (MQ=255)

                                                                          tCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa                >  1:2568361/1‑91 (MQ=37)

                                                                          tCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa                >  1:2569683/1‑91 (MQ=37)

                                                                           cTATGAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGAc               >  1:3280748/1‑91 (MQ=255)

                                                                            tATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa              >  1:3389478/1‑91 (MQ=255)

                                                                            tATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa              >  1:1407739/1‑91 (MQ=255)

                                                                             aTTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACat             <  1:2894763/91‑1 (MQ=255)

                                                                             aTTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACat             >  1:1962470/1‑91 (MQ=255)

                                                                                 aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc         <  1:2713603/91‑1 (MQ=255)

                                                                                 aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc         <  1:2713659/91‑1 (MQ=255)

                                                                                 aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAAGATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc         <  1:4142480/91‑1 (MQ=255)

                                                                                   ttttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAt       <  1:2823676/91‑1 (MQ=255)

                                                                                    tttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAtt      <  1:555072/91‑1 (MQ=255)

                                                                                     ttattaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTa     >  1:3153805/1‑91 (MQ=255)

                                                                                       attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc   >  1:868016/1‑91 (MQ=255)

                                                                                       attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc   <  1:472097/91‑1 (MQ=255)

                                                                                       attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc   >  1:4020751/1‑91 (MQ=255)

                                                                                       attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc   >  1:3774276/1‑91 (MQ=255)

                                                                                       attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc   <  1:216493/91‑1 (MQ=255)

                                                                                        ttaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCt  >  1:1196559/1‑91 (MQ=255)

                                                                                        ttaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCt  <  1:1416989/91‑1 (MQ=255)

AACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCT  >  USA300TCH1516_ALE/950010‑950188

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