BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,559,831	G→A	6.3%	intergenic (+25/+98)	bcr_2 → / ← USA300HOU_RS13190	Bicyclomycin resistance protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,559,831	0	G	A	6.3%	404.2 / 8.0	120	intergenic (+25/+98)	bcr_2/USA300HOU_RS13190	Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (61/51); new base A (4/4); total (65/55)
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

Only 100 of 123 total aligned reads displayed.

GCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTT  >  USA300TCH1516_ALE/2559741‑2559920

gCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAg                                                                                           >  1:3266496/1‑91 (MQ=255)

gCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAg                                                                                           >  1:3493893/1‑91 (MQ=255)

gCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAg                                                                                           >  1:422572/1‑91 (MQ=255)

     tCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAg                                                                                      <  1:2930022/91‑1 (MQ=255)

       tATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGat                                                                                    >  1:2920730/1‑91 (MQ=255)

       tATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGat                                                                                    <  1:741146/91‑1 (MQ=255)

        aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                   >  1:3936099/1‑91 (MQ=255)

        aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                   >  1:255924/1‑91 (MQ=255)

        aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                   >  1:256231/1‑91 (MQ=255)

        aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                   >  1:977515/1‑91 (MQ=255)

         ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                  >  1:166407/1‑91 (MQ=255)

            gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                               <  1:1903672/91‑1 (MQ=255)

             tcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctat                                                                              >  1:1375195/1‑91 (MQ=255)

               agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATatctatct                                                                            <  1:1237257/91‑1 (MQ=255)

                gtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCtt                                                                           >  1:1308850/1‑91 (MQ=255)

                  cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                         <  1:109030/91‑1 (MQ=38)

                  cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                         <  1:2852803/91‑1 (MQ=38)

                  cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTg                                                                         >  1:903506/1‑91 (MQ=255)

                    aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                       <  1:2595816/91‑1 (MQ=255)

                    aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                       <  1:2595842/91‑1 (MQ=255)

                     cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                      >  1:2588901/1‑91 (MQ=255)

                     cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                      >  1:2855206/1‑91 (MQ=255)

                      aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgctag                                                                     <  1:1880889/91‑5 (MQ=25)

                        aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                   <  1:1717337/91‑1 (MQ=255)

                        aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                   <  1:3731884/91‑1 (MQ=255)

                         tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCGTAACACTTCAACACAATTAGAACCCTAACAAAGATATCTATCTTTGTCAGGGt                                                                  >  1:3862826/1‑91 (MQ=255)

                         tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                  <  1:1701036/91‑1 (MQ=255)

                         tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                  <  1:2317717/91‑1 (MQ=255)

                          cATTTACTTTAAAATGATTAAAAAGAAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGtt                                                                 <  1:727817/91‑1 (MQ=255)

                            tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttct                                                               <  1:1873387/91‑1 (MQ=255)

                             ttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttctt                                                              >  1:1330796/1‑91 (MQ=255)

                              tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                             >  1:1400304/1‑91 (MQ=255)

                               aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                            >  1:379714/1‑91 (MQ=255)

                               aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                            >  1:1622759/1‑91 (MQ=255)

                                cTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTcttctt                                                           >  1:3763043/1‑91 (MQ=255)

                                 tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                          <  1:2537752/91‑1 (MQ=255)

                                  ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                         <  1:53180/91‑1 (MQ=255)

                                  ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                         <  1:53181/91‑1 (MQ=255)

                                    aaaaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATg                                                       >  1:169639/1‑91 (MQ=255)

                                    aaaaTGAGTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATg                                                       <  1:2127990/91‑1 (MQ=255)

                                       aTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAAt                                                    <  1:3666430/91‑1 (MQ=255)

                                          aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                 >  1:662875/1‑91 (MQ=255)

                                          aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                 >  1:248902/1‑91 (MQ=255)

                                          aTTAAAAAGCAACATGTAGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                 <  1:2186037/91‑1 (MQ=255)

                                           ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                <  1:3767537/91‑1 (MQ=255)

                                           ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                <  1:2020683/91‑1 (MQ=255)

                                            tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                               <  1:1615854/91‑1 (MQ=255)

                                            tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                               <  1:1615785/91‑1 (MQ=255)

                                            tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                               <  1:1745605/91‑1 (MQ=255)

                                              aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                             >  1:646651/1‑91 (MQ=255)

                                              aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                             >  1:979312/1‑91 (MQ=255)

                                              aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                             >  1:1616026/1‑91 (MQ=255)

                                                 aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                          <  1:1053193/91‑1 (MQ=255)

                                                 aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                          <  1:1838645/91‑1 (MQ=255)

                                                 aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                          <  1:3716348/91‑1 (MQ=255)

                                                  gCAACATGTCGCATAACACTTCAACATAATTAGTACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                         <  1:3396984/91‑1 (MQ=255)

                                                  gCAACATGTCGCATAACACTTCAACATAATTAGCACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                         <  1:3397144/91‑1 (MQ=255)

                                                  gCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                         >  1:232402/1‑91 (MQ=255)

                                                   cAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGaa                                        >  1:3209330/1‑91 (MQ=255)

                                                   cAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGaa                                        >  1:2160799/1‑91 (MQ=255)

                                                    aaCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAAt                                       >  1:94419/1‑91 (MQ=255)

                                                     aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                      >  1:133210/1‑91 (MQ=255)

                                                     aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                      >  1:1898618/1‑91 (MQ=255)

                                                     aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                      >  1:711171/1‑91 (MQ=255)

                                                     aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                      >  1:732015/1‑91 (MQ=255)

                                                       aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                    >  1:2465150/1‑91 (MQ=255)

                                                       aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                    >  1:347502/1‑91 (MQ=255)

                                                        tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                   >  1:3501701/1‑91 (MQ=255)

                                                        tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                   >  1:3446748/1‑91 (MQ=255)

                                                         gTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCa                                  >  1:324363/1‑91 (MQ=255)

                                                           cGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                <  1:1694184/91‑1 (MQ=255)

                                                            gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                               >  1:4003877/1‑91 (MQ=255)

                                                              aTAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaa                             <  1:1995811/91‑1 (MQ=255)

                                                               tAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaa                            >  1:1707978/1‑91 (MQ=255)

                                                                aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                           <  1:4048773/91‑1 (MQ=255)

                                                                  cacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAtt                         >  1:853930/1‑91 (MQ=255)

                                                                   acTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                        <  1:115746/91‑1 (MQ=255)

                                                                    cTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTAc                       >  1:45465/1‑91 (MQ=255)

                                                                     ttCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACg                      >  1:1736496/1‑91 (MQ=255)

                                                                       cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                    <  1:2176415/91‑1 (MQ=255)

                                                                        aaCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCt                   >  1:4034131/1‑91 (MQ=255)

                                                                         aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                  <  1:758898/91‑1 (MQ=255)

                                                                         aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                  <  1:758979/91‑1 (MQ=255)

                                                                            tAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCAt               <  1:341020/91‑1 (MQ=255)

                                                                            tAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCAt               >  1:3347304/1‑91 (MQ=255)

                                                                              aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg             <  1:2790018/91‑1 (MQ=255)

                                                                              aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg             <  1:4027266/91‑1 (MQ=255)

                                                                                tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa           <  1:1473803/91‑1 (MQ=255)

                                                                                aaGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa           >  1:4070440/2‑91 (MQ=21)

                                                                                 aGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa          >  1:2864163/1‑91 (MQ=255)

                                                                                  gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg         >  1:2573269/1‑91 (MQ=255)

                                                                                    aCCCTGACAAAGGTAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa       >  1:3184522/1‑91 (MQ=25)

                                                                                    aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa       >  1:54458/1‑91 (MQ=255)

                                                                                    aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa       >  1:489315/1‑91 (MQ=255)

                                                                                    aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa       >  1:3885314/1‑91 (MQ=255)

                                                                                    aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa       <  1:2576115/91‑1 (MQ=255)

                                                                                     cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc      >  1:3221577/1‑91 (MQ=255)

                                                                                      ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa     >  1:1974007/1‑91 (MQ=255)

                                                                                      ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa     <  1:2230861/91‑1 (MQ=255)

                                                                                         aGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt  <  1:3734147/91‑1 (MQ=255)

GCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTT  >  USA300TCH1516_ALE/2559741‑2559920

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