BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,038,122	G→T	6.0%	G237C (GGT→TGT)	USA300TCH1516_00964 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,038,122	0	G	T	6.0%	471.0 / 15.5	132	G237C (GGT→TGT)	USA300TCH1516_00964	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (70/54); new base T (3/5); total (73/59)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

AATTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCTTAA  >  USA300TCH1516_ALE/1038032‑1038212

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

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

 aTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACggg                                                                                           >  1:2394608/1‑91 (MQ=255)

 aTTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACggg                                                                                           <  1:1844071/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)

    tttCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACCCGTGTaa                                                                                        <  1:2630765/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)

      tCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAGATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAAtt                                                                                      >  1:3591740/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)

           gAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGAt                                                                                 <  1:785046/91‑1 (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)

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

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

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

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

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

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

                          ggTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGa                                                                  >  1:571930/1‑91 (MQ=255)

                           gTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGaa                                                                 >  1:764791/1‑91 (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:3121800/1‑91 (MQ=255)

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

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

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

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

                                      cATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCtt                                                      >  1:1908865/1‑91 (MQ=255)

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

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

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

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

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

                                                aTTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCttt                                            <  1:1659092/91‑1 (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:1908129/1‑91 (MQ=255)

                                                   cTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCa                                         <  1:3366924/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                            aCAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGa                                <  1:2178090/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                               tAAATTACCCGTGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATg             >  1:176767/1‑91 (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:3330301/1‑91 (MQ=255)

                                                                                        cGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCtt    <  1:210284/91‑1 (MQ=255)

                                                                                        cGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTG‑TATTAATTTATGGACTTTCTTa   >  1:2717686/1‑91 (MQ=255)

                                                                                         gggTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCTTa   <  1:741610/91‑1 (MQ=255)

                                                                                          ggTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCTTaa  >  1:708639/1‑91 (MQ=255)

                                                                                          ggTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCTTaa  >  1:2106612/1‑91 (MQ=255)

AATTTTTCAAAGAAATGCACGTATTGGGTAACGTTAATCATAACTCAAATTCTTTTAACTACAAGTGTACTTTGGATATTAAATTACACGGGTAATTTAGATTTATTCATATTAAGAATAGTTCTGCTTTCATGTCTTTTCATTTTAACGATTAGTGTTATTAATTTATGGACTTTCTTAA  >  USA300TCH1516_ALE/1038032‑1038212

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