BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	110,341	A→T	D67V (GAT→GTT)	USA300HOU_RS00515 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	110,341	0	A	T	99.5%	1566.5 / ‑4.3	419	D67V (GAT→GTT)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base T (176/240); total (177/242)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 419 total aligned reads displayed.

GAAAAAACATTGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGATAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCTATATATGA  >  USA300TCH1516_ALE/110253‑110427

gAAAAAACATTGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTa                                                                                      >  2:3590312/1‑91 (MQ=255)

     aaCATTGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGt                                                                                 <  2:2462535/91‑1 (MQ=255)

       cATTGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTAc                                                                               >  1:5364478/1‑91 (MQ=255)

       cATTGAGTATAAAAACAATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTAc                                                                               <  2:5167226/91‑1 (MQ=17)

        aTTGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACa                                                                              >  2:2823480/1‑91 (MQ=255)

         ttGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACAt                                                                             <  1:435130/91‑1 (MQ=255)

          tGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATg                                                                            >  1:5979513/1‑91 (MQ=255)

          tGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATg                                                                            >  1:3235998/1‑91 (MQ=255)

           gAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATgg                                                                           <  1:5939103/91‑1 (MQ=255)

           gAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATgg                                                                           >  2:4327852/1‑91 (MQ=255)

            aGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGa                                                                          >  1:2243594/1‑91 (MQ=255)

            aGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGa                                                                          <  2:4256274/91‑1 (MQ=255)

             gTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGAt                                                                         >  1:505779/1‑91 (MQ=255)

             gTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGAt                                                                         <  1:6388400/91‑1 (MQ=255)

             gTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGAt                                                                         <  1:2472544/91‑1 (MQ=255)

              tATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGAtt                                                                        <  1:6957522/91‑1 (MQ=255)

               aTGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATta                                                                       <  1:6476384/91‑1 (MQ=255)

                tGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATtat                                                                      <  1:6365772/91‑1 (MQ=255)

                  tACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATaa                                                                    <  2:6660199/91‑1 (MQ=255)

                    cccTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAAt                                                                  <  1:3705645/91‑1 (MQ=255)

                     ccTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAAtt                                                                 <  1:3468209/91‑1 (MQ=255)

                     ccTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAAtt                                                                 >  1:6893419/1‑91 (MQ=255)

                     ccTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAAtt                                                                 <  1:1917202/91‑1 (MQ=255)

                      cTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTc                                                                <  1:662073/91‑1 (MQ=255)

                       tATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCt                                                               <  2:6449864/91‑1 (MQ=255)

                         ttAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGa                                                             >  2:1938924/1‑91 (MQ=255)

                         ttAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGa                                                             <  1:4812405/91‑1 (MQ=255)

                         ttAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGa                                                             <  2:6089455/91‑1 (MQ=255)

                           aaaaaTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGaaa                                                           >  2:2888850/1‑91 (MQ=255)

                            aaaaTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAAt                                                          >  1:2079009/1‑91 (MQ=255)

                             aaaTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATg                                                         >  2:2333270/1‑91 (MQ=255)

                             aaaTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATg                                                         >  2:4436729/1‑91 (MQ=255)

                             aaaTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATg                                                         <  2:4985548/91‑1 (MQ=255)

                              aaTGTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATgg                                                        <  1:1359753/91‑1 (MQ=255)

                              aaTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATgg                                                        <  2:5283328/91‑1 (MQ=255)

                              aaTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATgg                                                        <  2:3169629/91‑1 (MQ=255)

                               aTCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGt                                                       <  1:2365165/91‑1 (MQ=255)

                                tCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGtt                                                      <  2:4952301/91‑1 (MQ=255)

                                  tAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTAt                                                    >  2:4784505/1‑91 (MQ=255)

                                   agagGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTAtt                                                   <  1:1055433/91‑1 (MQ=255)

                                    gagGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTc                                                  >  1:3607497/1‑91 (MQ=255)

                                    gagGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTc                                                  <  1:682564/91‑1 (MQ=255)

                                     agGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCa                                                 >  1:618915/1‑91 (MQ=255)

                                         tttATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCt                                             <  1:4157251/91‑1 (MQ=255)

                                         tttATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCt                                             <  1:4662881/91‑1 (MQ=255)

                                            ataCGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCtaat                                          >  1:919740/1‑91 (MQ=255)

                                             taCGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCtaata                                         >  2:717699/1‑91 (MQ=255)

                                              aCGATAAGGAAGGCTATCGTGATGATCAGTTTGATCAAAATGTTACAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCtaataa                                        >  2:6779102/1‑91 (MQ=255)

                                              aCGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCtaataa                                        <  1:305686/91‑1 (MQ=255)

                                              aCGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCtaataa                                        >  2:451189/1‑91 (MQ=255)

                                               cGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTaataat                                       <  1:388082/91‑1 (MQ=255)

                                                gATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATg                                      <  1:5533742/91‑1 (MQ=255)

                                                 aTAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGa                                     <  2:1326791/91‑1 (MQ=255)

                                                  tAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGaa                                    >  2:802092/1‑91 (MQ=255)

                                                   aaggaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAg                                   >  2:30158/1‑91 (MQ=255)

                                                    aggaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGa                                  <  2:4020398/91‑1 (MQ=255)

                                                     ggaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGat                                 <  2:3785229/91‑1 (MQ=255)

                                                     ggaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGat                                 >  1:6690461/1‑91 (MQ=255)

                                                     ggaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGat                                 >  2:3088054/1‑91 (MQ=255)

                                                     ggaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGat                                 >  1:713367/1‑91 (MQ=255)

                                                      gaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGata                                >  1:4266851/1‑91 (MQ=255)

                                                      gaaggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGata                                >  2:5965036/1‑91 (MQ=255)

                                                        aggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATg                              >  2:5396277/1‑91 (MQ=255)

                                                         ggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATgg                             <  1:82187/91‑1 (MQ=255)

                                                         ggCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATgg                             <  1:330968/91‑1 (MQ=255)

                                                         ggCTAACGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATgg                             >  2:5468746/1‑91 (MQ=255)

                                                          gCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGt                            >  2:3272705/1‑91 (MQ=255)

                                                          gCTATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGt                            >  2:1689583/1‑91 (MQ=255)

                                                           cgaTCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTa                           <  2:2779267/89‑1 (MQ=255)

                                                            tATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAg                          <  1:4003156/91‑1 (MQ=255)

                                                            tATCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAg                          <  1:5958007/91‑1 (MQ=255)

                                                             aTCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGc                         <  2:1418133/91‑1 (MQ=255)

                                                             aTCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGc                         <  1:5642048/91‑1 (MQ=255)

                                                              tCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGCTATGGTAGCt                        >  2:5567482/1‑91 (MQ=255)

                                                              tCGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTCAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCt                        >  1:5805119/1‑91 (MQ=255)

                                                               cGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTa                       <  2:1062560/91‑1 (MQ=255)

                                                               cGTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTa                       <  2:4173010/91‑1 (MQ=255)

                                                                gTGATGATCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTaa                      >  2:3453296/1‑91 (MQ=255)

                                                                 tgatgaTCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTaaa                     <  2:1831429/91‑1 (MQ=255)

                                                                  gatgatCAGTTTGATAAAAATGTTAAAGGTACATGGATTAT‑AATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAgg                   >  2:2285490/1‑91 (MQ=255)

                                                                   atgatCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAgg                   <  2:3721579/91‑1 (MQ=255)

                                                                    tgatCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGc                  <  1:5261593/91‑1 (MQ=255)

                                                                    tgatCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGc                  <  1:3270852/91‑1 (MQ=255)

                                                                     gatCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCa                 <  2:1554667/91‑1 (MQ=255)

                                                                      atCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCAt                <  1:2657125/91‑1 (MQ=255)

                                                                      atCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCAt                <  2:5822196/91‑1 (MQ=255)

                                                                       tCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATg               >  2:4621193/1‑91 (MQ=255)

                                                                       tCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATg               <  1:1764842/91‑1 (MQ=255)

                                                                       tCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATg               >  1:6064363/1‑91 (MQ=255)

                                                                       tCAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATg               >  2:5780276/1‑91 (MQ=255)

                                                                        cAGTTTGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATgg              <  2:6320739/91‑1 (MQ=255)

                                                                           tttGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTc           <  1:4912345/91‑1 (MQ=255)

                                                                           tttGATAAAAATGATAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTc           <  2:5343326/91‑1 (MQ=255)

                                                                            ttGATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCt          <  1:6271476/91‑1 (MQ=255)

                                                                              gATAAAAATGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCtat        >  1:286648/1‑91 (MQ=255)

                                                                                tAAAAATGTTAAAGGTACATGGATATTAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGTATGGTTCtatat      <  2:6522886/91‑1 (MQ=37)

                                                                                    aaTGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCtat        >  1:4075463/1‑85 (MQ=255)

                                                                                    aaTGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCTATATATGa  >  2:2048372/1‑91 (MQ=255)

                                                                                    aaTGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCTATATATGa  >  2:694769/1‑91 (MQ=255)

                                                                                    aaTGTTAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCTATATATGa  <  1:1527704/91‑1 (MQ=255)

GAAAAAACATTGAGTATGTACCCTATTAAAAATCTAGAGGATTTATACGATAAGGAAGGCTATCGTGATGATCAGTTTGATAAAAATGATAAAGGTACATGGATTATAAATTCTGAAATGGTTATTCAACCTAATAATGAAGATATGGTAGCTAAAGGCATGGTTCTATATATGA  >  USA300TCH1516_ALE/110253‑110427

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