BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	494,244	G→T	8.2%	R253R (CGC→CGA)	USA300HOU_RS02360 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	494,244	0	G	T	8.2%	224.5 / 2.1	74	R253R (CGC→CGA)	USA300HOU_RS02360	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (47/20); new base T (5/2); total (52/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAGGTGCTTTTTGTTATAGATCATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCAAATGA  >  USA300TCH1516_ALE/494104‑494384

aTAAAGGGGCTGTTTGTTATAGATCATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTg                                                                                                                                              <  2:1047154/141‑1 (MQ=255)

 aaaaTGTGATTTTTGTTATAGATCATCGAACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTgc                                                                                                                                             <  2:501248/140‑1 (MQ=255)

       gtcTTTTTGTTATAGAACATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTc                                                                                                                                       <  2:1049118/139‑1 (MQ=255)

            tgtgggataGATCATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGAt                                                                                                                                  >  2:1098602/7‑141 (MQ=255)

               gTTATAGATCATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGTGAGTTGCGATTCATTACCAGACTCATTATTTTTGCGCTGTCGAGATGGt                                                                                                                               >  1:418414/1‑141 (MQ=255)

                 tatgGATCATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCt                                                                                                                             <  2:128967/141‑1 (MQ=255)

                 tataGATCATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATACCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCt                                                                                                                             <  1:446432/141‑1 (MQ=255)

                 tataGAGTAACGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCt                                                                                                                             <  2:675527/141‑1 (MQ=255)

                         atcGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTttatta                                                                                                                     <  2:842836/141‑1 (MQ=255)

                          tcGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGGTGGTCTTTTTATTaa                                                                                                                    >  2:758032/1‑141 (MQ=255)

                           cGGACGATTTACTATAGTAAACAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTGGCCCCGTCGAGATGGTCTTTTTATTaaa                                                                                                                   >  2:774895/1‑141 (MQ=255)

                           cGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATTGTCTTTTTATTgaa                                                                                                                   >  1:967818/1‑141 (MQ=255)

                           cGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTaaa                                                                                                                   >  1:723615/1‑141 (MQ=255)

                             gACAATTTACTATAGTAACAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTCGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGGCTTCTTATTcaaac                                                                                                                 >  2:830020/1‑140 (MQ=255)

                             gACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTaaaaa                                                                                                                 >  1:1056704/1‑141 (MQ=255)

                             gACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTaaaaa                                                                                                                 <  2:12151/141‑1 (MQ=255)

                             gACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTaaaaa                                                                                                                 >  2:94349/1‑141 (MQ=255)

                             gACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTGTTGCGCTGTCGAGCTGGTCTTTTTATTaagaa                                                                                                                 >  1:737962/1‑141 (MQ=255)

                               cAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTACAAATg                                                                                                               >  1:993740/1‑141 (MQ=255)

                               cAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATg                                                                                                               >  1:436188/1‑141 (MQ=255)

                               cAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTAGAGATGGTCTTTTTATTAAAAATg                                                                                                               >  2:545453/1‑141 (MQ=255)

                                      cTATAGTAAAAAGCGACCTAGTGAACAATTGTCATATATCCGCAGGTCGCTTACCTTGAGTTATATTGCTAGTTGCGTTTAATTGTTGGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCGGTatt                                                                                                        >  2:103204/1‑139 (MQ=255)

                                      cTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAAt                                                                                                        <  1:364958/141‑1 (MQ=255)

                                            tAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGc                                                                                                  <  2:846486/141‑1 (MQ=255)

                                             aaaaaGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTCTTCCTGAAAAATGCCGTAATCCAAGcc                                                                                                 <  1:1098813/141‑1 (MQ=255)

                                              aaaaGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCg                                                                                                <  2:819240/141‑1 (MQ=255)

                                                aaGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTa                                                                                              <  1:901514/141‑1 (MQ=255)

                                                  gCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCA‑TTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATc                                                                                           >  1:637174/1‑141 (MQ=255)

                                                  gCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATAGGCCATCTTATCTTCTTTGCGCTGTCGAGATGGTCTTTTTATTCAAAATGCCGTAATCTAAGACGTact                                                                                            >  1:1088261/1‑139 (MQ=255)

                                                       cTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGaa                                                                                       >  1:424606/1‑141 (MQ=255)

                                                          gTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCTGTAATCGGAATAc                                                                                    >  1:11287/1‑141 (MQ=255)

                                                           tGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACt                                                                                   >  2:601413/1‑141 (MQ=255)

                                                           tGAACAATTGACATATATCCACAGGTCGCTTAACTGAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTACTCGGAGTAcg                                                                                   >  2:1017696/1‑140 (MQ=255)

                                                                aaTTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTg                                                                              <  1:601413/141‑1 (MQ=255)

                                                                 aTTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTTTTTTTTTTTTATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAAGCGGAATACTGATTGc                                                                             >  2:681038/1‑141 (MQ=255)

                                                                  tggACATAGATGCACAGGTCGCTTAACTTAAGTTTGAATGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGGTGGTCTTTTTATTAAAAATGCCGGAATCCAAGCCGTTATCGGAATACTGATTGCa                                                                            <  2:94333/139‑1 (MQ=255)

                                                                  atGACATATATCCACAGGTCGTTTAACTGAAGTGGTCTTTCGAGTGGCGGGGAATCGTGAGAGTGACCATTTTGTCGCTGTCTAGATGTTCTTTTTATTAAAGATGCCGTATTCCAAGCCTTGATCGTAATACTGATTGCa                                                                            <  2:118855/140‑1 (MQ=255)

                                                                   tGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCaa                                                                           >  2:1045943/1‑141 (MQ=255)

                                                                    gccATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAAc                                                                          <  2:798879/139‑1 (MQ=255)

                                                                     aCATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACg                                                                         >  1:870162/1‑141 (MQ=255)

                                                                     aCATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACg                                                                         >  1:868922/1‑141 (MQ=255)

                                                                                caGGTCGCTTAATTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGcc                                                              <  1:425597/141‑1 (MQ=255)

                                                                                caGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCACTACCGcc                                                              >  1:274366/1‑141 (MQ=255)

                                                                                  ggTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTCATCCAAGCCGTAATCGGAATACTGATTGCAACGGCCATACCGCCTa                                                            >  1:506684/1‑141 (MQ=255)

                                                                                    tCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTgaa                                                          >  1:420827/1‑141 (MQ=255)

                                                                                     cGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGGAACGGCAATACCGCCTaaaa                                                         >  1:622201/1‑141 (MQ=255)

                                                                                     cGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTaaaa                                                         >  1:897860/1‑141 (MQ=255)

                                                                                     cGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGGCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTaaaa                                                         >  2:336954/1‑141 (MQ=255)

                                                                                             ttaaCTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACTGCAATACCGCCTAAAATAATAGaa                                                 >  1:750549/1‑141 (MQ=255)

                                                                                              taaGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGaaa                                                <  1:305985/141‑1 (MQ=255)

                                                                                              taaGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGaaa                                                <  1:515825/141‑1 (MQ=255)

                                                                                              taaGTTATATTGCTAGTTCCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGaaa                                                <  1:1022644/141‑1 (MQ=255)

                                                                                               aaGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCTCCTAAAATAATAGaaat                                               >  1:684833/1‑141 (MQ=255)

                                                                                                     tatTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTGGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATGCTGATTGCAACGGCAATACCGGCTAAAGTAATAGaaaaaaatgc                                         >  1:76799/1‑139 (MQ=255)

                                                                                                        tGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATCCTGATTGCAACGACAATACCGCCTAAAATAATAGAAATAAATTCTTg                                      >  2:485188/1‑141 (MQ=255)

                                                                                                        tGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAAAAAATTCTTg                                      >  2:495269/1‑141 (MQ=255)

                                                                                                                gCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGAAGGTCTTTTTATTAAAACTGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCaaaga                              >  1:368601/1‑139 (MQ=255)

                                                                                                                 cGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAatat                             <  2:74343/141‑1 (MQ=255)

                                                                                                                 cGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAatat                             <  2:707611/141‑1 (MQ=255)

                                                                                                                  gATTAATTGATAGACTCATCATTTTTTTTTTTTTTTTATGATCTTTTTAGTAAAAATGCCGTAATCCAGGTCGTAATGGGAATACTGATTGCCGCGGCAGTACCGCCTAACATCATAAAAATAAATCCTTGCGCaattttt                            >  2:243949/1‑136 (MQ=255)

                                                                                                                  gATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATAtt                            >  2:82156/1‑141 (MQ=255)

                                                                                                                       aTTGATAGACTCATCATTTTTTTTTTTTTTTTTTAGTCTTCTTATTAGAAATGCCGTAATCCAAGCTGTAGTCGGAATACTTATTGCGACGTCAGTACCGCCTACAATAATAGTAATCAATTTTTTGTGAACTTTTTTCGa                       >  2:716492/1‑141 (MQ=255)

                                                                                                                       aTTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCAGAATACTTATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAACTCTTTGGCACATAtttacca                       >  2:172480/1‑137 (MQ=255)

                                                                                                                        ttGATAGACTCATCATTTTTGCGCTGTGGAGATGGTCTTTTTATTAAAAATGCCGTAATACAAGCGGTAATCGGAATACTGATTGCAACGGCAATACCGCCTATAATAATCGAAATAAATTCTTGGGCAAATATTTTCggg                      >  1:365276/1‑139 (MQ=255)

                                                                                                                         tGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGt                     >  1:781954/1‑141 (MQ=255)

                                                                                                                          gATAGACTCATCATTTTTTTTTTTTTTTTATGGTCTTTTTAGCAAAAATGACGTAATCCAAGCCGTGATCGGAATACTGATTGCAACTGCAATACCGCCTAAAATAATAGAAATACATTTTTTGGCTTTTATTTTCGAGtt                    >  2:12953/1‑141 (MQ=255)

                                                                                                                          gATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGtt                    >  2:651988/1‑141 (MQ=255)

                                                                                                                            tAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTta                  >  1:234344/1‑141 (MQ=255)

                                                                                                                               aCTCATCATTTTTTTTTTTTTTAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACGTATTGCATCGGCAATACCGCCTAAAATAATAGAAATAAATTTTTTGGCAAATTTTTTCGCGTTTataa               >  2:556542/1‑141 (MQ=255)

                                                                                                                                cTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTGTTGGGCAAATATGTCGGAGTTTataat              >  1:902731/1‑141 (MQ=255)

                                                                                                                                  catcatTTTGTCTCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAatat            <  1:1026234/141‑1 (MQ=255)

                                                                                                                                       tttttGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCa       >  1:784975/1‑141 (MQ=255)

                                                                                                                                       tttttGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCa       >  1:179971/1‑141 (MQ=255)

                                                                                                                                       tttttGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCa       >  1:117348/1‑141 (MQ=255)

                                                                                                                                        ttttGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCaa      >  2:881938/1‑141 (MQ=255)

                                                                                                                                          ttGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGTCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGAACAAAt    >  2:217435/1‑141 (MQ=255)

                                                                                                                                            gcgcTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTACTCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCAAATGa  >  1:969948/1‑141 (MQ=255)

                                                                                                                                            gcgcTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCAAATGa  >  1:934569/1‑141 (MQ=255)

                                                                                                                                            gcgcTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCAAATGa  >  1:916547/1‑141 (MQ=255)

                                                                                                                                            gcgcTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCAAATGa  >  1:755301/1‑141 (MQ=255)

                                                                                                                                            gcgcTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCAAATGa  >  1:271125/1‑141 (MQ=255)

ATAAAGGTGCTTTTTGTTATAGATCATCGGACAATTTACTATAGTAAAAAGCGACCTAGTGAACAATTGACATATATCCACAGGTCGCTTAACTTAAGTTATATTGCTAGTTGCGATTAATTGATAGACTCATCATTTTTGCGCTGTCGAGATGGTCTTTTTATTAAAAATGCCGTAATCCAAGCCGTAATCGGAATACTGATTGCAACGGCAATACCGCCTAAAATAATAGAAATAAATTCTTGGGCAAATATTTTCGAGTTTATAATATGACCAAATGA  >  USA300TCH1516_ALE/494104‑494384

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