BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,192,759:1	+G	5.1%	coding (378/480 nt)	USA300HOU_RS11200 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,192,759	1	.	G	5.1%	924.6 / 44.7	216	coding (378/480 nt)	USA300HOU_RS11200	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (118/87); new base G (8/3); total (126/90)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 232 total aligned reads displayed.

TTCAACCCCTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGACTCGTT  >  USA300TCH1516_ALE/2192669‑2192849

 tCAACCCCTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCGTGAGACATCATCGGTAAATAAATAGTATGA‑CTGCA‑gg                                                                                          <  2:652310/91‑2 (MQ=255)

 tCAACCCCTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑g                                                                                           >  2:944500/1‑91 (MQ=255)

 tCAACCCCTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑g                                                                                           >  1:31550/1‑91 (MQ=255)

 tCAACCCCTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑g                                                                                           <  2:496432/91‑1 (MQ=255)

  cAACCCCTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGA‑CTGCA‑ggt                                                                                         <  1:2541343/91‑3 (MQ=255)

   aaCCCCTCTCAAATTCGACATATTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTc                                                                                         >  1:149637/1‑91 (MQ=255)

     ccccTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCAt                                                                                       >  1:2951379/1‑91 (MQ=255)

      cccTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATa                                                                                      <  1:2215135/91‑1 (MQ=255)

        ctctCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATaaa                                                                                    >  2:2823162/1‑91 (MQ=255)

        ctctCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATaaa                                                                                    >  1:1328002/1‑91 (MQ=255)

         tctcAAATTCGGCATAGTGCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGA‑CTGCAGGTCATacta                                                                                   <  2:1739817/91‑4 (MQ=255)

         tctcAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAAt                                                                                   <  2:1039901/91‑1 (MQ=255)

          ctcAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATc                                                                                  >  2:3210479/1‑91 (MQ=255)

            cAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCa                                                                                <  1:2726572/91‑1 (MQ=255)

                 tCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACttt                                                                           >  1:2534317/1‑91 (MQ=255)

                  cGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTta                                                                          <  2:3155302/91‑1 (MQ=255)

                  cGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTta                                                                          >  1:1415819/1‑91 (MQ=255)

                    aCATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTtata                                                                        <  2:2075257/91‑1 (MQ=255)

                    aCATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTtata                                                                        <  2:3261260/91‑1 (MQ=255)

                    aCATAGTTCTCTCTTCAATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCATCTTtata                                                                        <  2:3079234/91‑1 (MQ=255)

                       tAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACaa                                                                     <  1:3093430/91‑1 (MQ=255)

                        aGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACaaa                                                                    >  1:295761/1‑91 (MQ=255)

                         gTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAAt                                                                   <  2:1953830/91‑1 (MQ=255)

                           tctctcTTCTATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTa                                                                 <  1:3798861/91‑1 (MQ=255)

                            ctctctTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTaa                                                                <  1:2973755/91‑1 (MQ=255)

                                ctTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCAc                                                            >  1:685401/1‑91 (MQ=255)

                                ctTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCAc                                                            >  2:1364023/1‑91 (MQ=255)

                                ctTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCAc                                                            >  1:61806/1‑91 (MQ=255)

                                  tCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACtt                                                          <  2:3508537/91‑1 (MQ=255)

                                  tCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACtt                                                          <  2:1578937/91‑1 (MQ=255)

                                      ttattTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACt                                                      <  1:2409751/91‑1 (MQ=255)

                                       tattTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTa                                                     >  1:2503783/1‑91 (MQ=255)

                                        attTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTaa                                                    >  1:2993905/1‑91 (MQ=255)

                                        attTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTaa                                                    >  1:1832006/1‑91 (MQ=255)

                                          tttAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATTCAACTTTATACAAATTAAGCACTTTACTAAtt                                                  >  1:3355312/1‑91 (MQ=255)

                                          tttAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAAtt                                                  >  2:1523881/1‑91 (MQ=255)

                                          tttAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAAtt                                                  >  2:1544986/1‑91 (MQ=255)

                                           ttAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTg                                                 <  1:3395279/91‑1 (MQ=255)

                                            tAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTgg                                                >  1:2866741/1‑91 (MQ=255)

                                             aaCATCGTCATGAGACATCATCGGTAAATAAATATTATGACCTGCA‑GTCATAACTCCAACTTTATACAAATTAAGCACTTTACTACTTGGa                                               >  2:1412876/1‑91 (MQ=255)

                                             aaCATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGa                                               >  1:1318394/1‑91 (MQ=255)

                                             aaCATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGa                                               >  1:450941/1‑91 (MQ=255)

                                             aaCATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGa                                               >  1:1669846/1‑91 (MQ=255)

                                               cATCGTCATGAGACATCATCGGTAAATTAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGAtt                                             >  2:723676/1‑91 (MQ=255)

                                                 tCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAg                                           >  1:2087542/1‑91 (MQ=255)

                                                     cATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGAttt                                       <  2:2301792/91‑1 (MQ=255)

                                                      aTGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTa                                      >  2:3634606/1‑91 (MQ=255)

                                                        gagaCATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAAt                                    <  1:521983/91‑1 (MQ=255)

                                                        gagaCATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAAt                                    <  2:116291/91‑1 (MQ=255)

                                                        gagaCATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAAt                                    <  1:3267994/91‑1 (MQ=255)

                                                        gagaCATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGAATTAAt                                    >  2:3645557/1‑91 (MQ=255)

                                                          gaCATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCg                                  >  1:2340034/1‑91 (MQ=255)

                                                               catcGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAg                             <  1:3772105/91‑1 (MQ=255)

                                                                 tcGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTa                           >  2:2518786/1‑91 (MQ=255)

                                                                 tcGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTa                           >  1:301762/1‑91 (MQ=255)

                                                                  cGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTAt                          >  2:1552229/1‑91 (MQ=255)

                                                                   ggTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTAtt                         >  1:1142370/1‑91 (MQ=255)

                                                                    gTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTg                        >  2:2175288/1‑91 (MQ=255)

                                                                     taaataaaTAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGt                       >  1:1131008/1‑91 (MQ=255)

                                                                      aaataaatAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTa                      >  2:239510/1‑91 (MQ=255)

                                                                       aataaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTaa                     >  2:1199660/1‑91 (MQ=255)

                                                                       aataaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTaa                     <  2:2977103/91‑1 (MQ=255)

                                                                        ataaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCGCTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTaaa                    <  2:2654192/91‑1 (MQ=255)

                                                                        ataaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTaaa                    >  2:1706681/1‑91 (MQ=255)

                                                                        ataaataGTATGACCTGCA‑GTAATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTaaa                    >  1:1173480/1‑91 (MQ=255)

                                                                         tacataGTATGACCTGCC‑GTCATAAATCCAACTCTATACCAATTCAGCACTCTACTAATTGGATCCGATTTAATCGACACGTATTGTAAAc                   >  2:1053679/4‑91 (MQ=255)

                                                                         taaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAAc                   >  2:967878/1‑91 (MQ=255)

                                                                         taaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAAc                   >  1:1600630/1‑91 (MQ=255)

                                                                         taaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAAc                   >  1:1842715/1‑91 (MQ=255)

                                                                          aaataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACg                  >  2:3783894/1‑91 (MQ=255)

                                                                           aataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGt                 <  2:3740081/91‑1 (MQ=255)

                                                                           aataGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGt                 <  2:332891/91‑1 (MQ=255)

                                                                             taGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTc               >  1:1006247/1‑91 (MQ=255)

                                                                             taGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTc               >  1:939195/1‑91 (MQ=255)

                                                                             taGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACCTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTc               >  1:3133978/1‑91 (MQ=255)

                                                                              aGTATGACCTGCA‑GTCATAAATCCAGCTTTATACAGATTAAGCACTTTACTAAGTGGATTAGATTTAATCGACAAGTATGGTAAACGTTCa              <  2:2690249/91‑1 (MQ=255)

                                                                              aGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCa              <  1:949315/91‑1 (MQ=255)

                                                                              aGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCa              <  2:1898942/91‑1 (MQ=255)

                                                                               gTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCaa             >  1:2332352/1‑91 (MQ=255)

                                                                               gTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCaa             >  2:3911043/1‑91 (MQ=255)

                                                                                tATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAAt            <  2:601225/91‑1 (MQ=255)

                                                                                tATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAAt            <  1:3204885/91‑1 (MQ=255)

                                                                                tATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAAt            <  1:1859573/91‑1 (MQ=255)

                                                                                 aTGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAAtt           >  1:796262/1‑91 (MQ=255)

                                                                                 aTGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAAtt           >  1:1954442/1‑91 (MQ=255)

                                                                                   tgacTGCAGGTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTc          >  1:3495529/4‑91 (MQ=255)

                                                                                   gACCTGCA‑GTCATAAATCCATCTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCg         <  1:703233/91‑1 (MQ=255)

                                                                                   gACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCg         <  1:2119005/91‑1 (MQ=255)

                                                                                   gACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCg         >  2:454367/1‑91 (MQ=255)

                                                                                     acTGCAGGTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGa        >  2:3369490/2‑91 (MQ=255)

                                                                                     acTGCAGGTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGa        >  1:2454813/2‑91 (MQ=255)

                                                                                      cggCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGACt      <  2:783084/89‑1 (MQ=255)

                                                                                      cTGCAGGTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGAc       >  1:924410/1‑91 (MQ=255)

                                                                                      cTGCAGGTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGAc       >  2:2832895/1‑91 (MQ=255)

                                                                                      cTGCAGGTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGAc       >  2:2229508/1‑91 (MQ=255)

                                                                                      cTGCAGGTCATAAATCCAACTTTATAAAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATACGAc       >  2:2237183/1‑91 (MQ=255)

                                                                                      cTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGACt      <  2:2016765/91‑1 (MQ=255)

                                                                                       tGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGACTc     <  2:3607664/91‑1 (MQ=255)

                                                                                       tGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGACTc     <  2:272253/91‑1 (MQ=255)

                                                                                         cA‑GTCATAAATCCAACTTTATACAAATGAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGACTCGt   <  2:936945/91‑1 (MQ=255)

TTCAACCCCTCTCAAATTCGACATAGTTCTCTCTTCGATTATTTTAACATCGTCATGAGACATCATCGGTAAATAAATAGTATGACCTGCA‑GTCATAAATCCAACTTTATACAAATTAAGCACTTTACTAATTGGATTAGATTTAATCGACAAGTATTGTAAACGTTCAATTCGACTCGTT  >  USA300TCH1516_ALE/2192669‑2192849

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