BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,430,691:1	+A	9.7%	intergenic (‑115/‑270)	suhB_2 ← / → birA_2	Inositol‑1‑monophosphatase/Bifunctional ligase/repressor BirA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,430,691	1	.	A	9.7%	421.1 / 41.8	113	intergenic (‑115/‑270)	suhB_2/birA_2	Inositol‑1‑monophosphatase/Bifunctional ligase/repressor BirA
Reads supporting (aligned to +/- strand): ref base . (47/55); new base A (4/7); total (51/62)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 124 total aligned reads displayed.

TTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATA  >  USA300TCH1516_ALE/2430600‑2430778

ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTg                                                                                          <  1:2650078/91‑1 (MQ=255)

ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTg                                                                                          <  1:1626386/91‑1 (MQ=255)

ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTg                                                                                          <  1:2348879/91‑1 (MQ=255)

ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTg                                                                                          <  1:411031/91‑1 (MQ=255)

ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTg                                                                                          <  1:3059279/91‑1 (MQ=255)

ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTg                                                                                          <  1:1802829/91‑1 (MQ=255)

ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTg                                                                                          <  1:667150/91‑1 (MQ=255)

 tatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTga                                                                                         <  1:4093914/91‑2 (MQ=255)

    cattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTa                                                                                      >  1:1366682/1‑91 (MQ=255)

     attTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAg                                                                                     <  1:1339358/91‑1 (MQ=255)

     attTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAg                                                                                     <  1:1339359/91‑1 (MQ=255)

     attTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAg                                                                                     <  1:1339374/91‑1 (MQ=255)

     attTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAg                                                                                     <  1:282777/91‑1 (MQ=255)

       tttATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTGACTAGcc                                                                                   <  1:559151/91‑1 (MQ=255)

       tttATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTGACTAGcc                                                                                   <  1:2487807/91‑1 (MQ=255)

        ttattattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTGACTAGCCg                                                                                  <  1:2998982/91‑1 (MQ=255)

         tattattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTGACTAGCCGa                                                                                 <  1:2536071/91‑1 (MQ=255)

         tatgattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAG‑CACGTGACTAGCCGa                                                                                 <  1:3218653/91‑1 (MQ=255)

            tattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACaa                                                                              >  1:1964701/1‑91 (MQ=255)

               taCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATa                                                                           >  1:3021802/1‑91 (MQ=255)

               taCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATa                                                                           >  1:669951/1‑91 (MQ=255)

                aCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAg                                                                          <  1:2463302/91‑1 (MQ=255)

                aCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑‑AGAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAga                                                                         >  1:1736763/1‑91 (MQ=255)

                aCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑‑AAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAga                                                                         >  1:1736528/1‑91 (MQ=255)

                aCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑‑AAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAga                                                                         >  1:1736610/1‑91 (MQ=255)

                 cccTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAga                                                                         <  1:290144/91‑1 (MQ=255)

                   cTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAgaaa                                                                       <  1:718725/91‑1 (MQ=255)

                    tataTTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAgaaag                                                                      >  1:3628929/1‑91 (MQ=255)

                        ttAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAgaaaggaaa                                                                  <  1:2699974/91‑1 (MQ=255)

                        ttAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAgaaaggaaa                                                                  >  1:2973263/1‑91 (MQ=255)

                         tAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAgaaaggaaag                                                                 <  1:1141800/91‑1 (MQ=255)

                          aaaaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAaagt                                                                <  1:3107422/91‑1 (MQ=255)

                            aaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAaagtaa                                                              >  1:3539404/1‑91 (MQ=255)

                             aGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAaagtaag                                                             >  1:3695843/1‑91 (MQ=255)

                             aGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAaagtaag                                                             <  1:1491347/91‑1 (MQ=255)

                              gAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAaagtaagt                                                            <  1:1080082/91‑1 (MQ=255)

                              gAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAaagtaagt                                                            <  1:1229414/91‑1 (MQ=255)

                                aTCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAgtaagtaa                                                          <  1:2261170/91‑1 (MQ=255)

                                  caacaaTACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGtaata                                                        <  1:1407896/91‑1 (MQ=255)

                                  caacaaTACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGtaata                                                        <  1:2532055/91‑1 (MQ=255)

                                  caacaaTACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGtaata                                                        <  1:2532325/91‑1 (MQ=255)

                                      aaTACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTaataaata                                                    <  1:2883490/91‑1 (MQ=255)

                                        tACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATAtt                                                  <  1:3943685/91‑1 (MQ=255)

                                        tACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATAtt                                                  <  1:1168558/91‑1 (MQ=255)

                                        tACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATAtt                                                  <  1:1863285/91‑1 (MQ=255)

                                         aCAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTg                                                 >  1:2384861/1‑91 (MQ=255)

                                          cAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGa                                                <  1:3648264/91‑1 (MQ=255)

                                          cAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGa                                                >  1:4060319/1‑91 (MQ=255)

                                           aaCTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGaa                                               <  1:2647543/91‑1 (MQ=255)

                                           aaCTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGaa                                               <  1:2647761/91‑1 (MQ=255)

                                            aCTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAg                                              >  1:2852001/1‑91 (MQ=255)

                                                  gACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgttg                                        <  1:3809527/91‑1 (MQ=255)

                                                  gACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgttg                                        <  1:4065296/91‑1 (MQ=255)

                                                    cttcttCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATGTTGAAGATGTtgtg                                      >  1:1065013/1‑91 (MQ=255)

                                                    cttcttCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTtgtg                                      >  1:905404/1‑91 (MQ=255)

                                                    cttcttCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTtgtg                                      >  1:853919/1‑91 (MQ=255)

                                                     ttcttcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGa                                     <  1:1170758/91‑1 (MQ=255)

                                                     ttcttcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGa                                     <  1:1172587/91‑1 (MQ=255)

                                                      tcttcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGAt                                    >  1:410774/1‑91 (MQ=255)

                                                        ttcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGt                                  >  1:162613/1‑91 (MQ=255)

                                                        ttcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGt                                  >  1:1716033/1‑91 (MQ=255)

                                                         tcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTa                                 >  1:2458873/1‑91 (MQ=255)

                                                         tcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTa                                 >  1:2251870/1‑91 (MQ=255)

                                                         tcATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTa                                 >  1:4091879/1‑91 (MQ=255)

                                                           aTTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAAc                               <  1:3627194/91‑1 (MQ=255)

                                                            ttttATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACt                              <  1:2581715/91‑1 (MQ=255)

                                                               tATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGa                           <  1:2374785/91‑1 (MQ=255)

                                                               tATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGa                           <  1:734873/91‑1 (MQ=255)

                                                               tATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGa                           <  1:543370/91‑1 (MQ=255)

                                                                aTGCATAAAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGa                           >  1:932371/1‑91 (MQ=255)

                                                                aTGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGaa                          >  1:3436734/1‑91 (MQ=255)

                                                                aTGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGaa                          >  1:3436603/1‑91 (MQ=255)

                                                                 tGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGaat                         >  1:838750/1‑90 (MQ=255)

                                                                 tGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                         >  1:2752469/1‑91 (MQ=255)

                                                                 tGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                         <  1:3884106/91‑1 (MQ=255)

                                                                 tGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                         <  1:365423/91‑1 (MQ=255)

                                                                 tGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                         >  1:2037071/1‑91 (MQ=255)

                                                                 tGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAAc                         <  1:212673/91‑1 (MQ=255)

                                                                  gCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGCACg                        >  1:3348483/1‑91 (MQ=255)

                                                                  gCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACg                        <  1:1260687/91‑1 (MQ=255)

                                                                   cAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGa                       <  1:1033364/91‑1 (MQ=255)

                                                                   cAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGa                       <  1:1451626/91‑1 (MQ=255)

                                                                             gtCGGCTAG‑CACGTGACTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCt              >  1:1852753/2‑91 (MQ=39)

                                                                               cGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTAt            >  1:1684271/1‑91 (MQ=255)

                                                                                ggCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATc           >  1:1435012/1‑91 (MQ=255)

                                                                                 gCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCt          >  1:2179820/1‑91 (MQ=255)

                                                                                 gCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCt          >  1:2179358/1‑91 (MQ=255)

                                                                                 gCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCt          >  1:3218708/1‑91 (MQ=255)

                                                                                 gCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCt          >  1:2179332/1‑91 (MQ=255)

                                                                                 gCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCt          >  1:2179104/1‑91 (MQ=255)

                                                                                  cTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTg         >  1:3968621/1‑91 (MQ=255)

                                                                                   tAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGt        <  1:3351502/91‑1 (MQ=255)

                                                                                    aGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGtt       <  1:3123732/91‑1 (MQ=255)

                                                                                     gTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTta      <  1:2446151/91‑1 (MQ=255)

                                                                                      gcACGTGACTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTta      >  1:3998576/2‑91 (MQ=255)

                                                                                      gcACGTGACTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTta      >  1:1990364/2‑91 (MQ=255)

                                                                                       cACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtata    <  1:2329768/91‑1 (MQ=255)

                                                                                         tgTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtatata  <  1:504750/90‑1 (MQ=255)

                                                                                         cGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtatata  <  1:229208/91‑1 (MQ=255)

                                                                                         cGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtatata  <  1:3472006/91‑1 (MQ=255)

TTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCAT‑AAAAAAATCGGCTAGTCACGTG‑CTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATA  >  USA300TCH1516_ALE/2430600‑2430778

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