BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	747,521	A→T	5.5%	I92I (ATA→ATT)	USA300HOU_RS03630 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	747,521	0	A	T	5.5%	746.1 / 22.8	219	I92I (ATA→ATT)	USA300HOU_RS03630	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (107/99); new base T (8/4); total (115/104)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 236 total aligned reads displayed.

ATGGTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAGTTG  >  USA300TCH1516_ALE/747431‑747611

aTGGTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGAtata                                                                                             <  2:539561/91‑1 (MQ=255)

aTGGTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGAtata                                                                                             >  1:3205906/1‑91 (MQ=255)

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

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

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

    ttGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAata                                                                                         <  1:2689241/91‑1 (MQ=255)

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

     tGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAatat                                                                                        >  1:678903/1‑91 (MQ=255)

      gCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATc                                                                                       >  1:2532383/1‑91 (MQ=255)

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

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

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

          aTACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑aaaa                                                                                  >  2:3820751/1‑91 (MQ=255)

          aTACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑aaaa                                                                                  >  1:648278/1‑91 (MQ=255)

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

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

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

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

                ccTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATaa                                                                            >  1:2912331/1‑91 (MQ=255)

                   gAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACga                                                                         <  2:3264664/91‑2 (MQ=255)

                    aaaaGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCa                                                                        <  1:197618/91‑1 (MQ=255)

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

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

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

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

                         aaaCTTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAg                                                                   <  1:3225826/91‑1 (MQ=255)

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

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

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

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

                           aCGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGt                                                                 >  1:918426/1‑91 (MQ=255)

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

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

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

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

                                     ataACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTa                                                       <  2:623890/91‑1 (MQ=255)

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

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

                                          tGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTa                                                  >  2:3526184/1‑91 (MQ=255)

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

                                            tAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGa                                                <  2:963133/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                                          >  2:1801899/1‑91 (MQ=255)

                                                  ttGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATg                                          <  2:2345402/91‑1 (MQ=255)

                                                   tGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATGTC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGt                                         >  2:2994210/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                   ccATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCAt                         <  2:3092315/91‑1 (MQ=255)

                                                                   ccATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCAt                         >  1:1519918/1‑91 (MQ=255)

                                                                   ccATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCAt                         <  2:2430371/91‑1 (MQ=255)

                                                                   ccATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCAt                         >  1:1769721/1‑91 (MQ=255)

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

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

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

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

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

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

                                                                        aGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCa                    <  2:1121520/91‑1 (MQ=255)

                                                                        aGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCa                    <  2:1121684/91‑1 (MQ=255)

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

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

                                                                              ttattTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAgtggt              >  1:17564/1‑91 (MQ=255)

                                                                              ttattTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAgtggt              >  2:2463289/1‑91 (MQ=255)

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

                                                                               tattTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAgtggtg             <  2:1935255/91‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                  tttGATATAACTATC‑AAAAATATAACGCAATAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAat          <  1:539245/91‑1 (MQ=255)

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

                                                                                   ttGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCACGTACAGTAGATTTATGTTTAGATGCATTAGCATTACCAAGTGGTGAata         >  2:3607203/1‑91 (MQ=255)

                                                                                    tGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAatat        <  2:2140062/91‑1 (MQ=255)

                                                                                     gATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATg       <  1:2314543/91‑1 (MQ=255)

                                                                                      atatAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAata         >  1:989375/1‑88 (MQ=255)

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

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

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

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

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

                                                                                         taAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAGtt   >  1:1892962/1‑91 (MQ=255)

ATGGTTGCAAATACTACCTGAAAAGAAACGTTATAGTATAACTGTAATGTTTGATAATAAAGGCAATCCATTAGAATATTATTTTGATATAAATATC‑AAAAATATAACGCAAAAAGGTAATGCGCGTACAGTAGATTTATGTTTAGATGTTTTAGCATTACCAAGTGGTGAATATGAGTTG  >  USA300TCH1516_ALE/747431‑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: `‑`