BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,038,119	A→C	6.0%	T236P (ACG→CCG)	USA300TCH1516_00964 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,038,119	0	A	C	6.0%	464.6 / 17.2	132	T236P (ACG→CCG)	USA300TCH1516_00964	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (69/55); new base C (3/5); total (72/60)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

GAAAATTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCT  >  USA300TCH1516_ALE/1038029‑1038209

gAAAATTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTaca                                                                                            >  1:3785697/1‑91 (MQ=255)

 aaaaTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACATGTGTACTTTGGATATTAAATTacac                                                                                           >  1:3849427/1‑91 (MQ=255)

 aaaaTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTacac                                                                                           <  1:1632554/91‑1 (MQ=255)

 aaaaTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTacac                                                                                           <  1:1968161/91‑1 (MQ=255)

   aaTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACgg                                                                                         >  1:2380627/1‑91 (MQ=255)

    aTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACggg                                                                                        <  1:3430467/91‑1 (MQ=255)

    aTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACggg                                                                                        <  1:1843942/91‑1 (MQ=255)

     tttttCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGt                                                                                       >  1:979416/1‑91 (MQ=255)

     tttttCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGt                                                                                       >  1:977995/1‑91 (MQ=255)

       tttCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACCCGTGTaa                                                                                     <  1:2753409/91‑1 (MQ=255)

       tttCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTaa                                                                                     >  1:2866426/1‑91 (MQ=255)

        ttCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAAt                                                                                    >  1:381155/1‑91 (MQ=255)

        ttCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAAt                                                                                    >  1:237708/1‑91 (MQ=255)

         tCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACCCGTGTAAtt                                                                                   <  1:775168/91‑1 (MQ=255)

         tCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACCCGTGTAAtt                                                                                   <  1:450909/91‑1 (MQ=255)

          cAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACCCGTGTAAttt                                                                                  <  1:1150606/91‑1 (MQ=255)

          cAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAAttt                                                                                  >  1:1400163/1‑91 (MQ=255)

           aaaGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTa                                                                                 <  1:3095745/91‑1 (MQ=255)

            aaGAAATGCACGCATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAg                                                                                >  1:3200198/1‑91 (MQ=255)

              gAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGAt                                                                              <  1:785046/91‑1 (MQ=255)

                  tGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTAt                                                                          >  1:201598/1‑91 (MQ=255)

                  tGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTAt                                                                          <  1:275105/91‑1 (MQ=255)

                    cACGTATTGGGTCACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTc                                                                        >  1:3171979/1‑91 (MQ=255)

                    cACGTATTGGGTCACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTc                                                                        >  1:3171782/1‑91 (MQ=255)

                    cACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTc                                                                        <  1:1095049/91‑1 (MQ=255)

                    cACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTc                                                                        >  1:449213/1‑91 (MQ=255)

                     aCGTATTGGGTAACGTTAATCTTAACTCACATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCa                                                                       <  1:3596212/91‑1 (MQ=255)

                     aCGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCa                                                                       <  1:3391368/91‑1 (MQ=255)

                     aCGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCa                                                                       >  1:2669580/1‑91 (MQ=255)

                      cGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCat                                                                      <  1:3724225/91‑1 (MQ=255)

                         aTTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATAtt                                                                   <  1:2841936/91‑1 (MQ=255)

                          ttGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTa                                                                  <  1:2427390/91‑1 (MQ=255)

                            gggTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAg                                                                <  1:877037/91‑1 (MQ=255)

                                 aCGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAg                                                           >  1:1647649/1‑91 (MQ=255)

                                  cGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGt                                                          >  1:3631388/1‑91 (MQ=255)

                                   gTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGGTATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGtt                                                         <  1:4138168/91‑1 (MQ=255)

                                    ttAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTc                                                        <  1:2148992/91‑1 (MQ=255)

                                      aaTCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTg                                                      <  1:698961/91‑1 (MQ=255)

                                      aaTCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTg                                                      >  1:882045/1‑91 (MQ=255)

                                      aaTCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTg                                                      >  1:3121800/1‑91 (MQ=255)

                                      aaTCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTg                                                      >  1:1987644/1‑91 (MQ=255)

                                      aaTCATAACTCAAATTCTTTTAACTACAAGTATACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTg                                                      >  1:3682474/1‑91 (MQ=255)

                                           tAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTc                                                 >  1:503373/1‑91 (MQ=255)

                                           tAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTc                                                 <  1:1378395/91‑1 (MQ=255)

                                             aCTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCAt                                               <  1:621929/91‑1 (MQ=255)

                                             aCTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCAt                                               <  1:3214404/91‑1 (MQ=255)

                                               tCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGt                                             >  1:2940318/1‑91 (MQ=255)

                                                   aTTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCttt                                         >  1:2956640/1‑91 (MQ=255)

                                                   aTTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGCCttt                                         <  1:869650/91‑1 (MQ=255)

                                                    ttCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCtttt                                        <  1:3584233/91‑1 (MQ=255)

                                                    ttCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCtttt                                        <  1:233929/91‑1 (MQ=255)

                                                      cTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCa                                      >  1:389066/1‑91 (MQ=255)

                                                      cTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCa                                      >  1:1908129/1‑91 (MQ=255)

                                                        tttAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCAtt                                    >  1:2838529/1‑91 (MQ=255)

                                                          tAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCAtttt                                  <  1:2979038/91‑1 (MQ=255)

                                                          tAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCAtttt                                  <  1:2978818/91‑1 (MQ=255)

                                                          tAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCAtttt                                  >  1:2113236/1‑91 (MQ=255)

                                                           aaCTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTa                                 <  1:1376570/91‑1 (MQ=255)

                                                           aaCTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTa                                 <  1:957934/91‑1 (MQ=255)

                                                           aaCTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTa                                 <  1:3362276/91‑1 (MQ=255)

                                                            aCTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTaa                                >  1:2291203/1‑91 (MQ=255)

                                                            aCTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTaa                                >  1:1633083/1‑91 (MQ=255)

                                                               aCAAGTGTACTTTGGATATTAAATTACACGGATAATTTAGATTTATTCATATTAAGACTAGTTCTGCTTTCATGTCTTTTCATTTTAACGa                             >  1:3683461/1‑91 (MQ=255)

                                                                cAAGTGTACTTTGGATATTAAATTACACGGGTAATTTCGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGAt                            >  1:3507514/1‑91 (MQ=255)

                                                                cAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGAt                            <  1:1891115/91‑1 (MQ=255)

                                                                cAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGAt                            >  1:3507567/1‑91 (MQ=255)

                                                                 aaGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGAtt                           <  1:166319/91‑1 (MQ=255)

                                                                   gtgtACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAg                         >  1:2891941/1‑91 (MQ=255)

                                                                    tgtACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAgt                        >  1:3967444/1‑91 (MQ=255)

                                                                     gtACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAgtg                       >  1:74351/1‑91 (MQ=255)

                                                                      tACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAgtgt                      <  1:2237893/91‑1 (MQ=255)

                                                                      tACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAgtgt                      >  1:3346009/1‑91 (MQ=255)

                                                                      tACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAgtgt                      >  1:2425750/1‑91 (MQ=255)

                                                                       aCTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGtt                     >  1:2729254/1‑91 (MQ=255)

                                                                        cTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGtta                    <  1:1155722/91‑1 (MQ=255)

                                                                        cTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGtta                    >  1:1279625/1‑91 (MQ=255)

                                                                         tttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttat                   >  1:1975422/1‑91 (MQ=255)

                                                                         tttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttat                   <  1:2687039/91‑1 (MQ=255)

                                                                         tttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttat                   >  1:1551742/1‑91 (MQ=255)

                                                                         tttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttat                   >  1:1551829/1‑91 (MQ=255)

                                                                          ttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttatt                  >  1:12871/1‑91 (MQ=255)

                                                                          ttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttatt                  >  1:2594558/1‑91 (MQ=255)

                                                                          ttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttatt                  >  1:2656442/1‑91 (MQ=255)

                                                                          ttGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttatt                  >  1:2656457/1‑91 (MQ=255)

                                                                           tGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttatta                 >  1:891711/1‑91 (MQ=255)

                                                                           tGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGttatta                 <  1:2400871/91‑1 (MQ=255)

                                                                              aatctAAATTACCCGTGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAAtt              >  1:3771619/5‑91 (MQ=39)

                                                                                atTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTa            <  1:3998278/91‑1 (MQ=255)

                                                                                atTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTGCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTa            <  1:4003522/91‑1 (MQ=255)

                                                                                 ttAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTAt           <  1:294688/91‑1 (MQ=255)

                                                                                 ttAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTAt           <  1:3068244/91‑1 (MQ=255)

                                                                                  tAAATTACCCGTGTAATTTAGATTTATTCATATTAAGAATAGCTCTGCTTTCCTGTCTTTTCACTTTAACGACTAGTGTTATTAATGTATg          >  1:176467/1‑91 (MQ=255)

                                                                                   aaaTTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATgg         >  1:3844198/1‑91 (MQ=255)

                                                                                    aaTTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGa        >  1:547713/1‑91 (MQ=255)

                                                                                      ttACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACt      <  1:247496/91‑1 (MQ=255)

                                                                                       tACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACtt     >  1:2808647/1‑91 (MQ=255)

                                                                                         cacGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTc   >  1:1820317/1‑91 (MQ=255)

                                                                                          acGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCt  <  1:1793630/91‑1 (MQ=255)

                                                                                          acGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCt  >  1:1965903/1‑91 (MQ=255)

                                                                                          acGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCt  >  1:3330301/1‑91 (MQ=255)

GAAAATTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCT  >  USA300TCH1516_ALE/1038029‑1038209

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