BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,113,299	G→T	6.6%	D38E (GAC→GAA)	USA300HOU_RS10725 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,113,299	0	G	T	6.6%	697.9 / 6.7	230	D38E (GAC→GAA)	USA300HOU_RS10725	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (115/95); new base T (17/3); total (132/98)
Fisher's exact test for biased strand distribution p-value = 8.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 238 total aligned reads displayed.

CTTGTAGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGCT  >  USA300TCH1516_ALE/2113159‑2113439

cTTGTAGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGGAATTA‑TTGTTTTTTCg                                                                                                                                              >  1:940478/1‑141 (MQ=255)

  tGTAGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTc                                                                                                                                            >  1:1052008/1‑141 (MQ=255)

   gTAGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAACTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGAAATTA‑TTGTTTTTTCGTCa                                                                                                                                           >  2:434734/1‑141 (MQ=255)

   gTAGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCAAACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCa                                                                                                                                           >  1:367549/1‑141 (MQ=255)

    tcggATGTTAGTCAGTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCAt                                                                                                                                          <  1:175870/139‑1 (MQ=255)

    tatgATGTTAGTTATGTTCCTAGTTCTCCTGATCAATTCCCTCCAAGTCACACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACAGGAATACCTTTAACTCCTGGAATTCTGTAATTA‑TTGTTTTTTCGTCAt                                                                                                                                          <  2:692242/138‑1 (MQ=255)

    tacgATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCACCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCAt                                                                                                                                          <  2:987162/138‑1 (MQ=255)

      ggATGTTGGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATgg                                                                                                                                        <  2:308161/141‑1 (MQ=255)

           ttatctaTTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGAt                                                                                                                                   <  2:379366/135‑1 (MQ=255)

             aGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGt                                                                                                                                 >  1:966268/1‑141 (MQ=255)

              gTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTa                                                                                                                                >  2:797878/1‑141 (MQ=255)

                cagtttCCTGTTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGACTAACTTTAAAACATGATACGGGGGCATTA‑TTGTTTTTTCGTCATGGTCGATGTaaa                                                                                                                              >  2:499544/4‑141 (MQ=255)

                cATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTaaa                                                                                                                              <  2:492377/141‑1 (MQ=255)

                cATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAAAGCTACAAAGACTTTTATATATTCCTACCATTGACACGAATACCTTAAAATCATGTAATTCTGTTATTA‑TTGTTTTTTCGTCATGGTCGGGGTaaa                                                                                                                              >  1:808752/1‑141 (MQ=255)

                  ttttCCTTGCCCTCCCCATCAATTCCAGCCAATTCCTGAATAACGCTCTCCATCCTTTTATCTTTTAGATCTTCTGATACGAATACCTTTAAATCATTGAATTCTGTAATTA‑TTTTTTTTTCGTCATGGTCGATGTAAAtt                                                                                                                            <  2:299395/141‑1 (MQ=255)

                      ccTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTGCATAATTTTATCTTTTAAATCTTCTGACACGCATACCTTTAAATCATGGAATTCTGTATTTA‑TTGTTTTTTGGTCTTGGTTGATGTAAATTTCTa                                                                                                                        >  2:155239/1‑141 (MQ=255)

                      ccTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTACATTTCTa                                                                                                                        >  2:973075/1‑141 (MQ=255)

                      ccTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTa                                                                                                                        >  1:68963/1‑141 (MQ=255)

                      ccTTGTCCTCCTCATAAATTCCAGCAAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTa                                                                                                                        >  1:553786/1‑141 (MQ=255)

                       cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAAGGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATGCCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAGATTTCtgg                                                                                                                       >  2:598823/1‑139 (MQ=255)

                       cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTACTTA‑TTGTTTTTTCGTCATGGTCGATGTGAATCTATAg                                                                                                                       >  1:135740/1‑141 (MQ=255)

                       cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAg                                                                                                                       >  1:407224/1‑141 (MQ=255)

                       cTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAg                                                                                                                       <  1:623988/141‑1 (MQ=255)

                           tcctccTCATCTATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTa                                                                                                                   <  2:1024632/141‑1 (MQ=255)

                            cctcctCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTAc                                                                                                                  <  1:895149/141‑1 (MQ=255)

                                   atcaATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGc                                                                                                           >  1:435802/1‑141 (MQ=255)

                                    tcaATTCCAGCCAATTCATACATAACGCTCTCCATCCTTTTATCTTTTAACTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCCGCt                                                                                                          >  2:617754/1‑141 (MQ=255)

                                      aaTTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCCGCTaa                                                                                                        >  2:531362/1‑141 (MQ=255)

                                      aaTTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTaa                                                                                                        >  1:344442/1‑141 (MQ=255)

                                        ttCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCTTCTTTTTTGATGTAAATTTCTAGCGTACCATCAGCTAAAt                                                                                                      >  2:955227/1‑141 (MQ=255)

                                           cAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaa                                                                                                   >  1:327945/1‑141 (MQ=255)

                                           cAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaa                                                                                                   <  2:1095287/141‑1 (MQ=255)

                                            cgCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaaa                                                                                                  <  1:118236/140‑1 (MQ=255)

                                            aGCCAATTCATACATAACGATCTCCATTGTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaaa                                                                                                  <  1:197728/141‑1 (MQ=255)

                                             gCCAATTCATACATCACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTaaaa                                                                                                 <  1:906568/141‑1 (MQ=255)

                                             gCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCGACTAAATTaaaa                                                                                                 >  1:653207/1‑141 (MQ=255)

                                                    catacataACGCTCTCCATACTTTTATCTTTTAAATCTACTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCCGCTACATTAAAAATAGCtt                                                                                          <  1:719574/141‑1 (MQ=255)

                                                       acataACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTAt                                                                                       >  1:87561/1‑141 (MQ=255)

                                                                tctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAACGTACCATAAGCTAAATTAAAAATGGCTTTATCTTCAtttt                                                                              >  2:883980/1‑141 (MQ=255)

                                                                 ctcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAACTCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTc                                                                             <  2:130983/141‑1 (MQ=255)

                                                                  tcCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCt                                                                            >  1:482818/1‑141 (MQ=255)

                                                                   ccATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTTTTTTTTTTTGATGTAAATTTCTGGCGTACCGTCAGCTAAATTAAAATTAGCTTTATTTTAATTTTCtt                                                                           >  2:509553/1‑141 (MQ=255)

                                                                    cATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCTTTTTTTTTGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAATTGGCTTCATCTTCATTTTCttt                                                                          >  2:944165/1‑141 (MQ=255)

                                                                    cATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCttt                                                                          >  1:931966/1‑141 (MQ=255)

                                                                    cATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCttt                                                                          >  1:464417/1‑141 (MQ=255)

                                                                     aTACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCtttt                                                                         >  1:865033/1‑140 (MQ=255)

                                                                      tACTTTTATCTTTTAAATATTCTGAAACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTAGTCATGTTCGCTGTAAATTTATAGCGAACCATCAACTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGt                                                                        >  2:41683/1‑141 (MQ=255)

                                                                          tttATCTTTTAAATAGTCTGCATCGAATGCGTTTAAATTATGGAATTCTGTACTTA‑TTTTTTTTTCGTCATGGTGGATGTCAATTTCTAGCGTACCATCATCTAAATTAAAAATAGCTTTATCTTCCTTTTCTTTGTATcc                                                                    <  2:867454/141‑1 (MQ=255)

                                                                                  ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCTTTGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTCCATTTTCTTTGTATCCTTTTTCtt                                                            >  2:903223/1‑141 (MQ=255)

                                                                                  ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtt                                                            >  1:642370/1‑141 (MQ=255)

                                                                                  ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtt                                                            >  1:292469/1‑141 (MQ=255)

                                                                                  ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTGATCTTCCTTTTCTTTGTATCCTTTTTCtt                                                            >  1:385910/1‑141 (MQ=255)

                                                                                  ttAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGGCCTGGTCGATGTAAATTTCTCGCGTACCATCAGCTACATTAAAAATAGCTTTATCTTCATTTTCTTTGTTTCCGTTTTCtt                                                            >  1:58306/1‑141 (MQ=255)

                                                                                   tAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCttt                                                           <  1:880320/141‑1 (MQ=255)

                                                                                   tAAATCTTCTGACACGAATACCTTTAAATCAGGGAATTCT‑AATTCACCTACTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCttt                                                           <  1:520905/141‑1 (MQ=255)

                                                                                    aaaTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtttt                                                          >  1:299691/1‑141 (MQ=255)

                                                                                    aaaTCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGAACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCtttt                                                          <  1:664744/141‑1 (MQ=255)

                                                                                         ttctGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAAttt                                                     >  1:757999/1‑141 (MQ=255)

                                                                                          tctGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGCTGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCGTTTTCTTTGTATCCTTTTTCTTTTTAtttt                                                    >  2:828753/1‑140 (MQ=255)

                                                                                          tctGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTTAtttt                                                    >  2:551235/1‑140 (MQ=255)

                                                                                             gACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTTTTTTTTTTCTATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATttcttc                                                 >  2:281210/1‑141 (MQ=255)

                                                                                             gACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATttcttc                                                 >  1:785700/1‑141 (MQ=255)

                                                                                              acacGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTACTTTCTTCa                                                >  1:83238/1‑141 (MQ=255)

                                                                                              acacGAATACCTTTAAATCATGGAATTCGGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCa                                                >  1:578297/1‑141 (MQ=255)

                                                                                               cacGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAGTTAAAAATAGCTTTATCTTCATTTTCTTTGTTTCCTTTTTCTTTTATTTTGTTCaa                                               >  2:107887/1‑141 (MQ=255)

                                                                                                acGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTTTTTTTTTTTTATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGGTTTTTCTTCTTTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCatt                                              >  2:56532/1‑139 (MQ=255)

                                                                                                 cGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTATTTTCTTCCAAg                                             >  2:34281/1‑141 (MQ=255)

                                                                                                 cGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTAGCCTTTTTCTTTTAATTGCTTCaaaa                                             >  1:499792/1‑140 (MQ=255)

                                                                                                  gAATACCTTTAAATCATGGAGTTCTGTAATTA‑TTTTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGa                                            <  2:725770/141‑1 (MQ=255)

                                                                                                  gAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGa                                            <  2:847997/141‑1 (MQ=255)

                                                                                                       ccTTTAAATCAGGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAGATTTCTAGCGTACCATCAGCTGAATTAAAAGTAGCTTTATCATCATGTTCTTTGTGTCCTTTTTCTTTTAGTTTCTTCACAGATGCAt                                       >  2:24471/1‑141 (MQ=255)

                                                                                                         tttAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTa                                     >  2:1040480/1‑141 (MQ=255)

                                                                                                          ttAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTaa                                    >  2:366349/1‑141 (MQ=255)

                                                                                                           tAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTCCCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTCAATTTCTTCAAAGATTCATTAat                                   <  2:792079/141‑1 (MQ=255)

                                                                                                             aaTCATGGAATTCTGTAATTA‑TTGTTTTTTTTTTTTTTTCGATGTAAATTTCTAGCTTACCATCAGCTAAATTAAAATTAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGTTTCATTAattt                                 >  2:471590/1‑141 (MQ=255)

                                                                                                             aaTCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattt                                 <  2:456933/141‑1 (MQ=255)

                                                                                                              aTCATGCAATTCTGTAATTT‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAatttt                                <  2:149945/141‑1 (MQ=255)

                                                                                                                  tGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCATAGATTCATTAattttcatt                            >  2:620202/1‑141 (MQ=255)

                                                                                                                  tGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattttcatt                            >  1:142988/1‑141 (MQ=255)

                                                                                                                  tGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAAATTCATTAattttcatt                            >  2:75854/1‑141 (MQ=255)

                                                                                                                   ggAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAattttcattt                           >  2:638233/1‑141 (MQ=255)

                                                                                                                      aTTCTGTAATTA‑TTGTTTTTTTTTTTTTTTTGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAATGATTTATTTatttttattttct                        >  2:677443/1‑140 (MQ=255)

                                                                                                                       ttCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCt                       <  1:609144/141‑1 (MQ=255)

                                                                                                                        tCTTTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTg                      <  1:323777/141‑1 (MQ=255)

                                                                                                                        tCTGTAATTC‑TTGTTTTTTCGCCATGGGCGATGTAAAGTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTg                      <  1:310027/141‑1 (MQ=255)

                                                                                                                        tCTGTAATTA‑TTTGTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTg                      <  2:20146/141‑1 (MQ=255)

                                                                                                                         ttGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCCTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGt                     <  2:185155/140‑1 (MQ=255)

                                                                                                                            tAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAACATCGCTTTATCTACTTTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTATTTTTCATTTTCTTGtttc                  >  2:393371/1‑139 (MQ=255)

                                                                                                                              atta‑ttGTTTTTTCGTCATTGTCGATGTAAATTTCTAGCGTACCTTCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATTCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTtcctc                <  2:505045/141‑1 (MQ=255)

                                                                                                                                ta‑ttGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTcctcct              <  1:881674/141‑1 (MQ=255)

                                                                                                                                    tGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCat           <  2:872608/141‑1 (MQ=255)

                                                                                                                                    tGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCat           <  1:972846/141‑1 (MQ=255)

                                                                                                                                      ttttttCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCatat         <  1:1093875/141‑1 (MQ=255)

                                                                                                                                       tttttCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTGTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATAtt        <  2:667102/141‑1 (MQ=255)

                                                                                                                                       ttgttCGTCATGGTCTATGTAAATTTCTAGCGTCCCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATAtt        <  2:702485/138‑1 (MQ=255)

                                                                                                                                           tCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAg    <  2:746333/141‑1 (MQ=255)

                                                                                                                                            cGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGc   >  1:487486/1‑141 (MQ=255)

                                                                                                                                             gTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTccatct  >  1:236866/1‑136 (MQ=255)

                                                                                                                                             gTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGct  <  1:426744/141‑1 (MQ=255)

                                                                                                                                             gTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCAGTTCCTCCTCATATTCAAGct  <  2:940478/141‑1 (MQ=255)

CTTGTAGGATGTTAGTCATTTTCCTTGTCCTCCTCATCAATTCCAGCCAATTCATACATAACGCTCTCCATACTTTTATCTTTTAAATCTTCTGACACGAATACCTTTAAATCATGGAATTCTGTAATTA‑TTGTTTTTTCGTCATGGTCGATGTAAATTTCTAGCGTACCATCAGCTAAATTAAAAATAGCTTTATCTTCATTTTCTTTGTATCCTTTTTCTTTTAATTTCTTCAAAGATTCATTAATTTTCATTTTCCTGTTCCTCCTCATATTCAAGCT  >  USA300TCH1516_ALE/2113159‑2113439

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