BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,861,759	C→T	10.8%	A481V (GCA→GTA)	bceB_3 →	Bacitracin export permease protein BceB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,861,759	0	C	T	10.8%	353.7 / 24.4	120	A481V (GCA→GTA)	bceB_3	Bacitracin export permease protein BceB
Reads supporting (aligned to +/- strand): ref base C (55/52); new base T (10/3); total (65/55)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 122 total aligned reads displayed.

AAAAGTATCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATATTTGTTAC  >  USA300TCH1516_ALE/2861669‑2861849

   aGTATCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAAt                                                                                         <  1:584232/91‑1 (MQ=255)

    gTATCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGTAAtt                                                                                        >  1:521287/1‑91 (MQ=255)

     tATCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTg                                                                                       >  1:1190854/1‑91 (MQ=255)

      aTCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGc                                                                                      <  1:1443679/91‑1 (MQ=255)

      aTCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGc                                                                                      <  1:348675/91‑1 (MQ=255)

       tCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACCAAAGGAAAAAGCCCAAGCAATTGCa                                                                                     >  1:340836/1‑91 (MQ=255)

       tCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCa                                                                                     >  1:1355361/1‑91 (MQ=255)

       tCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCa                                                                                     >  1:341673/1‑91 (MQ=255)

        cagcagTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGTAATTGCaa                                                                                    >  1:1695786/1‑91 (MQ=255)

        cagcagTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCaa                                                                                    >  1:1038177/1‑91 (MQ=255)

        cagcagTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCaa                                                                                    <  1:3939684/91‑1 (MQ=255)

         agcagTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCaaa                                                                                   >  1:2412801/1‑91 (MQ=255)

           cagTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAAc                                                                                 <  1:2281351/91‑1 (MQ=255)

             gTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAg                                                                               >  1:3255165/1‑91 (MQ=255)

              ttGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGt                                                                              <  1:1308424/91‑1 (MQ=255)

              ttGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGt                                                                              >  1:3152720/1‑91 (MQ=255)

              ttGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGt                                                                              >  1:2138847/1‑91 (MQ=255)

               tGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGtt                                                                             <  1:596194/91‑1 (MQ=255)

               tGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGtt                                                                             >  1:1028688/1‑91 (MQ=255)

                gAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGttt                                                                            >  1:1422175/1‑91 (MQ=255)

                 aaaaCACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTg                                                                           <  1:281415/91‑1 (MQ=255)

                  aaaCACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGTAATTGCAAAACAGTTTgg                                                                          >  1:282585/1‑91 (MQ=255)

                  aaaCACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGTAATTGCAAAACAGTTTgg                                                                          >  1:282555/1‑91 (MQ=255)

                  aaaCACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTgg                                                                          >  1:514242/1‑91 (MQ=255)

                   aaCACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGga                                                                         <  1:3879917/91‑1 (MQ=255)

                   aaCACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGga                                                                         >  1:1837871/1‑91 (MQ=255)

                   aaCACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGga                                                                         >  1:1945957/1‑91 (MQ=255)

                      acaAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGAc                                                                      >  1:3040543/1‑91 (MQ=255)

                        aaaaTAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACaa                                                                    <  1:3524316/91‑1 (MQ=255)

                            tAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGtc                                                                <  1:3320924/91‑1 (MQ=255)

                             aaCGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGtca                                                               <  1:2489542/91‑1 (MQ=255)

                              aCGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGtcat                                                              >  1:1959431/1‑91 (MQ=255)

                                gTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGtcatca                                                            >  1:2889096/1‑91 (MQ=255)

                                 ttCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCAc                                                           >  1:1253099/1‑91 (MQ=255)

                                 ttCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCAc                                                           <  1:2825497/91‑1 (MQ=255)

                                  tCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACg                                                          <  1:2474157/91‑1 (MQ=255)

                                  tCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACg                                                          >  1:42892/1‑91 (MQ=255)

                                           ttATGGATATGATATTATACAAACATCACAAAGGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaa                                                 <  1:3544603/91‑1 (MQ=255)

                                           ttATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaa                                                 <  1:2515282/91‑1 (MQ=255)

                                           ttATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaa                                                 <  1:3679178/91‑1 (MQ=255)

                                            tATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaaa                                                <  1:3707881/91‑1 (MQ=255)

                                             aTGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaaat                                               >  1:2268270/1‑91 (MQ=255)

                                               ggATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaaatga                                             >  1:2394370/1‑91 (MQ=255)

                                               ggATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaaatga                                             >  1:1727015/1‑91 (MQ=255)

                                                gatatgatatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGatgaaatgaa                                            <  1:3057094/91‑1 (MQ=255)

                                                 atatgatatTAAACAAACATCACAAAAGGAAAAAGCTCAAGTAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGAtgaaatgaaa                                           >  1:3125560/1‑91 (MQ=255)

                                                 atatgatatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGAtgaaatgaaa                                           >  1:600732/1‑91 (MQ=255)

                                                  tatgatatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGaaaa                                          <  1:1214495/91‑1 (MQ=255)

                                                  tatgatatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGaaaa                                          <  1:1919589/91‑1 (MQ=255)

                                                  tatgatatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGaaaa                                          <  1:1214084/91‑1 (MQ=255)

                                                    tgatatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGaaaaaa                                        <  1:66009/91‑1 (MQ=255)

                                                      atatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAaga                                      <  1:2713790/91‑1 (MQ=255)

                                                      atatTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAaga                                      <  1:1641010/91‑1 (MQ=255)

                                                        atTAAGCAAACATCACAAAAGGAAATAGCTCAAGCAATTGCAAAACAGTTTGGGGACAAAGTCATCACGTATGATGAAATGAAAAAAGAgg                                    <  1:3860122/91‑1 (MQ=255)

                                                        atTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACTAAGTCATCACGTATGATGAAATGAAAAAAGAgg                                    >  1:3392714/1‑91 (MQ=255)

                                                        atTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAgg                                    <  1:1780872/91‑1 (MQ=255)

                                                        atTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAgg                                    <  1:334890/91‑1 (MQ=255)

                                                        atTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAgg                                    <  1:2133327/91‑1 (MQ=255)

                                                        atTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAgg                                    <  1:2564027/91‑1 (MQ=255)

                                                         ttAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGt                                   <  1:263460/91‑1 (MQ=255)

                                                          tAAACAAACATCACAAAAGGAAAAAGCTCAAGTAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTc                                  >  1:702505/1‑91 (MQ=255)

                                                                   aTCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGctact                         <  1:1328397/91‑1 (MQ=255)

                                                                    tCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGctacta                        >  1:2757843/1‑91 (MQ=255)

                                                                     cacaAAAGGAAAAAGCTCAAGTAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTaa                       >  1:3532/1‑91 (MQ=255)

                                                                     cacaAAAGGAAAAAGCTCAAGTAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTaa                       >  1:3938381/1‑91 (MQ=255)

                                                                     cacaAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTaa                       <  1:3445450/91‑1 (MQ=255)

                                                                     cacaAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTaa                       >  1:1101903/1‑91 (MQ=255)

                                                                     cacaAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTaa                       >  1:1102248/1‑91 (MQ=255)

                                                                     cacaAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTaa                       >  1:3750282/1‑91 (MQ=255)

                                                                     cacaAAAGGAAAAAGCTCAAGCAATGGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTaa                       <  1:3445417/91‑1 (MQ=255)

                                                                      acaAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAAc                      >  1:1167695/1‑91 (MQ=255)

                                                                       caAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACg                     >  1:1166050/1‑91 (MQ=255)

                                                                       caAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACg                     <  1:3131381/91‑1 (MQ=255)

                                                                       caAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACg                     >  1:3651292/1‑91 (MQ=255)

                                                                       caAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACg                     <  1:2259322/91‑1 (MQ=255)

                                                                         aaaGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGt                   >  1:1927602/1‑91 (MQ=255)

                                                                          aaGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTa                  >  1:2952943/1‑91 (MQ=255)

                                                                          aaGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTa                  >  1:2952898/1‑91 (MQ=255)

                                                                           aGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTAt                 <  1:3003310/91‑1 (MQ=255)

                                                                            ggAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTAtt                >  1:48055/1‑91 (MQ=255)

                                                                            ggAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTAtt                >  1:2788224/1‑91 (MQ=255)

                                                                             gAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTc               >  1:3552515/1‑91 (MQ=255)

                                                                             gAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTc               >  1:676129/1‑91 (MQ=255)

                                                                              aaaaaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCt              <  1:4127468/91‑1 (MQ=255)

                                                                              aaaaaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCt              <  1:1325871/91‑1 (MQ=255)

                                                                              aaaaaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCt              <  1:4096052/91‑1 (MQ=255)

                                                                               aaaaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTa             >  1:592429/1‑91 (MQ=255)

                                                                               aaaaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTa             >  1:3280488/1‑91 (MQ=255)

                                                                                aaaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTaa            >  1:208867/1‑91 (MQ=255)

                                                                                aaaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTaa            <  1:2636454/91‑1 (MQ=255)

                                                                                 aaGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAat           >  1:1156157/1‑91 (MQ=255)

                                                                                   gCTCAAGTAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAatat         >  1:1382582/1‑91 (MQ=255)

                                                                                   gCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAatat         >  1:778508/1‑91 (MQ=255)

                                                                                     tCAAGTAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATAttt       <  1:3317846/91‑1 (MQ=255)

                                                                                     ccAAGTACTTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATAttt       <  1:1393624/90‑1 (MQ=255)

                                                                                      cAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATATTTg      <  1:834301/91‑1 (MQ=255)

                                                                                       aaGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATATTTGt     <  1:3926105/91‑1 (MQ=255)

                                                                                        aGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATATTTGtt    <  1:3834662/91‑1 (MQ=255)

                                                                                        aGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATATTTGtt    <  1:3834484/91‑1 (MQ=255)

                                                                                          cAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATATTTGTTAc  <  1:2201911/91‑1 (MQ=255)

AAAAGTATCAGCAGTTGAAAACACAAAATAACGTTCATACTTTTTATGGATATGATATTAAACAAACATCACAAAAGGAAAAAGCTCAAGCAATTGCAAAACAGTTTGGAGACAAAGTCATCACGTATGATGAAATGAAAAAAGAGGTCGATGCTACTAACGGTATTCTAATATTTGTTAC  >  USA300TCH1516_ALE/2861669‑2861849

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