BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	341,450	G→C	14.2%	D134H (GAT→CAT)	USA300HOU_RS01525 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	341,450	0	G	C	14.2%	630.0 / 83.1	219	D134H (GAT→CAT)	USA300HOU_RS01525	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (105/82); new base C (16/15); total (121/98)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 231 total aligned reads displayed.

GGGTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGAAA  >  USA300TCH1516_ALE/341360‑341540

gggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATg                                                                                            >  1:1592677/1‑91 (MQ=255)

 ggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGa                                                                                           >  1:992270/1‑91 (MQ=255)

   tttCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATCAtt                                                                                         <  2:1926578/91‑1 (MQ=255)

     tCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTgg                                                                                       <  2:2161567/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                       <  1:3760113/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                       <  1:899930/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTggg                                                                                      >  2:1134047/1‑91 (MQ=255)

        aaTACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTggac                                                                                     <  1:1896339/91‑3 (MQ=255)

         aTACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTggact                                                                                    <  1:829789/91‑4 (MQ=255)

           aCGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACttt                                                                                 <  2:745481/91‑1 (MQ=255)

            cGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTa                                                                                <  1:723052/91‑1 (MQ=255)

             gAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATT‑GGACttttc                                                                               <  1:1211292/91‑3 (MQ=38)

             gAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaa                                                                               <  2:608307/91‑1 (MQ=255)

             gAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaa                                                                               <  1:1883514/91‑1 (MQ=255)

             gAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaa                                                                               >  1:2709836/1‑91 (MQ=255)

             gAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaa                                                                               >  2:1268284/1‑91 (MQ=255)

              aaTTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACttttct                                                                              <  1:3455834/91‑4 (MQ=255)

               aTTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAAc                                                                             >  2:2773457/1‑91 (MQ=255)

               aTTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAAAGATTGGGACTTTAAAc                                                                             >  2:3717173/1‑91 (MQ=255)

                 tATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACaa                                                                           <  2:3011079/91‑1 (MQ=255)

                      aTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGa                                                                      >  2:1799118/1‑91 (MQ=255)

                      aTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGa                                                                      >  1:3381958/1‑91 (MQ=255)

                       tAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGAc                                                                     >  2:1387834/1‑91 (MQ=255)

                         aaGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTAGTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGa                                                                   <  2:3421447/91‑1 (MQ=255)

                           gACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATg                                                                 >  1:2840352/1‑91 (MQ=255)

                            aCATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGa                                                                >  1:886068/1‑91 (MQ=255)

                            aCATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGa                                                                >  1:1834412/1‑91 (MQ=255)

                            aCATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGa                                                                <  2:3184843/91‑1 (MQ=255)

                              aTTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGaaa                                                              >  1:2082270/1‑91 (MQ=255)

                              aTTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGaaa                                                              >  2:3391218/1‑91 (MQ=255)

                               ttGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAAc                                                             >  2:2313715/1‑91 (MQ=255)

                               ttGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAAc                                                             >  1:3828948/1‑91 (MQ=255)

                                tGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACa                                                            >  2:2558416/1‑91 (MQ=255)

                                 gCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAt                                                           <  1:1032071/91‑1 (MQ=255)

                                 gCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAt                                                           >  1:1435651/1‑91 (MQ=255)

                                  cATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAtt                                                          >  2:1748113/1‑91 (MQ=255)

                                  cATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAtt                                                          <  1:678270/91‑1 (MQ=255)

                                      gtggtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGgag                                                      <  1:1316692/91‑1 (MQ=255)

                                           cGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACtt                                                 >  1:2471877/1‑91 (MQ=255)

                                           cGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACtt                                                 >  1:55094/1‑91 (MQ=255)

                                            gATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTa                                                >  1:1595213/1‑91 (MQ=255)

                                              ttGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaaa                                              >  2:280244/1‑91 (MQ=255)

                                              ttGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaaa                                              >  2:2080033/1‑91 (MQ=255)

                                                gTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaata                                            <  1:3055027/91‑1 (MQ=255)

                                                 tACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaataa                                           >  2:1693292/1‑91 (MQ=255)

                                                  aCAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaataat                                          >  2:2730584/1‑91 (MQ=255)

                                                    aaTCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGa                                        >  2:2334314/1‑91 (MQ=255)

                                                     aTCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGaa                                       <  2:849530/91‑1 (MQ=255)

                                                       cATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAg                                     <  2:3677436/91‑1 (MQ=255)

                                                          atCCTTTATTGCTGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGCAAGTag                                  >  2:3861267/1‑90 (MQ=255)

                                                          atCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTAt                                  >  2:2196219/1‑91 (MQ=255)

                                                           tCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATc                                 >  2:2834913/1‑91 (MQ=255)

                                                            ccTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTa                                                >  1:2352567/1‑75 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               <  2:903519/91‑1 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               <  2:505967/91‑1 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               <  2:1477810/91‑1 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               <  1:3288447/91‑1 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               <  2:932659/91‑1 (MQ=255)

                                                                tATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATAc                            <  1:505463/91‑1 (MQ=255)

                                                                 aTTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACg                           <  1:2759944/91‑1 (MQ=255)

                                                                  ttGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGa                          >  1:3723895/1‑91 (MQ=255)

                                                                    gatgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGAtt                        >  1:1381932/1‑91 (MQ=255)

                                                                      tgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTgg                      >  1:3192541/1‑91 (MQ=255)

                                                                      tgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTgg                      >  2:3039417/1‑91 (MQ=255)

                                                                      tgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTgg                      >  2:2467213/1‑91 (MQ=255)

                                                                       gagaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTggg                     <  1:1420585/91‑1 (MQ=255)

                                                                       gagaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTggg                     >  1:3857930/1‑91 (MQ=255)

                                                                        agaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGc                    >  1:642549/1‑91 (MQ=255)

                                                                        agaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGc                    <  2:1558025/91‑1 (MQ=255)

                                                                         gaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCg                   >  1:2673040/1‑91 (MQ=255)

                                                                         gaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCg                   >  2:3374542/1‑91 (MQ=255)

                                                                           cAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCgtg                 <  1:3031284/91‑1 (MQ=255)

                                                                            aaGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACAATTGGGCgtgt                <  2:1024632/91‑1 (MQ=255)

                                                                             aGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGtt               <  1:2207886/91‑1 (MQ=255)

                                                                              gAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGttt              >  1:496134/1‑91 (MQ=255)

                                                                                aaaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTa            >  1:2234385/1‑91 (MQ=255)

                                                                                 aaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAg           <  1:418307/91‑1 (MQ=255)

                                                                                 aaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAg           <  2:2182701/91‑1 (MQ=255)

                                                                                 aaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAg           >  1:2106421/1‑91 (MQ=255)

                                                                                  agtcccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGa          >  1:3504530/5‑91 (MQ=255)

                                                                                  aGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGa          >  1:3689035/1‑91 (MQ=255)

                                                                                   gT‑CCAATGATTGGGACTTTAAACAAATAGACGAAGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGaa         >  1:1708660/1‑91 (MQ=255)

                                                                                      ccAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAt       <  2:2891659/91‑1 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  2:726908/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  2:522098/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  2:959450/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  2:244898/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  1:561593/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  1:3333052/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  1:1500689/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  1:1358195/1‑91 (MQ=255)

                                                                                       cAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      <  1:2890372/91‑1 (MQ=255)

                                                                                       cAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      <  1:2593919/91‑1 (MQ=255)

                                                                                        aaTGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTg     >  2:1942815/1‑91 (MQ=255)

                                                                                         aTGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGa    >  2:2911964/1‑91 (MQ=255)

                                                                                         aTGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGa    <  1:1133047/91‑1 (MQ=255)

                                                                                          tGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaa   <  2:2399206/91‑1 (MQ=255)

                                                                                           gATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaaa  >  2:3337626/1‑91 (MQ=255)

                                                                                           gATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaaa  >  1:1529611/1‑91 (MQ=255)

                                                                                           gATTGGGACTTAAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaaa  <  2:3499400/91‑1 (MQ=255)

GGGTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGAAA  >  USA300TCH1516_ALE/341360‑341540

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