BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,129,372	A→G	5.4%	intergenic (+205/+33)	hlb_2 → / ← USA300HOU_RS10870	Phospholipase C/putative leukocidin‑like protein 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,129,372	0	A	G	5.4%	363.1 / 6.6	110	intergenic (+205/+33)	hlb_2/USA300HOU_RS10870	Phospholipase C/putative leukocidin‑like protein 1
Reads supporting (aligned to +/- strand): ref base A (54/50); new base G (2/4); total (56/54)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 123 total aligned reads displayed.

TAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCC  >  USA300TCH1516_ALE/2129282‑2129460

tAATACAAGTCTTTACTCTACACTCGAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACcta                                                                                          <  1:3850465/91‑1 (MQ=255)

tAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACcta                                                                                          <  1:623623/91‑1 (MQ=255)

 aaTACAAGTCTATACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctac                                                                                         <  1:3760613/91‑1 (MQ=255)

  aTACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctact                                                                                        >  1:1095848/1‑91 (MQ=255)

  aTACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctact                                                                                        <  1:3886468/91‑1 (MQ=255)

  aTACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctact                                                                                        <  1:3404303/91‑1 (MQ=255)

   tACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctacta                                                                                       >  1:3315570/1‑91 (MQ=255)

   tACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctacta                                                                                       <  1:758654/91‑1 (MQ=255)

   tACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctacta                                                                                       >  1:451132/1‑91 (MQ=255)

    aCAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCtactac                                                                                      <  1:1033574/91‑1 (MQ=255)

      aaGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACAt                                                                                    <  1:1561448/91‑1 (MQ=255)

      aaGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACAt                                                                                    <  1:3408763/91‑1 (MQ=255)

                 ctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                         >  1:356764/1‑91 (MQ=255)

                 ctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                         <  1:771787/91‑1 (MQ=255)

                 ctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                         >  1:193177/1‑91 (MQ=255)

                   acacTCATACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAgg                                                                       <  1:2069206/91‑1 (MQ=255)

                   acacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAgg                                                                       <  1:2586562/91‑1 (MQ=255)

                   acacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAgg                                                                       <  1:2425699/91‑1 (MQ=255)

                    cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                      >  1:3505322/1‑91 (MQ=255)

                    cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                      >  1:955271/1‑91 (MQ=255)

                    cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                      >  1:1295139/1‑91 (MQ=255)

                     acTCAGACAAGATTCATACACTGCACGTCATAATAAATCGATCTATTCAAATAGAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                     <  1:3860216/91‑1 (MQ=255)

                     acTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                     <  1:1635355/91‑1 (MQ=255)

                     acTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                     <  1:1483908/91‑1 (MQ=255)

                       tCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAAc                                                                   <  1:1458351/91‑1 (MQ=21)

                       tCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAAc                                                                   <  1:1399893/91‑1 (MQ=255)

                       ccAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAAc                                                                   <  1:2167243/90‑1 (MQ=255)

                        cAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACt                                                                  >  1:633441/1‑91 (MQ=255)

                           aCAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACtttt                                                               >  1:3549304/1‑91 (MQ=255)

                             aaGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTttat                                                             <  1:289507/91‑1 (MQ=255)

                                aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTttattta                                                          <  1:1269915/91‑3 (MQ=17)

                                 ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTttatttat                                                         <  1:3883448/91‑4 (MQ=17)

                                 ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACt                                                         <  1:3355806/91‑1 (MQ=255)

                                 ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACt                                                         <  1:3412227/91‑1 (MQ=255)

                                 ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACt                                                         <  1:2547868/91‑1 (MQ=255)

                                  tCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACtt                                                        <  1:1981610/91‑1 (MQ=255)

                                  tCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACtt                                                        <  1:944800/91‑1 (MQ=255)

                                  tCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACtt                                                        <  1:944595/91‑1 (MQ=255)

                                    aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACttct                                                      >  1:2444349/1‑89 (MQ=255)

                                    aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                      >  1:2444166/1‑91 (MQ=255)

                                    aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                      >  1:3842783/1‑91 (MQ=255)

                                    aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                      >  1:3241489/1‑91 (MQ=255)

                                    aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTATATGTAGCAGGCAACTTTTATTACTTAt                                                      >  1:3241251/1‑91 (MQ=255)

                                     tACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAtt                                                     >  1:4054459/1‑91 (MQ=255)

                                       cacTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTc                                                   <  1:3105643/91‑1 (MQ=255)

                                            cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                              <  1:1952613/91‑1 (MQ=255)

                                            cACGTCATAATAAATCTATCTATGCAAATATAAATAAAAGTTACCTACTACATGCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                              <  1:1800818/91‑1 (MQ=255)

                                              cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                            <  1:2638457/91‑1 (MQ=255)

                                              cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                            >  1:716604/1‑91 (MQ=255)

                                              cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                            >  1:3157006/1‑91 (MQ=255)

                                               gTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcatt                                           >  1:1152807/1‑91 (MQ=255)

                                                  ataataAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAANTTTNATTANTTATTTCTTTtcattatc                                        <  1:3501616/91‑1 (MQ=255)

                                                  ataataAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatc                                        <  1:3500987/91‑1 (MQ=255)

                                                     ataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatt                                     <  1:1111367/91‑1 (MQ=255)

                                                      taaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatta                                    <  1:1298463/91‑1 (MQ=255)

                                                      taaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatta                                    <  1:3604504/91‑1 (MQ=255)

                                                        aaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAg                                  >  1:719600/1‑91 (MQ=255)

                                                         atctatctATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGt                                 >  1:929615/1‑91 (MQ=255)

                                                          tctatctaTTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTa                                <  1:635771/91‑1 (MQ=255)

                                                          tctatctaTTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTa                                >  1:1293146/1‑91 (MQ=255)

                                                             atctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtt                             <  1:2866113/91‑1 (MQ=255)

                                                                 atTCAAATCTAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTAc                         >  1:637764/1‑91 (MQ=255)

                                                                   tCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTaaaa                       >  1:3743214/1‑91 (MQ=255)

                                                                    cAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACaaa                      >  1:4105823/1‑91 (MQ=255)

                                                                      aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCGTTATCATTAAGTACTTTTACAAACt                    >  1:3860751/1‑91 (MQ=255)

                                                                      aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACt                    >  1:789354/1‑91 (MQ=255)

                                                                      aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACt                    <  1:178691/91‑1 (MQ=255)

                                                                      aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACt                    <  1:178745/91‑1 (MQ=255)

                                                                        tataAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTc                  <  1:2105378/91‑1 (MQ=255)

                                                                          taaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac                >  1:770211/1‑91 (MQ=255)

                                                                          taaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac                >  1:769587/1‑91 (MQ=255)

                                                                          gtaataaaAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac                >  1:1751998/3‑91 (MQ=17)

                                                                          gtaataaaAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac                >  1:1777754/3‑91 (MQ=17)

                                                                           aaataaaAGTTGCCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcaca               >  1:4050860/1‑91 (MQ=255)

                                                                           aaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcaca               >  1:2913336/1‑91 (MQ=255)

                                                                            aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt              >  1:1320985/1‑91 (MQ=255)

                                                                            aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt              >  1:1824287/1‑91 (MQ=255)

                                                                            aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt              >  1:1283651/1‑91 (MQ=255)

                                                                            aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt              >  1:1285919/1‑91 (MQ=255)

                                                                            aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt              >  1:1287206/1‑91 (MQ=255)

                                                                             ataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAtt             >  1:1270633/1‑91 (MQ=255)

                                                                             ataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAtt             >  1:1270610/1‑91 (MQ=255)

                                                                             ataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAtt             <  1:2690648/91‑1 (MQ=255)

                                                                              taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa            >  1:76410/1‑91 (MQ=255)

                                                                              taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa            >  1:3341058/1‑91 (MQ=255)

                                                                              taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa            >  1:1115612/1‑91 (MQ=255)

                                                                              taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa            >  1:1329836/1‑91 (MQ=255)

                                                                              taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa            <  1:3800225/91‑1 (MQ=255)

                                                                                aaaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtg          >  1:3385513/1‑91 (MQ=255)

                                                                                  aGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACACACTTCACATTAtgtg        <  1:3546681/91‑1 (MQ=255)

                                                                                  aGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgtg        >  1:3602691/1‑91 (MQ=255)

                                                                                   gTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATgtgt       >  1:2372795/1‑91 (MQ=255)

                                                                                   gTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATgtgt       <  1:2536104/91‑1 (MQ=255)

                                                                                   gTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATgtgt       <  1:3645616/91‑1 (MQ=255)

                                                                                     tACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCt     <  1:3129438/91‑1 (MQ=255)

                                                                                      aCCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCtt    >  1:714540/1‑91 (MQ=255)

                                                                                       ccTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTc   >  1:3254587/1‑91 (MQ=255)

                                                                                        ctactaCATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTcc  >  1:1973215/1‑91 (MQ=255)

                                                                                        ctactaCATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTcc  <  1:2285541/91‑1 (MQ=255)

                                                                                        ctactaCATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTcc  >  1:2908778/1‑91 (MQ=255)

TAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCC  >  USA300TCH1516_ALE/2129282‑2129460

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