BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,231,527	A→T	6.8%	F17I (TTT→ATT)	USA300HOU_RS11415 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,231,527	0	A	T	6.8%	594.4 / 22.8	176	F17I (TTT→ATT)	USA300HOU_RS11415	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (81/83); new base T (9/3); total (90/86)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 192 total aligned reads displayed.

CTGCAACTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAATATAGCAAA  >  USA300TCH1516_ALE/2231437‑2231617

cTGCAACTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTaaa                                                                                            >  2:2142877/1‑91 (MQ=255)

cTGCAACTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTAATTaaa                                                                                            <  2:3114629/91‑1 (MQ=255)

   cAACTTCCCTCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACCTCCAATTAAACATATATCTCCTTTTTTaataaa                                                                                         <  2:1516520/91‑1 (MQ=255)

   cAACTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTaaaaaa                                                                                         <  1:1156721/91‑1 (MQ=255)

   cAACTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTaaaaaa                                                                                         >  2:2605115/1‑91 (MQ=255)

    aaCTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTTTTAATAAAg                                                                                        <  2:2971788/91‑1 (MQ=255)

    aaCTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAg                                                                                        >  1:1458027/1‑91 (MQ=255)

     aCTTCCATCGATCCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAgg                                                                                       <  2:2969166/91‑1 (MQ=255)

     aCTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATGTCCTTTATTAAAAAAgg                                                                                       >  1:3325733/1‑91 (MQ=255)

     aCTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAgg                                                                                       <  1:1067744/91‑1 (MQ=255)

       ttCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGAc                                                                                     <  2:2655181/91‑1 (MQ=255)

       ttCCATCGAACCGACAGTCCCTATGTTTTCGCCAATGACTTCCGATGTTGGACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGAc                                                                                     <  2:3602554/91‑1 (MQ=255)

         ccATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACat                                                                                   <  1:704263/91‑1 (MQ=255)

          cATCGAACCGACAGTCCCTATTTTTTCTCCTATTACTTCCGTTGTTGAACATCCAATTAAACCTATATCTCCTTTATTAAAAAAGGACata                                                                                  >  2:1250437/1‑91 (MQ=255)

          cATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACata                                                                                  >  1:291552/1‑91 (MQ=255)

          cATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACata                                                                                  >  1:295727/1‑91 (MQ=255)

            tCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATc                                                                                >  2:1948365/1‑91 (MQ=255)

              gAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCtt                                                                              <  1:1699790/91‑1 (MQ=255)

               aaCCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCttt                                                                             >  1:3717719/1‑91 (MQ=255)

                aCCGACAGTCCCTATTTTTTTTCCTATTACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTTAAAAAGGACATATCtttt                                                                            >  2:726341/1‑91 (MQ=255)

                aCCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCtttt                                                                            <  2:3868230/91‑1 (MQ=255)

                  cGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTaa                                                                          >  1:1796862/1‑91 (MQ=255)

                  cGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTaa                                                                          <  1:3854989/91‑1 (MQ=255)

                   gACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTGAt                                                                         >  2:3220309/1‑91 (MQ=255)

                   gACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAAt                                                                         >  1:1558385/1‑91 (MQ=255)

                   gACAGTCCCTATTTTTTCGCCAATGACTTCCGAGGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAAt                                                                         <  1:3173574/91‑1 (MQ=255)

                     cAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTc                                                                       <  1:873430/91‑1 (MQ=255)

                       gTCCCTATTTTTTCTCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGt                                                                     >  2:2012458/1‑91 (MQ=255)

                       gTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTTTTAAAAAAGGACATATCTTTTAATTCGt                                                                     >  2:3775290/1‑91 (MQ=255)

                       gTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGt                                                                     >  2:3187262/1‑91 (MQ=255)

                        tCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTc                                                                    >  1:2668727/1‑91 (MQ=255)

                        tCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTc                                                                    <  2:840803/91‑1 (MQ=255)

                         cccTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCt                                                                   <  2:3299962/91‑1 (MQ=255)

                         cccTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCt                                                                   <  2:695982/91‑1 (MQ=255)

                           cTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCtaa                                                                 <  2:1060873/91‑1 (MQ=255)

                             aTTTTTTTGCTAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCtaata                                                               >  2:2218874/1‑91 (MQ=255)

                             aTTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCtaata                                                               >  2:182867/1‑91 (MQ=255)

                                  ttCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACAtt                                                          <  2:623869/91‑1 (MQ=255)

                                   tCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTg                                                         >  1:2837672/1‑91 (MQ=255)

                                   tCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTg                                                         >  1:2746457/1‑91 (MQ=255)

                                   tCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTg                                                         <  2:3454119/91‑1 (MQ=255)

                                     gCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTc                                                       >  2:2382148/1‑91 (MQ=255)

                                      ccAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCa                                                      <  2:163343/91‑1 (MQ=255)

                                      ccAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCa                                                      <  1:643972/91‑1 (MQ=255)

                                      ccAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCa                                                      >  1:3393152/1‑91 (MQ=255)

                                        aaTGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCaaa                                                    >  1:3473625/1‑91 (MQ=255)

                                          tGACCTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATc                                                  >  2:586473/1‑91 (MQ=255)

                                               tCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCa                                             <  1:1897587/91‑1 (MQ=255)

                                               tCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCa                                             >  1:558245/1‑91 (MQ=255)

                                                 cGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATa                                           >  1:3276865/1‑91 (MQ=255)

                                                  gATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATaa                                          >  2:187694/1‑91 (MQ=255)

                                                  gATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATaa                                          <  1:406886/91‑1 (MQ=255)

                                                      ttGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACGTATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGc                                      <  1:3074951/91‑1 (MQ=255)

                                                      ttGAACATCCAATTAAACATATATCTCCTTTATT‑AAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGcc                                     >  2:2565884/1‑91 (MQ=255)

                                                       tGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGcc                                     <  1:754811/91‑1 (MQ=255)

                                                       tGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGcc                                     <  1:753209/91‑1 (MQ=255)

                                                       tGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGcc                                     <  2:2357087/91‑1 (MQ=255)

                                                         aaCATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCa                                   >  1:3230159/1‑91 (MQ=255)

                                                          aCATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACGATGTCAAATCGTTCATAAAAGCCCAc                                  >  2:3355565/1‑91 (MQ=255)

                                                          aCATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCAc                                  >  1:3129989/1‑91 (MQ=255)

                                                          aCATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCAc                                  <  2:1051464/91‑1 (MQ=255)

                                                          aCATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCAc                                  >  1:3129403/1‑91 (MQ=255)

                                                            aTCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCAccc                                <  1:3484267/91‑1 (MQ=255)

                                                             tCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCt                               <  1:2559687/91‑1 (MQ=255)

                                                              ccAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCtt                              <  1:162753/91‑1 (MQ=255)

                                                               cAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTc                             >  2:2323515/1‑91 (MQ=255)

                                                                 aTTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCt                           >  2:3503403/1‑91 (MQ=255)

                                                                 aTTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCt                           >  2:1405095/1‑91 (MQ=255)

                                                                 aTTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCt                           >  2:2466839/1‑91 (MQ=255)

                                                                  ttAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCta                          >  2:1210616/1‑91 (MQ=255)

                                                                  ttAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCGTCCta                          >  2:3139296/1‑91 (MQ=255)

                                                                   tAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAAT‑ACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaa                        >  2:2873512/1‑91 (MQ=255)

                                                                    aaaCATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaa                        <  1:896610/91‑1 (MQ=255)

                                                                     aaCATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATaaa                                         <  1:1070537/73‑1 (MQ=255)

                                                                     aaCATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATaaa                                         >  2:1070537/1‑73 (MQ=255)

                                                                     aaCATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaa                       >  2:1175198/1‑91 (MQ=255)

                                                                     aaCATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaa                       <  1:3430128/91‑1 (MQ=255)

                                                                     aaCATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaa                       <  2:2176138/91‑1 (MQ=255)

                                                                      aCATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCACATCTTTCATAAAAGCCCACCCTTCCtaaaaa                      >  2:3807737/1‑91 (MQ=255)

                                                                       aagatatGTCCTTTTTTAATAAAGGAGATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCTACCCTTCCtaaaaat                     >  2:3389305/4‑91 (MQ=25)

                                                                         gatatGTCCTTTTTTAATAAAGGAGATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaaataa                   >  1:115856/2‑91 (MQ=37)

                                                                         gatatGTCCTTTTTTAATAAAGGAGATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaaataa                   >  2:1786555/2‑91 (MQ=37)

                                                                          atatCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaaataaa                  <  2:1182785/91‑1 (MQ=255)

                                                                           tatgtcCTTTTTTAATAAAGGAGATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaaataaaa                 >  2:538234/5‑91 (MQ=38)

                                                                           tatCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaaataaaa                 <  1:223487/91‑1 (MQ=255)

                                                                            atgtcCTTTTTTAATAAAGGAGATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaaataaaaa                <  1:2377149/88‑1 (MQ=38)

                                                                            atCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCtaaaaataaaaa                <  1:1405095/91‑1 (MQ=255)

                                                                             tctcCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCATCCTTCCTAAAAATaaaaaa               <  2:3070227/91‑1 (MQ=255)

                                                                              ctcCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAg              <  1:3374380/91‑1 (MQ=255)

                                                                              ctcCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAg              >  1:856490/1‑91 (MQ=255)

                                                                              ctcCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAg              >  1:6689/1‑91 (MQ=255)

                                                                               tcCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAgg             >  2:3047720/1‑91 (MQ=255)

                                                                                 cTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGaa           >  1:2869285/1‑91 (MQ=255)

                                                                                   tttttaATAAAGGAGATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCATAAAAATAAAAAAGGAata         >  2:683527/4‑91 (MQ=255)

                                                                                    tattaAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAatat        >  1:2412854/1‑91 (MQ=255)

                                                                                      ttaAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAATATAg      <  2:1004232/91‑1 (MQ=255)

                                                                                       taAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAATATAGc     <  2:3473625/91‑1 (MQ=255)

                                                                                       taAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAATATAGc     <  1:1250437/91‑1 (MQ=255)

                                                                                         aaaaaGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAATATAGCaa   >  1:2282509/1‑91 (MQ=255)

                                                                                         aaaaaGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAATATAGCaa   <  1:2498987/91‑1 (MQ=255)

CTGCAACTTCCATCGAACCGACAGTCCCTATTTTTTCGCCAATGACTTCCGATGTTGAACATCCAATTAAACATATATCTCCTTTATTAAAAAAGGACATATCTTTTAATTCGTCTAATAACATTGTCAAATCTTTCATAAAAGCCCACCCTTCCTAAAAATAAAAAAGGAATATAGCAAA  >  USA300TCH1516_ALE/2231437‑2231617

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