BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	692,551	Δ1 bp	5.2%	coding (423/1092 nt)	USA300HOU_RS03340 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	692,551	0	T	.	5.2%	740.6 / 28.8	213	coding (423/1092 nt)	USA300HOU_RS03340	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (97/105); new base . (9/2); total (106/107)
Fisher's exact test for biased strand distribution p-value = 3.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 244 total aligned reads displayed.

CTTCTTCCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCAAAAAA  >  USA300TCH1516_ALE/692462‑692640

cttcttCCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATtaata                                                                                          >  2:3626820/1‑91 (MQ=255)

 ttcttcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATtaataa                                                                                         >  1:3584695/1‑91 (MQ=255)

  tcttcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATaaa                                                                                        <  1:998704/91‑1 (MQ=255)

   cttcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAg                                                                                       <  2:3169270/91‑1 (MQ=255)

   cttcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAg                                                                                       <  1:1420048/91‑1 (MQ=255)

    ttcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAgg                                                                                      <  2:830946/91‑1 (MQ=255)

    ttcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAgg                                                                                      <  1:2609868/91‑1 (MQ=255)

    ttcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAgg                                                                                      <  2:2068614/91‑1 (MQ=255)

     tcCGTCTGACAAATTACAACCGTTGAGTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGc                                                                                     <  1:2891079/91‑1 (MQ=255)

     tcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTTATTAATTAA‑AAAGGc                                                                                     <  1:1308089/91‑1 (MQ=255)

     tcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGc                                                                                     >  1:1225716/1‑91 (MQ=255)

     tcCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGc                                                                                     >  1:2593120/1‑91 (MQ=255)

      ccGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCa                                                                                    <  2:3040825/91‑1 (MQ=255)

         tCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACtt                                                                                 <  2:721679/91‑1 (MQ=255)

         tCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACtt                                                                                 >  1:505431/1‑91 (MQ=255)

             aCAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTaa                                                                             >  2:2284447/1‑91 (MQ=255)

             aCAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTaa                                                                             <  2:2597497/91‑1 (MQ=255)

             aCAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTaa                                                                             >  1:2270141/1‑91 (MQ=255)

               aaaTTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAat                                                                           <  1:3224248/91‑1 (MQ=255)

                  ttAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATg                                                                        >  1:59731/1‑91 (MQ=255)

                  ttAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATg                                                                        <  2:782009/91‑1 (MQ=255)

                   tAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGt                                                                       <  2:3330036/91‑1 (MQ=255)

                   tAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGt                                                                       >  1:3664231/1‑91 (MQ=255)

                     aaaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGttt                                                                     <  1:2760898/91‑1 (MQ=255)

                     aaaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGttt                                                                     <  1:316588/91‑1 (MQ=255)

                      aaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGtttt                                                                    >  1:1404619/1‑91 (MQ=255)

                      aaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGtttt                                                                    >  1:2196479/1‑91 (MQ=255)

                      aaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGtttt                                                                    >  1:569959/1‑91 (MQ=255)

                      aaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGtttt                                                                    >  1:767945/1‑91 (MQ=255)

                      aaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGtttt                                                                    >  1:1210166/1‑91 (MQ=255)

                      aaCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGtttt                                                                    >  2:1057563/1‑91 (MQ=255)

                       aCCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTta                                                                   >  1:2302162/1‑91 (MQ=255)

                        ccGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTtat                                                                  >  1:3771051/1‑91 (MQ=255)

                        ccGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTtat                                                                  <  1:2431105/91‑1 (MQ=255)

                          gTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATaa                                                                <  2:3589659/91‑1 (MQ=255)

                          gTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATaa                                                                <  2:2808736/91‑1 (MQ=255)

                          gTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATaa                                                                >  1:3387768/1‑91 (MQ=255)

                            ttATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAatc                                                              <  2:689033/91‑1 (MQ=255)

                              aTTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAatcat                                                            <  2:1817169/91‑1 (MQ=255)

                              aTTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAatcat                                                            <  1:1579492/91‑1 (MQ=255)

                               ttGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAatcatc                                                           >  2:231590/1‑91 (MQ=255)

                                tGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAAtcatca                                                          >  2:481074/1‑91 (MQ=255)

                                 gcgcTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGTCACTTTTAAATATGTTTTATAATCATCaa                                                         <  2:3387768/91‑1 (MQ=255)

                                 gcgcTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCaa                                                         >  1:1284921/1‑91 (MQ=255)

                                  cgcTTCCAGCTCTAAAACGTGAAACACACTCCGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAAt                                                        >  2:128635/1‑91 (MQ=255)

                                   gcTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAAtt                                                       >  1:1418537/1‑91 (MQ=255)

                                        cAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAat                                                  >  2:2235262/1‑91 (MQ=255)

                                         aGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAata                                                 <  1:3600533/91‑1 (MQ=255)

                                          gCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAatat                                                <  2:526704/91‑1 (MQ=255)

                                          gCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAatat                                                <  2:2710250/91‑1 (MQ=255)

                                           ctctAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATc                                               >  1:2667037/1‑91 (MQ=255)

                                           ctctAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATc                                               >  2:2167017/1‑91 (MQ=255)

                                           ctctAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATc                                               >  1:1850100/1‑91 (MQ=255)

                                            tctAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCa                                              >  1:2173938/1‑91 (MQ=255)

                                            tctAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCa                                              >  2:146342/1‑91 (MQ=255)

                                              tAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAAt                                            <  1:1712484/91‑1 (MQ=255)

                                                aaaCGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTc                                          <  2:3881488/91‑1 (MQ=255)

                                                 aaCGTGAAACACACCCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCt                                         >  2:2994830/1‑91 (MQ=255)

                                                  aCGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCtt                                        <  1:2063152/91‑1 (MQ=255)

                                                   cGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTg                                       <  1:1147826/91‑1 (MQ=255)

                                                   cGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTg                                       <  1:3622954/91‑1 (MQ=255)

                                                    gTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGt                                      <  1:1262765/91‑1 (MQ=255)

                                                    gTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGt                                      >  1:148344/1‑91 (MQ=255)

                                                    gTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGt                                      >  1:2351499/1‑91 (MQ=255)

                                                      gAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTaa                                    >  1:3373410/1‑91 (MQ=255)

                                                      gAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTaa                                    >  2:1505247/1‑91 (MQ=255)

                                                       taaCACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTaaa                                   >  2:2002263/2‑91 (MQ=255)

                                                       aaaCACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTaaa                                   >  1:3832511/1‑91 (MQ=255)

                                                        aaCACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTaaaa                                  >  1:2213114/1‑91 (MQ=255)

                                                         acacacTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAAt                                 >  1:3695311/1‑91 (MQ=255)

                                                         acacacTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTCTATAATCATCAATTGCAATATCAATTCTTGTAAAAt                                 >  2:3057676/1‑91 (MQ=255)

                                                          cacacTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAAtt                                <  2:3173719/91‑1 (MQ=255)

                                                          cacacTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAAtt                                <  2:2907309/91‑1 (MQ=255)

                                                           acacTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAAttt                               >  1:2476705/1‑91 (MQ=255)

                                                           acacTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAACATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAAttt                               >  1:3072294/1‑91 (MQ=255)

                                                             acTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGc                             <  1:2029604/91‑1 (MQ=255)

                                                             acTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGc                             >  2:2098109/1‑91 (MQ=255)

                                                               tCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCAt                           <  2:1662407/91‑1 (MQ=255)

                                                               tCAGCTTATTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCAt                           <  2:3442497/91‑1 (MQ=255)

                                                                  gCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTaa                        >  2:2486179/1‑91 (MQ=255)

                                                                  gCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTaa                        <  2:1927204/91‑1 (MQ=255)

                                                                      tttCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGa                    <  1:729099/91‑1 (MQ=255)

                                                                      tttCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGa                    <  2:1114803/91‑1 (MQ=255)

                                                                      tttCAGCCTTT‑TCAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGa                    >  1:2969828/1‑91 (MQ=255)

                                                                       ttCAGCCTTT‑TTTATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCATTTGCAATATCAATTCGTGTAAAATTTGCATTAACTGAt                   <  2:3747844/91‑1 (MQ=255)

                                                                         gtgCCTTTATTAATTAA‑AAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGAtga                 >  1:66074/3‑89 (MQ=255)

                                                                         gtgCCTTTATTAATTAA‑AAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATaa                 >  1:2622321/3‑91 (MQ=255)

                                                                         cAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATaa                 <  2:1171158/91‑1 (MQ=255)

                                                                           gCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAAc               >  2:352521/1‑91 (MQ=255)

                                                                              tttattaattaa‑AAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACAtt            >  1:3439154/5‑91 (MQ=255)

                                                                                  ttaattaaTAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGt         <  1:977792/91‑1 (MQ=255)

                                                                                   taattaatAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTc        >  1:1115974/1‑91 (MQ=255)

                                                                                    aattaatAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCa       >  2:23347/1‑91 (MQ=255)

                                                                                    aattaatAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCa       >  2:1795068/1‑91 (MQ=255)

                                                                                     attaatAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATa                  >  2:1359596/1‑79 (MQ=255)

                                                                                       taataaAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCaaaa    <  2:3795381/91‑1 (MQ=255)

                                                                                       taataaAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCaaaa    <  2:2393249/91‑1 (MQ=255)

                                                                                       taataaAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTAGTAaaaa    <  2:2438487/91‑1 (MQ=255)

                                                                                        aataaAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCaaaaa   >  2:1238740/1‑91 (MQ=255)

                                                                                         ataaAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCaaaaaa  >  2:3161979/1‑91 (MQ=255)

CTTCTTCCGTCTGACAAATTAAAACCGTTTATTGCGCTTCCAGCTCTAAAACGTGAAACACACTCAGCTTTTTCAGCCTTT‑TTAATTAATAAAGGCACTTTTAAATATGTTTTATAATCATCAATTGCAATATCAATTCTTGTAAAATTTGCATTAACTGATAAACATTTGTCAAAAAA  >  USA300TCH1516_ALE/692462‑692640

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