BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,375,929	T→C	6.3%	N32S (AAT→AGT)	topB ←	DNA topoisomerase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,375,929	0	T	C	6.3%	687.3 / 25.5	220	N32S (AAT→AGT)	topB	DNA topoisomerase 3
Reads supporting (aligned to +/- strand): ref base T (113/93); new base C (2/12); total (115/105)
Fisher's exact test for biased strand distribution p-value = 4.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 229 total aligned reads displayed.

AATCGCCATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATT  >  USA300TCH1516_ALE/2375839‑2376019

  tCGCCATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTAttt                                                                                          >  1:3453000/1‑91 (MQ=255)

  tCGCCATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTAttt                                                                                          >  2:1828564/1‑91 (MQ=255)

   cGCCATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTAtttt                                                                                         <  1:3390060/91‑1 (MQ=255)

   cGCCATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTAtttt                                                                                         >  1:1914144/1‑91 (MQ=255)

     ccATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCa                                                                                       <  1:639636/91‑1 (MQ=255)

      cATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCATGTGACAATATATTGGTTATTTTCaa                                                                                      >  1:3616199/1‑91 (MQ=255)

      cATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCaa                                                                                      <  2:1620844/91‑1 (MQ=255)

       aTTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCaaa                                                                                     >  2:2112349/1‑91 (MQ=255)

       aTTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCaaa                                                                                     <  1:1811705/91‑1 (MQ=255)

          ccTTTAATTTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                  <  1:1937275/91‑1 (MQ=255)

          ccTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                  >  1:3034842/1‑91 (MQ=255)

          ccTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                  >  2:1937242/1‑91 (MQ=255)

          ccTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                  >  2:3547406/1‑91 (MQ=255)

           cTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATaa                                                                                 <  1:3171673/91‑1 (MQ=255)

            tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAAc                                                                                >  2:3742610/1‑91 (MQ=255)

            tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAAc                                                                                <  1:1742584/91‑1 (MQ=255)

            tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAAc                                                                                <  2:1847877/91‑1 (MQ=255)

             ttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAAcc                                                                               >  2:2121361/1‑91 (MQ=255)

              tAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAACCa                                                                              <  2:1030466/91‑1 (MQ=255)

              tAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAACCa                                                                              <  1:1260377/91‑1 (MQ=255)

               aaaTTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCAt                                                                             >  1:1525479/1‑91 (MQ=255)

                aaTTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCAtt                                                                            <  1:3778193/91‑1 (MQ=255)

                  tttttATCGTATTGTTCAGGTGTCGCATTTGTCACTGGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTAc                                                                          <  2:3061770/91‑1 (MQ=255)

                  tttttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTAc                                                                          <  2:2099635/91‑1 (MQ=255)

                  tttttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTAc                                                                          <  1:112541/91‑1 (MQ=255)

                   ttttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACg                                                                         >  2:1635481/1‑91 (MQ=255)

                     ttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCt                                                                       >  2:233502/1‑91 (MQ=255)

                     ttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCt                                                                       >  2:247096/1‑91 (MQ=255)

                     ttATCGTATTGTTCAGGTGTCACATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCt                                                                       <  2:2777448/91‑1 (MQ=255)

                      tATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCtt                                                                      <  2:2924002/91‑1 (MQ=255)

                        tCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCt                                                                    <  1:2567679/91‑1 (MQ=255)

                        tCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCt                                                                    <  2:3746524/91‑1 (MQ=255)

                         cGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTg                                                                   <  1:705316/91‑1 (MQ=255)

                         cGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTg                                                                   >  1:3653295/1‑91 (MQ=255)

                          gTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGa                                                                  >  2:1562375/1‑91 (MQ=255)

                          gTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGa                                                                  >  1:425452/1‑91 (MQ=255)

                            aTTGTTCCGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGatt                                                                <  1:1942796/91‑1 (MQ=255)

                             ttGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGattt                                                               <  1:1705866/91‑1 (MQ=255)

                             ttGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGattt                                                               <  2:340461/91‑1 (MQ=255)

                              tGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGattta                                                              <  2:1195562/91‑1 (MQ=255)

                              tGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGattta                                                              <  2:1808740/91‑1 (MQ=255)

                               gTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGatttat                                                             <  2:2798412/91‑1 (MQ=255)

                                   aGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGt                                                         >  1:2681263/1‑91 (MQ=255)

                                   aGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGt                                                         >  2:2699832/1‑91 (MQ=255)

                                    ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                        >  1:2402919/1‑91 (MQ=255)

                                    ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                        <  2:21968/91‑1 (MQ=255)

                                    ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                        >  1:3292492/1‑91 (MQ=255)

                                    ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                        >  1:2076223/1‑91 (MQ=255)

                                    ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                        >  2:405177/1‑91 (MQ=255)

                                     gtgtCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTaa                                                       <  2:3699468/91‑1 (MQ=255)

                                       gtCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAg                                                     >  2:1404669/1‑91 (MQ=255)

                                        tCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGc                                                    <  1:2913853/91‑1 (MQ=255)

                                         cGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCa                                                   >  1:2008562/1‑91 (MQ=255)

                                         cGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCa                                                   >  1:2213222/1‑91 (MQ=255)

                                          gCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCAt                                                  <  1:85548/91‑1 (MQ=255)

                                           cATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATc                                                 >  2:3106704/1‑91 (MQ=255)

                                            aTTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCa                                                <  2:580041/91‑1 (MQ=255)

                                              ttGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGc                                              >  2:959213/1‑91 (MQ=255)

                                              ttGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGc                                              >  1:1687314/1‑91 (MQ=255)

                                               tGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCa                                             <  2:2844632/91‑1 (MQ=255)

                                                gTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCaa                                            <  1:1616913/91‑1 (MQ=255)

                                                gTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCaa                                            <  1:3251585/91‑1 (MQ=255)

                                                 tCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAat                                           <  1:1552067/91‑1 (MQ=255)

                                                 tCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAat                                           >  1:3200910/1‑91 (MQ=255)

                                                 tCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAat                                           <  2:713252/91‑1 (MQ=255)

                                                   aCTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAatat                                         >  2:2087208/1‑91 (MQ=255)

                                                   aCTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAatat                                         >  2:2412158/1‑91 (MQ=255)

                                                   aCTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAatat                                         >  2:2719195/1‑91 (MQ=255)

                                                    cTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGatttattt                                                           >  1:965254/1‑72 (MQ=255)

                                                    cTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATAtc                                        <  2:3292492/91‑1 (MQ=255)

                                                    cTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATAtc                                        >  1:296675/1‑91 (MQ=255)

                                                       gATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATctct                                     <  1:2341733/91‑1 (MQ=255)

                                                        aTGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCtt                                    <  2:3671903/91‑1 (MQ=255)

                                                         tGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTg                                   >  2:1931982/1‑91 (MQ=255)

                                                         tGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTg                                   <  2:3427547/91‑1 (MQ=255)

                                                             cTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATATGTAAAGCATCAGCAATATCTCTTGCGAc                               >  1:3108655/1‑91 (MQ=255)

                                                              tAACGCCCACGTGACAATATATTGATTATTTTCAAAGTAACCATTACGCTGCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACt                              <  2:3597749/91‑1 (MQ=255)

                                                                aCGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGa                            >  2:1835989/1‑91 (MQ=255)

                                                                aCGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGa                            >  2:3682195/1‑91 (MQ=255)

                                                                aCGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGa                            <  2:2008562/91‑1 (MQ=255)

                                                                 cGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGAt                           <  2:3738895/91‑1 (MQ=255)

                                                                 cGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGAt                           >  2:3524545/1‑91 (MQ=255)

                                                                  gCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATg                          >  1:2506761/1‑91 (MQ=255)

                                                                    ccACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGt                        <  1:3578412/91‑1 (MQ=255)

                                                                     cACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtt                       >  1:3213036/1‑91 (MQ=255)

                                                                       cGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtttt                     >  2:2627894/1‑91 (MQ=255)

                                                                       cGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtttt                     >  2:2587211/1‑91 (MQ=255)

                                                                         tGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTc                   >  1:1001967/1‑91 (MQ=255)

                                                                          gACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCa                  >  1:198585/1‑91 (MQ=255)

                                                                           aCAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAg                 >  2:229115/1‑91 (MQ=255)

                                                                            cAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGc                <  2:2143975/91‑1 (MQ=255)

                                                                             aaTATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCt               <  2:3112609/91‑1 (MQ=255)

                                                                              atatatTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATCGTTTTTCAGCTa              >  2:3223597/1‑91 (MQ=255)

                                                                               tatatTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTaa             >  2:2933031/1‑91 (MQ=255)

                                                                               tatatTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTaa             >  2:2931892/1‑91 (MQ=255)

                                                                                  atTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAatat          >  2:301060/1‑91 (MQ=255)

                                                                                   atGGTTACTTTGAAAATAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTCATAtt         >  2:1283075/2‑91 (MQ=255)

                                                                                    tGGTTACTTTGAAAATAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTa        >  1:2285982/1‑91 (MQ=255)

                                                                                       ttATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAg     >  1:1989996/1‑91 (MQ=255)

                                                                                       ttATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAg     <  1:1719814/91‑1 (MQ=255)

AATCGCCATTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATT  >  USA300TCH1516_ALE/2375839‑2376019

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