BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	479,553	G→A	A194T (GCA→ACA)	USA300HOU_RS02280 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	479,553	0	G	A	94.0%	528.7 / 13.6	167	A194T (GCA→ACA)	USA300HOU_RS02280	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (7/3); new base A (83/74); total (90/77)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

ATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATGCATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGA  >  USA300TCH1516_ALE/479464‑479640

aTCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAG‑TTTTACcataca                                                                                       <  2:1137624/91‑1 (MQ=255)

  cGTCAAATGATTAAATGGACGGAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACcatacat                                                                                      <  2:1244474/91‑1 (MQ=255)

  cGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACcatacat                                                                                      <  2:279165/91‑1 (MQ=255)

  cGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACcatacat                                                                                      <  2:449827/91‑1 (MQ=255)

   gTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACcatacatg                                                                                     <  2:2344092/91‑1 (MQ=255)

     cAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACcatacatgta                                                                                   <  1:612427/91‑3 (MQ=255)

     cAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGt                                                                                   <  1:1624101/91‑1 (MQ=255)

      aaaTGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGtt                                                                                  <  1:3436369/91‑1 (MQ=255)

      aaaTGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGtt                                                                                  <  1:3438705/91‑1 (MQ=255)

        aTGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGc                                                                                <  1:880183/91‑1 (MQ=255)

        aTGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGc                                                                                >  2:2065747/1‑91 (MQ=255)

          gATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATGCATGGTTGCAt                                                                              <  1:1637244/91‑1 (MQ=255)

          gATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCAt                                                                              <  1:906719/91‑1 (MQ=255)

           aTTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCAtt                                                                             <  2:1008418/91‑1 (MQ=255)

            ttAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCAttt                                                                            >  2:2485441/1‑91 (MQ=255)

              aaaTGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAg                                                                          <  2:2919942/91‑1 (MQ=255)

               aaTGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAgc                                                                         >  1:2135403/1‑91 (MQ=255)

               aaTGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAgc                                                                         >  2:1566474/1‑91 (MQ=255)

               aaTGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAgc                                                                         <  1:2403517/91‑1 (MQ=255)

                aTGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAgcg                                                                        <  2:577050/91‑1 (MQ=255)

                 tGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAgcgc                                                                       <  1:864289/91‑1 (MQ=255)

                    aCGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAAc                                                                    <  2:2249239/91‑1 (MQ=255)

                        aaCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGa                                                                >  2:1932804/1‑91 (MQ=255)

                        aaCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGa                                                                <  2:748296/91‑1 (MQ=255)

                        aaCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGa                                                                <  2:748490/91‑1 (MQ=255)

                         aCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGAc                                                               >  1:1455415/1‑91 (MQ=255)

                         aCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGAc                                                               >  2:294217/1‑91 (MQ=255)

                          cTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGAcc                                                              <  2:2722184/91‑1 (MQ=255)

                          cTTTATATCGATCAAGTTCTATCGATCTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGAcc                                                              <  2:3147892/91‑1 (MQ=255)

                           tttATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCa                                                             <  2:3123982/91‑1 (MQ=255)

                           ttgatatCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCa                                                             <  2:3523079/88‑1 (MQ=255)

                            ttATATCGATCAATTTCTATCGAGTTAGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCAt                                                            >  2:487567/1‑91 (MQ=255)

                            ttAGATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCAt                                                            <  2:1404479/91‑1 (MQ=255)

                               tatCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGt                                                         <  2:1121441/91‑1 (MQ=255)

                               tatCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGt                                                         >  1:1456757/1‑91 (MQ=255)

                                 tCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGttt                                                       >  2:2123832/1‑91 (MQ=255)

                                    aTCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTAc                                                    >  1:139532/1‑91 (MQ=255)

                                     tCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATGCATGGTTGCATTTAGCGCAACATGACCATAGTTTTACt                                                   >  1:576486/1‑91 (MQ=255)

                                      cAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTa                                                  <  2:295900/91‑1 (MQ=255)

                                       aaTTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTaa                                                 >  1:2452001/1‑91 (MQ=255)

                                        atatCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTaaa                                                <  1:966970/88‑1 (MQ=255)

                                        aTTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTaaa                                                >  1:2612417/1‑91 (MQ=255)

                                         tttCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAg                                               <  1:2425626/91‑1 (MQ=255)

                                           tCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGca                                             <  2:1763245/91‑1 (MQ=255)

                                           tCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGca                                             >  1:2815439/1‑91 (MQ=255)

                                           tCTATCGAGTTGGACAATGCCGAAGCCAGAGCAAAGTTTTTACCATGCATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGca                                             >  2:2064916/1‑91 (MQ=255)

                                            ctaccGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGcac                                            <  1:1931958/87‑1 (MQ=255)

                                            cTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGcac                                            >  2:2384581/1‑91 (MQ=255)

                                            cTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGcac                                            <  2:379794/91‑1 (MQ=255)

                                             tATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGcaca                                           >  2:578113/1‑91 (MQ=255)

                                             tATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACCAAGAcaca                                           >  2:3336988/1‑91 (MQ=255)

                                              aTCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAg                                          <  2:2852121/91‑1 (MQ=255)

                                              aTCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAg                                          >  1:1542344/1‑91 (MQ=255)

                                                cGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAgcg                                        <  2:2676940/91‑1 (MQ=255)

                                                 gAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAgcgc                                       >  2:1588814/1‑91 (MQ=255)

                                                 gAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAgcgc                                       >  1:3073603/1‑91 (MQ=255)

                                                  aGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGcc                                      >  1:3028340/1‑91 (MQ=255)

                                                    ttGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCaa                                    <  2:1186506/91‑1 (MQ=255)

                                                     tGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAg                                   <  1:477276/91‑1 (MQ=255)

                                                      ggACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGt                                  <  2:3442004/91‑1 (MQ=255)

                                                       gACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTg                                 <  2:1485933/91‑1 (MQ=255)

                                                        aCAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGAGCCAAGTGa                                >  1:3468894/1‑91 (MQ=255)

                                                         cAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGAt                               <  1:2335904/91‑1 (MQ=255)

                                                         cAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGAt                               <  2:2194242/91‑1 (MQ=255)

                                                          aaTGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGAtt                              >  1:666437/1‑91 (MQ=255)

                                                          aaTGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGAtt                              <  1:3542785/91‑1 (MQ=255)

                                                           aTGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTa                             >  1:2733894/1‑91 (MQ=255)

                                                           aTGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTa                             >  2:2010164/1‑91 (MQ=255)

                                                           aTGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTa                             <  1:3063648/91‑1 (MQ=255)

                                                           aTGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTa                             >  2:702622/1‑91 (MQ=255)

                                                            tGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTaa                            <  2:3252719/91‑1 (MQ=255)

                                                            tGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTaa                            >  2:1669066/1‑91 (MQ=255)

                                                            tGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTaa                            <  2:1557210/91‑1 (MQ=255)

                                                            tGCCGAAGCGTGAGCAAAGTTTTGCCCATACATGGTTGCATTTGGCGCAATATGGCCATAGTTTTACTAAAGCACAGCGCCAACTGATTaa                            <  2:3376948/91‑1 (MQ=255)

                                                            agCCGAAGCGTGAGCAAAGTTTTTACCATGCATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTaa                            <  2:3072277/90‑1 (MQ=255)

                                                             gCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTaaa                           >  1:1684743/1‑91 (MQ=255)

                                                             gCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTaaa                           <  1:3249026/91‑1 (MQ=255)

                                                               cGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAgg                         <  1:698199/91‑1 (MQ=255)

                                                                gAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGc                        >  1:3324317/1‑91 (MQ=255)

                                                                gAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGc                        >  1:1856762/1‑91 (MQ=255)

                                                                 aaGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCt                       >  2:2271632/1‑91 (MQ=255)

                                                                  aGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTATAGCACAGCGCCAAGTGATTAAAGGCtt                      >  2:944085/1‑91 (MQ=255)

                                                                    cGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTAc                    >  2:630936/1‑91 (MQ=255)

                                                                      tGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTAccc                  >  1:1610526/1‑91 (MQ=255)

                                                                      tGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTAccc                  >  1:1611465/1‑91 (MQ=255)

                                                                         gCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAAt               <  2:886233/91‑1 (MQ=255)

                                                                          cAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATg              <  2:2045618/91‑1 (MQ=255)

                                                                          cAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATg              <  2:799031/91‑1 (MQ=255)

                                                                           aaaGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGa             <  1:1001077/91‑1 (MQ=255)

                                                                            aaGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGCTCAAAGGCTTACCCAATGac            >  1:1009326/1‑90 (MQ=255)

                                                                             aGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCTAGTGATTAAAGGCTTACCCAATGATc           <  2:2064622/91‑1 (MQ=255)

                                                                             aGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATc           <  1:819825/91‑1 (MQ=255)

                                                                             aGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATc           <  1:1588814/91‑1 (MQ=255)

                                                                              gTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATcc          >  1:3417989/1‑91 (MQ=255)

                                                                              gTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATcc          >  2:242107/1‑91 (MQ=255)

                                                                                  ttACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGaa      >  1:1217097/1‑91 (MQ=255)

                                                                                   tACCATACATGGCTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGaaa     >  1:2788210/1‑91 (MQ=255)

                                                                                     ccATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATg   >  1:3000786/1‑91 (MQ=255)

                                                                                     ccATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATg   <  1:714893/91‑1 (MQ=255)

                                                                                      catacatgGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGa  <  2:1082216/91‑1 (MQ=255)

ATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATGCATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGA  >  USA300TCH1516_ALE/479464‑479640

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