BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,905,775	G→C	5.1%	E128Q (GAA→CAA)	sigS →	RNA polymerase sigma factor SigS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,905,775	0	G	C	5.1%	895.1 / 25.0	249	E128Q (GAA→CAA)	sigS	RNA polymerase sigma factor SigS
Reads supporting (aligned to +/- strand): ref base G (117/119); new base C (8/5); total (125/124)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

CATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTAGT  >  USA300TCH1516_ALE/1905685‑1905865

cATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTg                                                                                            >  1:2859615/1‑91 (MQ=255)

cATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTg                                                                                            >  1:407436/1‑91 (MQ=255)

    aCCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAAt                                                                                        >  1:3530631/1‑91 (MQ=255)

     ccTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAAtt                                                                                       <  1:3410785/91‑1 (MQ=255)

       ttCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGc                                                                                     <  1:1168784/91‑1 (MQ=255)

       ttCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGc                                                                                     >  2:3639732/1‑91 (MQ=255)

       ttCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGc                                                                                     <  2:3853620/91‑1 (MQ=255)

       ttCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTCAAATTGc                                                                                     >  1:454033/1‑91 (MQ=255)

        tCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCt                                                                                    >  2:434471/1‑91 (MQ=255)

         cGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTg                                                                                   <  2:918307/91‑1 (MQ=255)

         cGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTg                                                                                   <  1:510606/91‑1 (MQ=255)

         cGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTg                                                                                   <  1:190406/91‑1 (MQ=255)

          gTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGa                                                                                  <  2:2457355/91‑1 (MQ=255)

            ttACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGaaa                                                                                >  2:2328552/1‑91 (MQ=255)

             tACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAAt                                                                               >  1:2100391/1‑91 (MQ=255)

              aCAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATc                                                                              >  1:80333/1‑91 (MQ=255)

              aCAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATc                                                                              <  2:533116/91‑1 (MQ=255)

               cAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCACACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCa                                                                             >  1:3082295/1‑91 (MQ=255)

               cAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCa                                                                             >  2:2976565/1‑91 (MQ=255)

               cAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCa                                                                             >  1:2187928/1‑91 (MQ=255)

                aaGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCAt                                                                            >  2:1785715/1‑91 (MQ=255)

                 aGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATg                                                                           >  1:3433041/1‑91 (MQ=255)

                      tCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCAtt                                                                      <  1:1108676/91‑1 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    >  1:3210963/1‑91 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    <  2:2159450/91‑1 (MQ=255)

                           aaGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGc                                                                 <  2:1083044/91‑1 (MQ=255)

                           aaGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGc                                                                 <  2:3575056/91‑1 (MQ=255)

                            aGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGTAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCt                                                                >  2:1077184/1‑91 (MQ=255)

                            aGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCt                                                                <  1:1226106/91‑1 (MQ=255)

                             gCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCtt                                                               >  2:2543286/1‑91 (MQ=255)

                              cTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCttt                                                              <  1:1684900/91‑1 (MQ=255)

                              cTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCttt                                                              <  2:3587112/91‑1 (MQ=255)

                              cTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCttt                                                              <  2:802119/91‑1 (MQ=255)

                                tttAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACATGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCa                                                            <  2:1367507/91‑1 (MQ=255)

                                 ttAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCaa                                                           >  2:3067293/1‑91 (MQ=255)

                                  tAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAAc                                                          >  1:1823523/1‑91 (MQ=255)

                                   aaaTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACg                                                         <  1:743280/91‑1 (MQ=255)

                                    aaTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGa                                                        <  2:1756515/91‑1 (MQ=255)

                                     aTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGat                                                       <  2:2681800/91‑1 (MQ=255)

                                     aTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGAGACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGat                                                       >  1:3357676/1‑91 (MQ=255)

                                      tCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGatt                                                      >  2:1843830/1‑91 (MQ=255)

                                       cAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGatta                                                     <  1:135048/91‑1 (MQ=255)

                                         aaGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaa                                                   <  2:2301894/91‑1 (MQ=255)

                                         aaGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaa                                                   >  2:71206/1‑91 (MQ=255)

                                         aaGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaa                                                   <  2:2637977/91‑1 (MQ=255)

                                          agagAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaat                                                  >  1:2677923/1‑91 (MQ=255)

                                           gagaAAGACTATGGCTCGAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaatt                                                 <  1:3328673/91‑1 (MQ=255)

                                           gagaAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaatt                                                 <  2:1980058/91‑1 (MQ=255)

                                            agaAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaatta                                                >  1:2325732/1‑91 (MQ=255)

                                             gaAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaattaa                                               >  2:3746035/1‑91 (MQ=255)

                                              aaaGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCCACGATTAaattaat                                              <  2:3300694/91‑1 (MQ=255)

                                               aaGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAaattaatt                                             >  1:917676/1‑91 (MQ=255)

                                               aaGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAaattaatt                                             <  2:1586236/91‑1 (MQ=255)

                                               aaGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAaattaatt                                             >  2:3341595/1‑91 (MQ=255)

                                               aaGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAaattaatt                                             <  1:3585620/91‑1 (MQ=255)

                                               aaGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAaattaatt                                             >  2:3660513/1‑91 (MQ=255)

                                               aaGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAaattaatt                                             >  2:3559067/1‑91 (MQ=255)

                                                aGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAattaatta                                            <  1:2034536/91‑1 (MQ=255)

                                                 gACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAAttaattaa                                           >  1:377699/1‑91 (MQ=255)

                                                 gACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAAttaattaa                                           >  1:3145206/1‑91 (MQ=255)

                                                 gACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAAttaattaa                                           <  1:1153351/91‑1 (MQ=255)

                                                  aCTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAg                                          <  2:1551216/91‑1 (MQ=255)

                                                    tATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTg                                                                                   <  1:2844804/48‑1 (MQ=255)

                                                    tATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGAt                                        <  2:2082697/91‑1 (MQ=255)

                                                         cTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCg                                   >  1:568235/1‑91 (MQ=255)

                                                            aaaCTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTa                                <  1:3159172/91‑1 (MQ=255)

                                                            aaaCTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTa                                <  1:1866693/91‑1 (MQ=255)

                                                              aCTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAg                              <  2:1325318/91‑1 (MQ=255)

                                                                 ataCCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGt                           >  1:2400139/1‑91 (MQ=255)

                                                                 ataCCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGt                           <  1:1035238/91‑1 (MQ=255)

                                                                  taCCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTa                          >  1:2869578/1‑91 (MQ=255)

                                                                  taCCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTa                          >  1:1522000/1‑91 (MQ=255)

                                                                   aCCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTaa                         >  1:853302/1‑91 (MQ=255)

                                                                      ttGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATg                      >  1:3007456/1‑91 (MQ=255)

                                                                        gAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGcc                    <  1:2146063/91‑1 (MQ=255)

                                                                        gAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGcc                    <  1:1052866/91‑1 (MQ=255)

                                                                        gAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGcc                    <  1:989537/91‑1 (MQ=255)

                                                                          aGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCaa                  <  2:2520353/91‑1 (MQ=255)

                                                                           ggATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAAc                 >  1:348267/1‑91 (MQ=255)

                                                                            gATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACa                <  1:935176/91‑1 (MQ=255)

                                                                             atatAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACAt               <  1:3492061/91‑1 (MQ=255)

                                                                             aTACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACAt               <  2:2835198/91‑1 (MQ=255)

                                                                             aTACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACAt               >  1:2016070/1‑91 (MQ=255)

                                                                             aTACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACAt               <  1:1182902/91‑1 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:1402624/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:2076534/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:2782047/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  2:3210404/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:493175/1‑91 (MQ=255)

                                                                                cAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAAt            >  1:2986425/1‑91 (MQ=255)

                                                                                cAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAAt            <  2:444741/91‑1 (MQ=255)

                                                                                   gCAATTTCAAATTGCTGAAATCATGTCATTATCCCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAAtttt         >  2:1272784/1‑91 (MQ=255)

                                                                                     aaTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTaa       >  2:1933182/1‑91 (MQ=255)

                                                                                      aTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTATGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      <  1:3265774/91‑1 (MQ=255)

                                                                                      aTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      >  2:2845342/1‑91 (MQ=255)

                                                                                      aTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      <  1:3407561/91‑1 (MQ=255)

                                                                                      aTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      <  1:3842212/91‑1 (MQ=255)

                                                                                      aTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      <  2:2011118/91‑1 (MQ=255)

                                                                                          gAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTAGt  >  1:2237440/1‑91 (MQ=255)

                                                                                          gAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTAGt  <  2:1241169/91‑1 (MQ=255)

CATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTAGT  >  USA300TCH1516_ALE/1905685‑1905865

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