BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,479,038	C→T	8.7%	I486I (ATC→ATT)	USA300HOU_RS12765 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,479,038	0	C	T	8.7%	651.9 / 35.3	207	I486I (ATC→ATT)	USA300HOU_RS12765	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (91/98); new base T (10/8); total (101/106)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 217 total aligned reads displayed.

TGATAGATATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATC  >  USA300TCH1516_ALE/2478948‑2479128

tGATAGATATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATc                                                                                            >  1:3470250/1‑91 (MQ=255)

tGATAGATATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATc                                                                                            >  2:1442148/1‑91 (MQ=255)

   tagataTGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCaaa                                                                                         <  2:1866683/91‑1 (MQ=255)

   tagataTGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCaaa                                                                                         >  1:1399585/1‑91 (MQ=255)

    agataTGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAaat                                                                                        >  2:1424565/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

               tAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                             >  1:2553208/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                                                                                 >  1:842488/1‑85 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                         ttCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAGATCAGCtttt                                                                   <  2:3009337/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                    tGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCt                                                        >  2:1792982/1‑91 (MQ=37)

                                     gAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTa                                                       >  2:3065751/1‑91 (MQ=35)

                                      aaTGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTaa                                                      <  2:2582710/91‑1 (MQ=34)

                                       aTGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAg                                                     >  1:1301129/1‑91 (MQ=33)

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

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

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

                                                         tataAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCt                                   >  2:184487/1‑91 (MQ=33)

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

                                                             aaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttt                               >  1:1419992/1‑91 (MQ=39)

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

                                                             aaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttt                               >  1:2526776/1‑91 (MQ=39)

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

                                                                 aCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTg                           >  2:621491/1‑91 (MQ=255)

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

                                                                     atCAAAAAAGCTGATTTCTATCAATTAATTAATTGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                       <  1:3794889/91‑1 (MQ=255)

                                                                     atCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                       >  1:845779/1‑91 (MQ=255)

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

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

                                                                          aaaaGCTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaa                  >  1:3353271/1‑91 (MQ=18)

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

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

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

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

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

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

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

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

                                                                                 gATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt           >  2:3367504/1‑91 (MQ=21)

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

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

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

                                                                                  aTTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAtt          >  1:3707024/1‑91 (MQ=21)

                                                                                   tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTCTGTAAACCCAAATCGTATTc         >  2:1696491/1‑91 (MQ=255)

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

                                                                                    ttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCt        <  1:1227726/91‑1 (MQ=255)

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

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

                                                                                       tATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGtt     <  1:2824296/91‑1 (MQ=255)

                                                                                       tATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGtt     <  1:1093483/91‑1 (MQ=255)

                                                                                        aTCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTa    <  1:1873969/91‑1 (MQ=255)

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

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

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

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

TGATAGATATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATC  >  USA300TCH1516_ALE/2478948‑2479128

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