BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,725,379	T→A	5.6%	intergenic (‑100/+240)	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,379	0	T	A	5.6%	602.2 / 18.0	177	intergenic (‑100/+240)	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 = 5.44e-01

Only 100 of 199 total aligned reads displayed.

CTATTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATC  >  USA300TCH1516_ALE/1725289‑1725468

cTATTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttat                                                                                            >  2:1360888/1‑91 (MQ=255)

cTATTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttat                                                                                            >  2:1595082/1‑91 (MQ=255)

 tATTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttatt                                                                                           <  2:495181/91‑1 (MQ=255)

  aTTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTttatta                                                                                          >  1:2463620/1‑91 (MQ=255)

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

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

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

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

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

       cTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaa                                                                                    >  1:2182490/1‑91 (MQ=255)

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

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

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

           ttttAAATTACCATTCCATAAAAATTGATGACACAGATAGTTGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaat                                                                                >  1:520678/1‑91 (MQ=255)

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

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

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

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

             ttAAATTACCATTCCATAAAAATTGATGACAGAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTAATAAt‑aaaaaatga                                                                              <  2:3551275/91‑3 (MQ=255)

              tAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTAATAAt‑aaaaaatgat                                                                             <  2:676514/91‑4 (MQ=255)

               aaaTTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaataaaa                                                                            <  2:235461/91‑1 (MQ=255)

                aaTTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaataaaaa                                                                           >  1:2604479/1‑91 (MQ=255)

                aaTTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATt‑aaaaaataaaaa                                                                           <  2:2476253/91‑1 (MQ=255)

                   tACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACg                                                                        >  1:3375594/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                                                                   >  1:2460588/1‑91 (MQ=255)

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

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

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

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

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

                           cATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTACAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCt                                                                <  1:3873220/91‑1 (MQ=255)

                               aaaaTTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCaaa                                                            <  1:2304140/91‑1 (MQ=255)

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

                                aaaTTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCaaaa                                                           <  2:3878560/91‑1 (MQ=255)

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

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

                                aaaTTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCaaaa                                                           <  1:2500149/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                          cacaGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACg                                                 <  1:848793/91‑1 (MQ=255)

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

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

                                                   tCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCgg                                        <  1:1991489/91‑1 (MQ=255)

                                                   tCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCgg                                        >  2:2910859/1‑91 (MQ=255)

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

                                                     gATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTa                                      <  2:3722895/91‑1 (MQ=255)

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

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

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

                                                       tttGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTacc                                    <  1:757152/91‑1 (MQ=255)

                                                       tttGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTacc                                    <  1:2104119/91‑1 (MQ=255)

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

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

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

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

                                                            aaaaCTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCt                               >  1:1791242/1‑91 (MQ=255)

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

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

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

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

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

                                                                 tAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTta                          <  2:3502909/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                       aaaTCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTTAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATg                    >  1:3867419/1‑91 (MQ=255)

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

                                                                          tCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCaa                 >  1:2230511/1‑91 (MQ=255)

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

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

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

                                                                               atCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCaa            <  1:889098/91‑1 (MQ=255)

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

                                                                                tCATTTTTTATTATT‑AAAAAAT‑AAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAaca          <  1:1925151/91‑1 (MQ=255)

                                                                                 cATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAGCGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAaca          <  2:2946804/91‑1 (MQ=255)

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

                                                                                   ttttttAATAAT‑AAAAAATAAAAAACGCCCATCCTCATAAGGGACGAACGTCATCGCGGTACCACACTAGTTATAAATGCAATTCTAcaca        >  1:3613887/1‑91 (MQ=255)

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

                                                                                    tttttATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACAt       >  1:470124/1‑91 (MQ=255)

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

                                                                                     ttttATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACAtt      <  1:1670924/91‑1 (MQ=255)

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

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

                                                                                        tattatt‑aAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTAt   >  2:3584275/1‑91 (MQ=255)

CTATTTACTGTTTTTAAATTACCATTCCATAAAAATTGATGACACAGATAGTCGATTTGCAAAACTAGTATAAATCAATATCATTTTTTATTATT‑AAAAAATAAAAAACGCCCATCCTCAAAAGGGACGAACGTTATCGCGGTACCACCCTAGTTATAAATGCAATTCAACACATTTATC  >  USA300TCH1516_ALE/1725289‑1725468

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