BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	747,522	A→T	5.4%	N93Y (AAT→TAT)	USA300HOU_RS03630 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	747,522	0	A	T	5.4%	763.1 / 23.9	224	N93Y (AAT→TAT)	USA300HOU_RS03630	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (111/101); new base T (8/4); total (119/105)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 228 total aligned reads displayed.

TGGTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAGTTG  >  USA300TCH1516_ALE/747432‑747611

tGGTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATaa                                                                                            >  2:3621575/1‑91 (MQ=255)

 ggTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATaaa                                                                                           <  2:2005683/91‑1 (MQ=255)

   ttGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAata                                                                                         >  2:2244152/1‑91 (MQ=255)

       aaaTACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑aa                                                                                    <  1:1888887/91‑1 (MQ=255)

        aaTACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑aaa                                                                                   >  2:710804/1‑91 (MQ=255)

        aaTACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑aaa                                                                                   <  1:3165190/91‑1 (MQ=255)

        aaTACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑aaa                                                                                   >  2:2333679/1‑91 (MQ=255)

         aTACTACCTGAAAAGAATGGTTATAGTATAACTGTAATGTTTGATAATAAAAGCATTGCATTAGAATATTATTTTGATATAAATATC‑aaag                                                                                  <  2:1555108/91‑2 (MQ=255)

            ctacCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATTTATATC‑aaaataa                                                                               <  2:2681684/91‑4 (MQ=255)

             tacCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATTTATATC‑aaaataat                                                                              <  1:3666350/91‑5 (MQ=39)

             tacCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATTTATATC‑aaaataat                                                                              <  1:2640620/91‑5 (MQ=39)

              acCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAAtata                                                                             >  2:1166150/1‑91 (MQ=255)

              acCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAAtata                                                                             >  2:1281209/1‑91 (MQ=255)

                    aaaGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCaa                                                                       <  1:1025165/91‑1 (MQ=255)

                     aaGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCaaa                                                                      >  2:3122160/1‑91 (MQ=255)

                      aGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCaaaa                                                                     <  1:1662688/91‑1 (MQ=255)

                      aGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCaaaa                                                                     <  1:3417626/91‑1 (MQ=255)

                        aaaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAg                                                                   >  2:1122088/1‑91 (MQ=255)

                        aaaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAg                                                                   <  2:750755/91‑1 (MQ=255)

                        aaaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAg                                                                   <  2:1664812/91‑1 (MQ=255)

                         aaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAgg                                                                  >  1:1979239/1‑91 (MQ=255)

                         aaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAgg                                                                  >  2:211234/1‑91 (MQ=255)

                         aaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAgg                                                                  >  2:209170/1‑91 (MQ=255)

                         aaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAgg                                                                  >  1:3210931/1‑91 (MQ=255)

                         aaCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATA‑AAAAATATAACGCAAAAAgg                                                                  <  2:2452727/91‑1 (MQ=255)

                          aCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGt                                                                 <  2:1226809/91‑1 (MQ=255)

                           cGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTa                                                                <  1:961420/91‑1 (MQ=255)

                           cGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTa                                                                >  1:67624/1‑91 (MQ=255)

                           cGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTa                                                                <  2:594850/91‑1 (MQ=255)

                            gTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTaa                                                               >  1:2000993/1‑91 (MQ=255)

                            gTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTaa                                                               >  1:2140062/1‑91 (MQ=255)

                              tataGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATg                                                             >  2:2537364/1‑91 (MQ=255)

                              tataGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATg                                                             >  2:442216/1‑91 (MQ=255)

                               ataGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATgc                                                            >  1:2534159/1‑91 (MQ=255)

                                 aGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATgcgc                                                          >  1:1956799/1‑91 (MQ=255)

                                 aGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATgcgc                                                          >  2:1595145/1‑91 (MQ=255)

                                   tataACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGt                                                        >  2:3229943/1‑91 (MQ=255)

                                   tataACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGt                                                        <  2:1121618/91‑1 (MQ=255)

                                     taACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTAc                                                      <  1:996084/91‑1 (MQ=255)

                                     taACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGATAATGCGCGTAc                                                      <  1:2992461/91‑1 (MQ=255)

                                       aCTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAg                                                    <  1:3771930/91‑1 (MQ=255)

                                       aCTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAg                                                    <  1:3223266/91‑1 (MQ=255)

                                       aCTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAg                                                    >  1:2087842/1‑91 (MQ=255)

                                       aCTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAg                                                    <  1:825292/91‑1 (MQ=255)

                                         tGTAATGTGTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATGTAACGCAAAAAGGTAATGCGCGTACAGTa                                                  <  2:3371541/91‑1 (MQ=255)

                                           tAATGTTTGATAATAGAGGCAATCCTGTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGa                                                <  2:3359524/91‑1 (MQ=255)

                                           tAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGa                                                <  2:963133/91‑1 (MQ=255)

                                           tAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGa                                                <  2:1844542/91‑1 (MQ=255)

                                           tAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGa                                                <  1:2747808/91‑1 (MQ=255)

                                           tAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGa                                                <  1:2360488/91‑1 (MQ=255)

                                           tAATGCTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGa                                                <  2:3572540/91‑1 (MQ=255)

                                            aaTGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGAt                                               >  2:894077/1‑91 (MQ=255)

                                               gTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTa                                            >  2:3218211/1‑91 (MQ=255)

                                                 ttGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATg                                          >  1:2479257/1‑91 (MQ=255)

                                                 ttGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATg                                          >  2:1801899/1‑91 (MQ=255)

                                                  tGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGt                                         >  1:469814/1‑91 (MQ=255)

                                                   gATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGtt                                        <  1:3189022/91‑1 (MQ=255)

                                                   gATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGtt                                        <  2:3878510/91‑1 (MQ=255)

                                                   gATAATAAAGGCAATCCATTAGAATATTATTTAGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGtt                                        >  2:3816707/1‑91 (MQ=255)

                                                    ataataAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGttt                                       >  1:1073084/1‑91 (MQ=255)

                                                      aataaAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAg                                     <  1:480316/91‑1 (MQ=255)

                                                         aaaGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATg                                  >  1:3551599/1‑91 (MQ=255)

                                                         aaaGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATg                                  >  1:1677294/1‑91 (MQ=255)

                                                         aaaGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATg                                  >  2:1914227/1‑91 (MQ=255)

                                                         aaaGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATg                                  >  1:626717/1‑91 (MQ=255)

                                                           aGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGtt                                <  1:1897785/91‑1 (MQ=255)

                                                           aGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGtt                                >  1:937524/1‑91 (MQ=255)

                                                            ggCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGttt                               <  2:1028762/91‑1 (MQ=255)

                                                            ggCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGttt                               >  2:792583/1‑91 (MQ=255)

                                                            ggCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGttt                               >  1:688711/1‑91 (MQ=255)

                                                             gCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGtttt                              <  1:1426592/91‑1 (MQ=255)

                                                             gCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGtttt                              <  2:2984251/91‑1 (MQ=255)

                                                             gCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGtttt                              >  1:1874692/1‑91 (MQ=255)

                                                               aaTCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAg                            >  1:2933494/1‑91 (MQ=255)

                                                                 tCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCa                          >  2:1782158/1‑91 (MQ=255)

                                                                  ccAATAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCAt                         >  2:3413745/1‑91 (MQ=255)

                                                                   cATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCAtt                        <  1:3816707/91‑1 (MQ=255)

                                                                    aTTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTa                       >  2:457691/1‑91 (MQ=255)

                                                                    aTTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCACGTACAGTAGATTTATGTTTAGATGTTTTAGCATTa                       >  2:1700120/1‑91 (MQ=255)

                                                                     ttAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTAc                      >  1:2836196/1‑91 (MQ=255)

                                                                      tAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTAcc                     <  1:2446559/91‑1 (MQ=255)

                                                                      tAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTAcc                     <  1:1751515/91‑1 (MQ=255)

                                                                         aaTATTATTTTGATATAAATATCAAAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCaa                   >  2:3559403/1‑91 (MQ=255)

                                                                          atatTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTTGCATTACCAAgt                 <  1:2984162/91‑1 (MQ=255)

                                                                             ttattTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTAAGCATTACCAAgtggt              >  2:3705071/1‑91 (MQ=255)

                                                                              tattTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCGTTACCAAgtggtg             <  1:1756489/91‑1 (MQ=255)

                                                                               attTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGa            <  2:1624873/91‑1 (MQ=255)

                                                                                ttttGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGaa           <  1:1281209/91‑1 (MQ=255)

                                                                                 tttGATATTTATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAat          >  1:3627787/1‑91 (MQ=255)

                                                                                 tttGATATTTATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAat          >  2:1130840/1‑91 (MQ=255)

                                                                                 tttGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAat          <  2:3056276/91‑1 (MQ=255)

                                                                                  ttGATATTTATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAata         >  1:1881227/1‑91 (MQ=255)

                                                                                  ttGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAata         <  2:839503/91‑1 (MQ=255)

                                                                                  ttGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCACGTACAGTAGATTTATGTTTAGATGGCTTAGCATTACCAAGTGGTGAata         >  2:3609669/1‑91 (MQ=255)

                                                                                     atatGAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGa      <  2:2655486/91‑1 (MQ=255)

                                                                                     atatAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGa      >  2:3528656/1‑91 (MQ=255)

                                                                                     atatAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGa      <  1:402576/91‑1 (MQ=255)

                                                                                      tataAATATC‑AAAAATATAACGCAAAAGGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAg     >  1:295182/1‑91 (MQ=255)

                                                                                       ataAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAGt    >  1:2666467/1‑91 (MQ=255)

                                                                                         aaaTATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAGTTg  >  2:2421635/1‑91 (MQ=255)

TGGTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAGTTG  >  USA300TCH1516_ALE/747432‑747611

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