BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,725,382	T→A	5.6%	intergenic (‑103/+237)	alaS ← / ← recD2	Alanine‑‑tRNA ligase/ATP‑dependent RecD‑like DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,725,382	0	T	A	5.6%	605.6 / 18.0	177	intergenic (‑103/+237)	alaS/recD2	Alanine‑‑tRNA ligase/ATP‑dependent RecD‑like DNA helicase
Reads supporting (aligned to +/- strand): ref base T (87/80); new base A (3/7); total (90/87)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 208 total aligned reads displayed.

TTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATCACTT  >  USA300TCH1516_ALE/1725292‑1725472

tttACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattat                                                                                             >  2:3718589/1‑91 (MQ=255)

tttACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattat                                                                                             <  1:469977/91‑1 (MQ=255)

 ttACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattatt                                                                                            <  1:177653/91‑1 (MQ=255)

 ttACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattatt                                                                                            >  1:1079234/1‑91 (MQ=255)

 ttACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattatt                                                                                            <  1:1203360/91‑1 (MQ=255)

  tACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattatt‑a                                                                                          <  1:2364655/91‑1 (MQ=255)

  tACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattatt‑a                                                                                          <  2:1890504/91‑1 (MQ=255)

  tACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttattatt‑a                                                                                          <  2:3011008/91‑1 (MQ=255)

   aCTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aa                                                                                         >  2:1459871/1‑91 (MQ=255)

   aCTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aa                                                                                         >  1:3576125/1‑91 (MQ=255)

      gTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaa                                                                                      >  1:663102/1‑91 (MQ=255)

      gTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTAATAAt‑aaaaa                                                                                      <  1:1851280/91‑1 (MQ=255)

       tttttAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaa                                                                                     <  2:275125/91‑1 (MQ=255)

        ttttAAACTACCATTCCATAAAAATTGATGACACAGATAGTAGATTTGAAAAACTAGTATAAATCAATATCATTTTTTAATAAt‑aaaaaat                                                                                    <  2:3437706/91‑1 (MQ=255)

         tttAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaata                                                                                   >  2:1638573/1‑91 (MQ=255)

         tttAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaata                                                                                   <  1:717123/91‑1 (MQ=255)

         tttAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTAATAAt‑aaaaaatg                                                                                   <  1:1495251/91‑2 (MQ=255)

          ttAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaataa                                                                                  >  2:3523104/1‑91 (MQ=255)

              aTTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaataaaaaa                                                                              <  1:475691/91‑1 (MQ=255)

              aTTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaataaaaaa                                                                              >  1:417881/1‑91 (MQ=255)

               ttACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAAc                                                                             <  2:356068/91‑1 (MQ=255)

                tACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACg                                                                            >  1:3375594/1‑91 (MQ=255)

                  ccATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCATTATCATTTTTTATTATT‑AAAAAATAAAAAACGcc                                                                          >  2:2371421/1‑91 (MQ=255)

                  ccATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGcc                                                                          >  1:1312579/1‑91 (MQ=255)

                  ccATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGcc                                                                          >  1:1218908/1‑91 (MQ=255)

                    aTTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCa                                                                        >  2:266789/1‑91 (MQ=255)

                     ttCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCAt                                                                       <  2:3492463/91‑1 (MQ=255)

                     ttCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCAt                                                                       >  1:231710/1‑91 (MQ=255)

                     ttCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCAt                                                                       >  2:3481870/1‑91 (MQ=255)

                      tCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATc                                                                      <  1:1155089/91‑1 (MQ=255)

                      tCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAACATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATc                                                                      <  1:3283600/91‑1 (MQ=255)

                       ccATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATcc                                                                     >  1:3058606/1‑91 (MQ=255)

                        cATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCt                                                                    <  2:1067209/91‑1 (MQ=255)

                        cATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCt                                                                    <  1:869290/91‑1 (MQ=255)

                        cATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCt                                                                    <  1:2537976/91‑1 (MQ=255)

                        cATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTAGAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCt                                                                    <  1:3873204/91‑1 (MQ=255)

                            aaaaTTGATGACACAGATAGTCGATTTGCAAAACTGGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCaaa                                                                >  1:2026657/1‑91 (MQ=255)

                            aaaaTTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCaaa                                                                >  2:2013466/1‑91 (MQ=255)

                             aaaTTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCaaaa                                                               >  2:2720195/1‑91 (MQ=255)

                             aaaTTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCaaaa                                                               >  2:53162/1‑91 (MQ=255)

                               aTTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAgg                                                             <  1:1673411/91‑1 (MQ=255)

                                ttGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAggg                                                            >  1:1484154/1‑91 (MQ=255)

                                 tgatgaCACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGa                                                           >  1:2791794/1‑91 (MQ=255)

                                  gatgaCACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGAc                                                          >  1:1553369/1‑91 (MQ=255)

                                  gatgaCACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGAc                                                          >  2:618872/1‑91 (MQ=255)

                                    tgaCACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGa                                                        <  2:740635/91‑1 (MQ=255)

                                    tgaCACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGa                                                        >  1:812293/1‑91 (MQ=255)

                                     gaCACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGaa                                                       <  1:2928426/91‑1 (MQ=255)

                                       cacaGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACg                                                     >  2:3418947/1‑91 (MQ=255)

                                        acaGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGt                                                    <  2:2395677/91‑1 (MQ=255)

                                         caGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGtt                                                   <  1:1668917/91‑1 (MQ=255)

                                          aGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTa                                                  <  1:2507404/91‑1 (MQ=255)

                                          aGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTa                                                  <  2:292429/91‑1 (MQ=255)

                                          aGATAGTCCATTTACAAAACTAGTAAAAATCAATATCAATTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTa                                                  <  2:2135440/91‑1 (MQ=255)

                                             tAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATcg                                               <  2:3738448/91‑1 (MQ=255)

                                              aGTCGATTTGCAAAACTAGTATAAGTCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATcgc                                              <  1:2796827/91‑1 (MQ=255)

                                                 cGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGt                                           >  1:1038699/1‑91 (MQ=255)

                                                 cGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGt                                           <  1:2479510/91‑1 (MQ=255)

                                                  gATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTa                                          >  1:939350/1‑91 (MQ=255)

                                                   aTTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTac                                         >  1:1073316/1‑91 (MQ=255)

                                                   aTTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTac                                         >  2:153185/1‑91 (MQ=255)

                                                     ttGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTacca                                       <  1:525516/91‑1 (MQ=255)

                                                      tGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAATAGGGACGAACGTTATCGCGGTaccac                                      >  2:2001452/1‑91 (MQ=255)

                                                      tGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTaccac                                      <  2:1429572/91‑1 (MQ=255)

                                                      tGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTaccac                                      <  1:3642978/91‑1 (MQ=255)

                                                         aaaaCTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCt                                   <  1:2177736/91‑1 (MQ=255)

                                                           aaCTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAg                                 >  2:3637380/1‑91 (MQ=255)

                                                            aCTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGt                                >  1:2116819/1‑91 (MQ=255)

                                                             cTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGtt                               >  2:2553127/1‑91 (MQ=255)

                                                              tAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTta                              >  2:1533402/1‑91 (MQ=255)

                                                                gTATAAATCAATATCATTTTTTATTATTAAAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTtat                             >  1:3403506/1‑91 (MQ=255)

                                                                gTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTtata                            >  1:2259532/1‑91 (MQ=255)

                                                                gTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTtata                            <  2:2725743/91‑1 (MQ=255)

                                                                 tataAATCAATATCATTTTATATTATT‑AAATAATAAAAAACGCCCATCCTCAACGGGGACGGACGTTATCGTGGTTCCACCCTAGTTATaa                           >  2:464799/1‑91 (MQ=255)

                                                                  ataAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATaaa                          >  2:1533073/1‑91 (MQ=255)

                                                                   taAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAAt                         <  1:1966793/91‑1 (MQ=255)

                                                                   taAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAAt                         <  1:2213554/91‑1 (MQ=255)

                                                                   taAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAAt                         <  1:3163543/91‑1 (MQ=255)

                                                                    aaaTCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATg                        >  1:376681/1‑91 (MQ=255)

                                                                      aTCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCa                      <  2:2266598/91‑1 (MQ=255)

                                                                         aaTATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAAtt                   >  2:3234010/1‑91 (MQ=255)

                                                                          atatCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTc                  <  2:553150/91‑1 (MQ=255)

                                                                           tatCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCa                 >  1:1489198/1‑91 (MQ=255)

                                                                            atCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGt                                <  1:3161494/75‑1 (MQ=255)

                                                                            atCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGt                                >  2:3161494/1‑75 (MQ=255)

                                                                            atCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCaa                >  1:2775070/1‑91 (MQ=255)

                                                                             tCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAac               >  1:312442/1‑91 (MQ=255)

                                                                              cATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAaca              >  1:1726080/1‑91 (MQ=255)

                                                                              cATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAaca              <  2:180751/91‑1 (MQ=255)

                                                                              cATTTTTTATTATT‑AAAAAATAAAAAACGCCCAGCCTCAAAAGGGACGAACGTGATCGCGGTACCACCCTAGTTATAAATGCAATTCAaca              <  1:3852524/91‑1 (MQ=255)

                                                                                 tttttATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACAt           <  1:1364954/91‑1 (MQ=255)

                                                                                  ttttAATAAT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACAtt          >  2:1911915/1‑91 (MQ=255)

                                                                                   tttATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACAttt         <  2:812293/91‑1 (MQ=255)

                                                                                    ttattatt‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTa        >  1:625643/1‑91 (MQ=255)

                                                                                        tatt‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTaaatg    <  2:3648137/91‑5 (MQ=255)

                                                                                         att‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATg                        >  2:2349088/1‑70 (MQ=255)

                                                                                         att‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATCACt   >  2:1287010/1‑91 (MQ=255)

                                                                                         att‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATCACt   <  1:1005583/91‑1 (MQ=255)

                                                                                          tt‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATCACtt  <  2:2804937/91‑1 (MQ=255)

                                                                                          tt‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATCACtt  <  1:3024144/91‑1 (MQ=255)

TTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATCACTT  >  USA300TCH1516_ALE/1725292‑1725472

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