BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,513,975	C→T	6.5%	R54Q (CGA→CAA)	USA300HOU_RS12940 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,513,975	0	C	T	6.5%	845.2 / 30.3	261	R54Q (CGA→CAA)	USA300HOU_RS12940	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (131/113); new base T (14/3); total (145/116)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 267 total aligned reads displayed.

AAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCATAAC  >  USA300TCH1516_ALE/2513885‑2514065

aaaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATc                                                                                            <  2:1155168/91‑1 (MQ=255)

aaaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATc                                                                                            <  2:3240759/91‑1 (MQ=255)

aaaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATc                                                                                            <  2:3239169/91‑1 (MQ=255)

 aaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCg                                                                                           >  2:1001590/1‑91 (MQ=255)

 aaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCg                                                                                           <  2:1750730/91‑1 (MQ=255)

 aaCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCg                                                                                           >  2:1678803/1‑91 (MQ=255)

  aCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGa                                                                                          >  2:1892010/1‑91 (MQ=255)

  aCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGa                                                                                          >  2:733269/1‑91 (MQ=255)

   ccAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCGATTTAAATTGat                                                                                         <  2:1213637/91‑1 (MQ=255)

   ccAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGat                                                                                         >  1:1340755/1‑91 (MQ=255)

    cAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGatt                                                                                        >  2:669368/1‑91 (MQ=255)

      aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                                      >  2:112798/1‑91 (MQ=255)

      aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                                      >  1:3031016/1‑91 (MQ=255)

      aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                                      <  2:2933781/91‑1 (MQ=255)

      aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                                      >  2:2202453/1‑91 (MQ=255)

       cATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCGATTTAAATTGattatt                                                                                     <  1:2314056/91‑1 (MQ=255)

       cATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattatt                                                                                     <  1:2486721/91‑1 (MQ=255)

         ttGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCGATTTAAATTGATTAtttt                                                                                   <  1:3580849/91‑1 (MQ=255)

          tGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAttttt                                                                                  >  1:1589743/1‑91 (MQ=255)

          tGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAttttt                                                                                  <  1:2876911/91‑1 (MQ=255)

              tATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTttatta                                                                              >  1:1330792/1‑91 (MQ=255)

               aTTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTAc                                                                             <  1:3751003/91‑1 (MQ=255)

                     aTACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGaaa                                                                       <  1:405718/91‑1 (MQ=255)

                     aTACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGaaa                                                                       >  1:3176803/1‑91 (MQ=255)

                     aTACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGaaa                                                                       >  2:1796846/1‑91 (MQ=255)

                      tACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAAc                                                                      <  1:2325025/91‑1 (MQ=255)

                       aCGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACa                                                                     >  2:2008187/1‑91 (MQ=255)

                        cGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACAt                                                                    <  2:510037/91‑1 (MQ=255)

                         gCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATa                                                                   >  2:3516092/1‑91 (MQ=255)

                            ttCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTg                                                                <  2:1492233/91‑1 (MQ=255)

                            ttCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTg                                                                <  2:232655/91‑1 (MQ=255)

                             tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                               <  2:448522/91‑1 (MQ=255)

                             tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                               <  2:1454015/91‑1 (MQ=255)

                             tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                               <  2:2727763/91‑1 (MQ=255)

                              ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                              <  1:843324/91‑1 (MQ=255)

                              ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                              <  2:1948011/91‑1 (MQ=255)

                              ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                              >  2:3631523/1‑91 (MQ=255)

                              ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                              <  1:1017613/91‑1 (MQ=255)

                              ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                              >  1:1051601/1‑91 (MQ=255)

                               cTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCa                                                             >  1:2773486/1‑91 (MQ=255)

                                 gTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCCATTTAAATCGATTATTTTTATTACCTGAAACATACTGCCAtt                                                           >  1:3861855/1‑91 (MQ=255)

                                    ttCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTa                                                        >  2:2128160/1‑91 (MQ=255)

                                       aTAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTg                                                     <  2:1311742/91‑1 (MQ=255)

                                        tAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTgg                                                    <  1:1482098/91‑1 (MQ=255)

                                         aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                                   <  1:669368/91‑1 (MQ=255)

                                         aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                                   >  2:1409581/1‑91 (MQ=255)

                                         aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                                   <  2:1099762/91‑1 (MQ=255)

                                         aaCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGga                                                   >  1:1885716/1‑91 (MQ=255)

                                           cAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgaga                                                 <  2:2647804/91‑1 (MQ=255)

                                           cAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgaga                                                 <  1:2964443/91‑1 (MQ=255)

                                           cAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgaga                                                 >  2:3357769/1‑91 (MQ=255)

                                            aGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgat                                                <  2:135190/91‑1 (MQ=255)

                                              ttCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                              <  2:3390361/91‑1 (MQ=255)

                                              ttCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                              >  1:3027909/1‑91 (MQ=255)

                                               tCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatgat                                             >  2:545284/1‑91 (MQ=255)

                                               tCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatgat                                             >  1:2551551/1‑91 (MQ=255)

                                                ccAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGatgatg                                            >  1:2107730/1‑91 (MQ=255)

                                                ccAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGatgatg                                            >  1:3055310/1‑91 (MQ=255)

                                                ccAGCTAATCCACGACCCTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACTTGAAACATAGTGCCATTTTATTGGAGatgatg                                            >  1:337791/1‑91 (MQ=255)

                                                 cAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATgg                                           <  2:2270659/91‑1 (MQ=255)

                                                  aGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGt                                          <  1:3256887/91‑1 (MQ=255)

                                                   gCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGtt                                         <  1:560715/91‑1 (MQ=255)

                                                    cTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGttt                                        >  2:2371968/1‑91 (MQ=255)

                                                    cTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGttt                                        >  1:891786/1‑91 (MQ=255)

                                                     tAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTg                                       >  1:2200539/1‑91 (MQ=255)

                                                     tAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTg                                       <  2:2432629/91‑1 (MQ=255)

                                                       aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                                     <  2:3876917/91‑1 (MQ=255)

                                                       aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                                     >  2:1094147/1‑91 (MQ=255)

                                                         ccACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAg                                   >  1:3505775/1‑91 (MQ=255)

                                                          cACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGa                                  >  1:3342611/1‑91 (MQ=255)

                                                            cGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGAtt                                >  2:81864/1‑91 (MQ=255)

                                                            cGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTATGTTAGAtt                                <  1:1954126/91‑1 (MQ=255)

                                                             gACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTc                               <  2:956127/91‑1 (MQ=255)

                                                              aCCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCa                              >  1:1710325/1‑91 (MQ=255)

                                                                  tttCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGa                          <  1:1265255/91‑1 (MQ=255)

                                                                    tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                        <  2:2158798/91‑1 (MQ=255)

                                                                    tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                        <  2:2012140/91‑1 (MQ=255)

                                                                    tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                        <  2:1219607/91‑1 (MQ=255)

                                                                     cTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAg                       >  1:3105064/1‑91 (MQ=255)

                                                                      ttAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGc                      <  1:230729/91‑1 (MQ=255)

                                                                      ttAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCAAAGAAAGc                      >  2:3442980/1‑91 (MQ=255)

                                                                       tAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCg                     <  2:3312644/91‑1 (MQ=255)

                                                                        aaaaTAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGa                    <  2:2152177/91‑1 (MQ=255)

                                                                         aaaTAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGAt                   >  2:3489855/1‑91 (MQ=255)

                                                                           ataatCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTg                 >  1:1639018/1‑91 (MQ=255)

                                                                           ataatCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTg                 >  2:2615881/1‑91 (MQ=255)

                                                                             aatCAATTTAAATCGCTTATTTTTACAACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGcc               <  2:3357194/91‑1 (MQ=255)

                                                                              atCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGcca              >  2:418820/1‑90 (MQ=255)

                                                                              atCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCg              <  2:2587711/91‑1 (MQ=255)

                                                                               tCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGc             >  1:2464326/1‑91 (MQ=255)

                                                                               tCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGc             >  1:727590/1‑91 (MQ=255)

                                                                                cAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCa            >  1:308880/1‑91 (MQ=255)

                                                                                 aaTTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCaa           >  1:926731/1‑91 (MQ=255)

                                                                                  aTTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCaaa          >  1:1108795/1‑91 (MQ=255)

                                                                                  aTTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCaaa          <  2:789152/91‑1 (MQ=255)

                                                                                  aTTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCaaa          >  1:1137327/1‑91 (MQ=255)

                                                                                  aTTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCaaa          <  2:90880/91‑1 (MQ=255)

                                                                                     tAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATc       <  2:3144210/91‑1 (MQ=255)

                                                                                       aaTCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCAt     <  2:3676178/91‑1 (MQ=255)

                                                                                       aaTCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCAt     >  2:2864935/1‑91 (MQ=255)

AAACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCATAAC  >  USA300TCH1516_ALE/2513885‑2514065

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