BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	110,660	T→A	7.9%	I173I (ATT→ATA)	USA300HOU_RS00515 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	110,660	0	T	A	7.9%	1609.6 / 103.1	443	I173I (ATT→ATA)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (155/253); new base A (17/18); total (172/271)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 472 total aligned reads displayed.

TGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGA  >  USA300TCH1516_ALE/110570‑110749

  aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTc                                                                                         >  1:607688/1‑91 (MQ=255)

   aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCa                                                                                        <  1:4336222/91‑1 (MQ=255)

   aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCa                                                                                        <  2:5096108/91‑1 (MQ=255)

   aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCa                                                                                        <  2:2257790/91‑1 (MQ=255)

    aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                       >  2:4959991/1‑91 (MQ=255)

     aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                      <  1:507493/91‑1 (MQ=255)

      aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCatat                                                                                     <  2:2927522/91‑1 (MQ=255)

          aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                 >  1:3538291/1‑91 (MQ=255)

          aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                 >  2:1691637/1‑91 (MQ=255)

           aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATcc                                                                                >  1:4323706/1‑91 (MQ=255)

            aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCa                                                                               >  1:1009775/1‑91 (MQ=255)

            aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCa                                                                               <  1:1732556/91‑1 (MQ=255)

             aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                              >  1:5193532/1‑91 (MQ=255)

               aGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCagag                                                                            <  2:899227/91‑1 (MQ=255)

                  aaTCGAAAACATTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTg                                                                         <  1:5834012/91‑1 (MQ=255)

                   aTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                        >  2:1927526/1‑91 (MQ=255)

                   aTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                        <  2:1273476/91‑1 (MQ=255)

                    tCGAAATCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                       <  2:5361407/91‑1 (MQ=255)

                    tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                       <  1:707143/91‑1 (MQ=255)

                    tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                       >  1:4745654/1‑91 (MQ=255)

                    tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                       <  1:4686704/91‑1 (MQ=255)

                     cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                      <  1:565215/91‑1 (MQ=255)

                     cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                      <  2:551680/91‑1 (MQ=255)

                     cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                      <  1:108156/91‑1 (MQ=255)

                     cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                      >  1:327191/1‑91 (MQ=255)

                      gAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                     <  2:2340200/91‑1 (MQ=255)

                      gAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                     >  2:3615275/1‑91 (MQ=255)

                       aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                    <  1:5317801/91‑1 (MQ=255)

                       aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                    <  2:1726869/91‑1 (MQ=255)

                       aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                    <  2:3006156/91‑1 (MQ=255)

                        aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                   <  2:4857285/91‑1 (MQ=255)

                          aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTa                                                                 >  1:4613016/1‑91 (MQ=255)

                          aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTa                                                                 >  2:2148480/1‑91 (MQ=255)

                           cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                <  2:5132984/91‑1 (MQ=255)

                             ttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTc                                                              <  1:2956101/91‑1 (MQ=255)

                               aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                            <  2:259479/91‑1 (MQ=255)

                               aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                            >  1:877144/1‑91 (MQ=255)

                               aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                            <  2:748463/91‑1 (MQ=255)

                                 gTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                          <  2:2757187/91‑1 (MQ=255)

                                    cTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                       <  2:3841150/91‑1 (MQ=255)

                                     tttGTTCCATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                      <  1:262583/91‑1 (MQ=255)

                                     tttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                      >  1:3606548/1‑91 (MQ=255)

                                      ttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                     <  1:2344500/91‑1 (MQ=255)

                                       tGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                    <  1:5036406/91‑1 (MQ=255)

                                          tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                 <  1:5481554/91‑1 (MQ=255)

                                           cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                >  1:4275562/1‑91 (MQ=255)

                                           cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                <  2:4584617/91‑1 (MQ=255)

                                             atatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                              >  2:1594589/1‑91 (MQ=255)

                                               atGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                            <  1:5138323/91‑1 (MQ=255)

                                               atGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                            <  2:4278878/91‑1 (MQ=255)

                                                 ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                          <  2:3894596/91‑1 (MQ=255)

                                                 ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                          <  2:1512065/91‑1 (MQ=255)

                                                  gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                         >  2:5699219/1‑91 (MQ=255)

                                                    gACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                       >  1:2699061/1‑91 (MQ=255)

                                                    gACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                       <  2:3173177/91‑1 (MQ=255)

                                                    gACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                       >  1:31651/1‑91 (MQ=255)

                                                     aCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                      <  1:5570424/91‑1 (MQ=255)

                                                      cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                     >  2:406879/1‑91 (MQ=255)

                                                       tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                    <  1:1380670/91‑1 (MQ=255)

                                                       tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                    >  1:4420853/1‑91 (MQ=255)

                                                        ttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                   >  1:999166/1‑91 (MQ=255)

                                                         tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTGACTAATgatgat                                  >  2:3739204/1‑91 (MQ=255)

                                                         tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                  >  2:400477/1‑91 (MQ=255)

                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                 <  1:1202391/91‑1 (MQ=255)

                                                            aaaTTTGAAAGATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                               <  2:4666809/91‑1 (MQ=255)

                                                            aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                               >  2:1338266/1‑91 (MQ=255)

                                                            aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                               >  1:2146537/1‑91 (MQ=255)

                                                               tttGAAAAATTCTAAAGACGGAGATATTTCATATAATCCAGTGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                            <  2:5849512/91‑1 (MQ=255)

                                                                ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                           <  1:2041734/91‑1 (MQ=255)

                                                                 tGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGt                          <  1:4064175/91‑1 (MQ=255)

                                                                  gAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                         <  1:4585668/91‑1 (MQ=255)

                                                                   aaaaaTTATAGAGATGGAGATATTTCATATAATGCAGAGGTGCCGAGTTATTCGGCTAAGTTTCAATTAACTAACGATGATCATAATGTaa                        <  1:5365657/91‑1 (MQ=14)

                                                                    aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                       <  2:2667201/91‑1 (MQ=40)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                      >  1:4148525/1‑91 (MQ=40)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                      >  1:1524442/1‑91 (MQ=39)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATGCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                      <  1:5794012/91‑1 (MQ=255)

                                                                         tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                  <  2:2839248/91‑1 (MQ=39)

                                                                         tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                  <  2:5464137/91‑1 (MQ=39)

                                                                         tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                  <  1:977612/91‑1 (MQ=39)

                                                                            aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc               >  1:238362/1‑91 (MQ=39)

                                                                            aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc               >  1:3791833/1‑91 (MQ=39)

                                                                            aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc               >  2:516859/1‑91 (MQ=39)

                                                                            aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc               >  1:1664986/1‑91 (MQ=39)

                                                                            aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc               >  1:1908257/1‑91 (MQ=39)

                                                                              agatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc             <  1:199414/87‑1 (MQ=9)

                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc             >  2:3419769/1‑91 (MQ=39)

                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc             <  1:3078835/91‑1 (MQ=40)

                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc             <  2:1110091/91‑1 (MQ=40)

                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc             <  2:2532302/91‑1 (MQ=40)

                                                                                aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           <  2:2031722/91‑1 (MQ=255)

                                                                                 tggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa          <  2:2096361/90‑1 (MQ=255)

                                                                                 tggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa          >  2:3251502/2‑91 (MQ=255)

                                                                                 cGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACaa                   <  1:572391/82‑1 (MQ=40)

                                                                                 cGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa          <  2:4122039/91‑1 (MQ=255)

                                                                                  ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa         <  1:2518374/91‑1 (MQ=255)

                                                                                     gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat      <  2:2402363/91‑1 (MQ=255)

                                                                                     gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat      >  1:3359540/1‑91 (MQ=255)

                                                                                     gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat      >  2:1831329/1‑91 (MQ=255)

                                                                                      atatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgata     >  1:3930958/1‑91 (MQ=255)

                                                                                         tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatga  <  2:525926/91‑1 (MQ=255)

TGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGA  >  USA300TCH1516_ALE/110570‑110749

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