BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	693,258	T→A	7.8%	I65F (ATT→TTT)	USA300HOU_RS03345 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	693,258	0	T	A	7.8%	812.5 / 44.3	245	I65F (ATT→TTT)	USA300HOU_RS03345	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (106/119); new base A (11/8); total (118/127)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 256 total aligned reads displayed.

TATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTT  >  USA300TCH1516_ALE/693168‑693348

tataAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAt                                                                                            >  1:3216300/1‑91 (MQ=255)

tataAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAt                                                                                            >  2:3034704/1‑91 (MQ=255)

  taAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAttt                                                                                          <  1:1086177/91‑1 (MQ=255)

  taAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAttt                                                                                          >  1:271859/1‑91 (MQ=255)

    aaTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTa                                                                                        >  1:1324984/1‑91 (MQ=255)

    aaTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTa                                                                                        >  1:664392/1‑91 (MQ=255)

    aaTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTa                                                                                        >  1:2329298/1‑91 (MQ=255)

      tACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAgg                                                                                      <  1:3673624/91‑1 (MQ=255)

       aCAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGc                                                                                     <  1:3554926/91‑1 (MQ=255)

        cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTacgct                                                                                    <  2:83550/91‑5 (MQ=38)

        cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTacgct                                                                                    <  2:1193985/91‑5 (MQ=38)

        cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTacgct                                                                                    <  1:2106360/91‑5 (MQ=38)

        cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTACGct                                                                                    <  1:3572758/91‑2 (MQ=38)

        cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGcc                                                                                    <  1:3534785/91‑1 (MQ=255)

        cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGcc                                                                                    >  1:961053/1‑91 (MQ=255)

        cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGcc                                                                                    <  2:70882/91‑1 (MQ=255)

            cTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGAcc                                                                                >  1:1931413/1‑91 (MQ=255)

             ttCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCt                                                                               <  1:3325679/91‑1 (MQ=255)

             ttCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCt                                                                               >  2:1830783/1‑91 (MQ=255)

              tCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTACGCTGACCtt                                                                              <  1:2031535/91‑1 (MQ=37)

               ccAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCttt                                                                             >  2:3772841/1‑91 (MQ=255)

                cAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCtttt                                                                            <  2:271859/91‑1 (MQ=255)

                   aTTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTAcc                                                                         <  2:1345409/91‑1 (MQ=255)

                    ttCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCg                                                                        <  1:768561/91‑1 (MQ=255)

                    ttCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCg                                                                        <  2:535973/91‑1 (MQ=255)

                    ttCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCg                                                                        <  1:2058036/91‑1 (MQ=255)

                      cACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTa                                                                      <  1:1517137/91‑1 (MQ=255)

                      cACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTa                                                                      <  2:2253114/91‑1 (MQ=255)

                       aCTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAt                                                                     <  2:3565580/91‑1 (MQ=255)

                       aCTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAt                                                                     <  2:761660/91‑1 (MQ=255)

                         ttGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAttt                                                                   <  2:1599821/91‑1 (MQ=255)

                          tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAtttt                                                                  >  1:722283/1‑91 (MQ=255)

                          tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAtttt                                                                  <  1:1297684/91‑1 (MQ=255)

                           gTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTc                                                                 >  2:2201061/1‑91 (MQ=255)

                           gTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTc                                                                 >  1:926467/1‑91 (MQ=255)

                               cATTCTGCCATGTCCTGCTTCACTCACGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTa                                                             >  2:3177268/1‑91 (MQ=255)

                                aTTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTaa                                                            <  2:722283/91‑1 (MQ=255)

                                aTTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTaa                                                            <  2:1284608/91‑1 (MQ=255)

                                  tCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACg                                                          >  1:1740912/1‑91 (MQ=255)

                                     gCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTaa                                                       >  2:3566675/1‑91 (MQ=255)

                                     gCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTaa                                                       >  1:3686940/1‑91 (MQ=255)

                                       cATGTCCTGCTGCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAGATTTTAGGCGGACCTTTTACCGTATTTTCTTTAACGTAACt                                                     <  2:2951392/91‑1 (MQ=255)

                                          gTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATc                                                  >  2:361975/1‑91 (MQ=255)

                                             cTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCAtt                                               <  1:1704011/91‑1 (MQ=255)

                                              tGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTa                                              <  2:2011534/91‑1 (MQ=255)

                                                 ttCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAAt                                           <  1:2501874/91‑1 (MQ=255)

                                                 ttCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAAt                                           >  1:3729709/1‑91 (MQ=255)

                                                    aCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTa                                        >  1:794077/1‑91 (MQ=255)

                                                    aCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTa                                        >  2:1694510/1‑91 (MQ=255)

                                                      tCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGa                                      >  2:1724130/1‑91 (MQ=255)

                                                       cAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGAc                                     <  2:1834847/91‑1 (MQ=255)

                                                       cAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGAc                                     >  2:2577405/1‑91 (MQ=255)

                                                        aaGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACt                                    >  2:2075169/1‑91 (MQ=255)

                                                        aaGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACt                                    >  2:2619806/1‑91 (MQ=255)

                                                         aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                   <  1:3773654/91‑1 (MQ=255)

                                                         aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                   <  1:1830783/91‑1 (MQ=255)

                                                         aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                   >  1:1015550/1‑91 (MQ=255)

                                                         aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                   <  1:347078/91‑1 (MQ=255)

                                                         aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                   <  1:457115/91‑1 (MQ=255)

                                                          gCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGt                                  <  2:2758815/91‑1 (MQ=255)

                                                           cTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTa                                 <  2:1162354/91‑1 (MQ=255)

                                                           cTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTa                                 <  2:2407425/91‑1 (MQ=255)

                                                            tCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTAt                                <  1:2863026/91‑1 (MQ=255)

                                                              aTACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTAttt                              >  2:1099553/1‑91 (MQ=255)

                                                                 cacaATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCtt                           >  1:1422644/1‑91 (MQ=255)

                                                                     aTAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGc                       <  2:2820786/91‑1 (MQ=255)

                                                                      tAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCt                      <  1:2214281/91‑1 (MQ=255)

                                                                       aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                     >  2:3655916/1‑91 (MQ=255)

                                                                       aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                     >  2:3200226/1‑91 (MQ=255)

                                                                       aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                     >  1:2773037/1‑91 (MQ=255)

                                                                       aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                     >  2:2928692/1‑91 (MQ=255)

                                                                       aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                     <  1:2562573/91‑1 (MQ=255)

                                                                        gACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaa                    >  1:2452237/1‑91 (MQ=255)

                                                                         aCAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaaa                   >  2:1303557/1‑91 (MQ=255)

                                                                         aCAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaaa                   >  1:405106/1‑91 (MQ=255)

                                                                          cAAAGGTCAGCGTAGATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGGCTGTGTTTCTTTAGCTagag                  <  2:58935/91‑1 (MQ=255)

                                                                          cAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAg                  <  1:2372920/91‑1 (MQ=255)

                                                                          cAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAg                  >  2:3831300/1‑91 (MQ=255)

                                                                          cAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCACTAAATTTAGACTGTATTTCTTTAGCTAAAg                  >  2:399727/1‑91 (MQ=255)

                                                                            aaGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTc                <  1:1375124/91‑1 (MQ=255)

                                                                             aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa               <  2:1611199/91‑1 (MQ=255)

                                                                             aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa               >  1:1196383/1‑91 (MQ=255)

                                                                             aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa               <  2:2051067/91‑1 (MQ=255)

                                                                              ggTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCAt              <  2:876636/91‑1 (MQ=255)

                                                                              ggTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCAt              <  2:2073619/91‑1 (MQ=255)

                                                                               gTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCAtt             >  1:2154386/1‑91 (MQ=255)

                                                                                 cAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTaa           >  2:3428305/1‑91 (MQ=255)

                                                                                 cAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTaa           <  1:1648135/91‑1 (MQ=255)

                                                                                  ggcctAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAAt          >  1:3162765/5‑91 (MQ=38)

                                                                                   gCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATg         <  2:3310655/91‑1 (MQ=255)

                                                                                   gCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATg         <  1:2928631/91‑1 (MQ=255)

                                                                                    cctAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGa        >  2:2365033/3‑91 (MQ=255)

                                                                                    cGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGa        <  1:1643339/91‑1 (MQ=255)

                                                                                       aaaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAAc     <  2:1931413/91‑1 (MQ=255)

                                                                                       aaaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAAc     <  1:241656/91‑1 (MQ=255)

                                                                                       aaaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAAc     <  1:1264151/91‑1 (MQ=255)

                                                                                        aaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACt    <  2:1421937/91‑1 (MQ=255)

                                                                                        aaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACt    <  1:760332/91‑1 (MQ=255)

                                                                                        aaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACt    >  1:3725825/1‑91 (MQ=255)

                                                                                          ttttAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACttt  >  1:198633/1‑91 (MQ=255)

TATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTT  >  USA300TCH1516_ALE/693168‑693348

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