BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,479,044	A→T	9.7%	*488Y (TAA→TAT)	USA300HOU_RS12765 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,479,044	0	A	T	9.7%	600.9 / 42.9	186	*488Y (TAA→TAT)	USA300HOU_RS12765	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (81/87); new base T (10/8); total (91/95)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 213 total aligned reads displayed.

TATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACTT  >  USA300TCH1516_ALE/2478955‑2479134

tatGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAaagaat                                                                                           >  1:3133060/1‑91 (MQ=255)

 atGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatt                                                                                          >  2:3759653/1‑91 (MQ=255)

  tGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                         >  1:3242408/1‑91 (MQ=255)

  tGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                         >  2:468294/1‑91 (MQ=255)

  tGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                         >  2:2779512/1‑91 (MQ=255)

   ggTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                        <  1:1803461/91‑1 (MQ=255)

   ggTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                        <  1:17918/91‑1 (MQ=255)

    gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                       >  1:3057986/1‑91 (MQ=255)

    gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                       <  2:1418935/91‑1 (MQ=255)

    gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                       <  1:1068952/91‑1 (MQ=255)

    gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                       >  1:260695/1‑91 (MQ=255)

        tAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                   >  1:2935991/1‑91 (MQ=255)

        tAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                   <  1:444368/91‑1 (MQ=255)

         aaaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaaa                                                                                  <  2:3637300/91‑1 (MQ=255)

          aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                 <  2:3157785/91‑1 (MQ=18)

          aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                 <  2:321250/91‑1 (MQ=18)

          aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                       <  2:842488/85‑1 (MQ=255)

          aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                                 >  2:2711356/1‑91 (MQ=255)

          aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                                 >  1:1060834/1‑91 (MQ=255)

           aCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                                <  2:1060834/91‑1 (MQ=255)

            cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCa                                                                               <  2:1705964/91‑1 (MQ=18)

            cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                               >  2:424681/1‑91 (MQ=255)

            cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                               >  2:1609794/1‑91 (MQ=255)

             gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                              <  1:2526831/91‑1 (MQ=255)

             gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                              >  2:1774074/1‑91 (MQ=255)

             gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                              >  1:1857662/1‑91 (MQ=255)

              tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                             >  1:746971/1‑91 (MQ=255)

              tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                             >  2:767315/1‑91 (MQ=255)

              tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                             >  1:1483230/1‑91 (MQ=255)

                aCTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                           <  1:3862655/91‑1 (MQ=255)

                 cTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCAGCttt                                                                          <  1:2879575/91‑1 (MQ=17)

                 cTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCAGCttt                                                                          >  2:1024368/1‑91 (MQ=17)

                 cTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                          <  1:3760924/91‑1 (MQ=255)

                   tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                        <  1:897620/91‑1 (MQ=255)

                   tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                        <  1:136110/91‑1 (MQ=255)

                   tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                        <  1:2711356/91‑1 (MQ=255)

                    cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                       >  2:3667571/1‑91 (MQ=255)

                    cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                       <  2:3720869/91‑1 (MQ=255)

                    cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttt                                                                         >  2:2258639/1‑89 (MQ=255)

                    cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttt                                                                         <  1:2258639/89‑1 (MQ=255)

                      aGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTAc                                                                     >  1:2007934/1‑91 (MQ=255)

                      aGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTAc                                                                     <  2:1320562/91‑1 (MQ=255)

                       gCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACa                                                                    >  2:2410939/1‑91 (MQ=255)

                          atatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTg                                                                 <  2:1527770/91‑1 (MQ=40)

                           tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                                >  2:1934909/1‑91 (MQ=39)

                           tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                                <  1:1807179/91‑1 (MQ=39)

                             tGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCt                                                              <  1:2578519/91‑1 (MQ=37)

                              gAATGCGACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTa                                                             >  2:1844795/1‑91 (MQ=18)

                              gAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTa                                                             <  1:367659/91‑1 (MQ=35)

                                                atgataAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTGTACATTGCCTAAGAACTTAATGTCCCAGc                                           <  2:292131/88‑1 (MQ=12)

                                                aTTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGc                                           <  2:1280667/91‑1 (MQ=32)

                                                aTTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGc                                           <  1:2063078/91‑1 (MQ=32)

                                                aTTATAAAATACATAACAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGc                                           <  2:1136334/91‑1 (MQ=17)

                                                 ttATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGcc                                          <  1:116819/91‑1 (MQ=33)

                                                  tataAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCt                                         <  1:1789356/91‑1 (MQ=34)

                                                    taAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGt                                       <  1:1914788/91‑1 (MQ=37)

                                                     aaaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtt                                      <  1:2369943/91‑1 (MQ=38)

                                                      aaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttt                                     <  1:1885215/91‑1 (MQ=39)

                                                       aaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtttt                                    >  2:1209786/1‑91 (MQ=40)

                                                       aaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtttt                                    <  2:3864051/91‑1 (MQ=39)

                                                        aTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttttt                                   <  2:3344124/91‑1 (MQ=255)

                                                           cATATCACAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGt                                <  1:1931728/91‑1 (MQ=255)

                                                            atatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTa                               >  2:3009654/1‑91 (MQ=255)

                                                             tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaa                              >  1:659296/1‑91 (MQ=255)

                                                             tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaa                              <  1:1272158/91‑1 (MQ=255)

                                                              atCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                             <  2:1458381/91‑1 (MQ=255)

                                                               tCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAc                            >  2:2727259/1‑91 (MQ=255)

                                                               tCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAc                            >  1:3877545/1‑91 (MQ=255)

                                                                 aaaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAccc                          >  2:1871591/1‑91 (MQ=255)

                                                                   aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaa                        <  1:251066/91‑1 (MQ=255)

                                                                   aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaa                        >  2:2497034/1‑91 (MQ=255)

                                                                     aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAAt                      >  1:2130376/1‑91 (MQ=255)

                                                                     aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAAt                      >  1:3794521/1‑91 (MQ=255)

                                                                     aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAAt                      <  2:1101366/91‑1 (MQ=255)

                                                                      aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                     <  2:3057986/91‑1 (MQ=255)

                                                                      aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                     >  1:2117739/1‑91 (MQ=255)

                                                                      aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                     >  1:2841998/1‑91 (MQ=255)

                                                                       gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                    >  2:660603/1‑91 (MQ=255)

                                                                       gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                    <  1:1329139/91‑1 (MQ=255)

                                                                       gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                    <  1:1013434/91‑1 (MQ=255)

                                                                       gCTGATTTCTATCAAATAATTAATAGAAATCAGATTTTTTACATTGCCTAAGAACTTAAGGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                    <  1:3444983/91‑1 (MQ=255)

                                                                          gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt                 <  2:785046/91‑1 (MQ=255)

                                                                          gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt                 >  2:351461/1‑91 (MQ=255)

                                                                           aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAtt                <  2:1697955/91‑1 (MQ=255)

                                                                            tttCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc               >  1:845754/1‑91 (MQ=21)

                                                                            tttCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc               >  1:1798408/1‑91 (MQ=21)

                                                                            tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc               >  2:2559874/1‑91 (MQ=255)

                                                                            tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTCAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc               >  2:918040/1‑91 (MQ=255)

                                                                              tCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTg             <  2:675889/91‑1 (MQ=255)

                                                                               cTATCAAATAATTAATAGAAATCAGCTTTTTTGCATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGt            >  2:3816206/1‑91 (MQ=255)

                                                                                 aTCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTa          >  2:702726/1‑91 (MQ=255)

                                                                                  tCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTAt         <  2:1499454/91‑1 (MQ=255)

                                                                                   cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc               >  2:2680331/1‑84 (MQ=255)

                                                                                   cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc               <  1:2680331/84‑1 (MQ=255)

                                                                                   cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATc        <  2:116846/91‑1 (MQ=255)

                                                                                   cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATc        >  1:3346609/1‑91 (MQ=255)

                                                                                    aaaTAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCt       >  1:1520569/1‑91 (MQ=255)

                                                                                         attaatAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTCTTCTGTTATCTAACtt  >  2:2083061/1‑91 (MQ=255)

                                                                                         attaatAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACtt  >  1:1418289/1‑91 (MQ=255)

                                                                                         attaatAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACtt  >  2:3138442/1‑91 (MQ=255)

TATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACTT  >  USA300TCH1516_ALE/2478955‑2479134

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