BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	837,540:1	+T	11.3%	intergenic (+34/‑235)	USA300HOU_RS04100 → / → uvrB	hypothetical protein/UvrABC system protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	837,540	1	.	T	11.3%	831.5 / 107.5	230	intergenic (+34/‑235)	USA300HOU_RS04100/uvrB	hypothetical protein/UvrABC system protein B
Reads supporting (aligned to +/- strand): ref base . (118/86); new base T (14/12); total (132/98)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

Only 100 of 251 total aligned reads displayed.

GGTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGCAC  >  USA300TCH1516_ALE/837450‑837630

 gttgttTAATTATAGATCAACGCGAATTTCATAAACAGCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATttg‑tt                                                                                          <  1:3678065/91‑1 (MQ=255)

 gTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttg‑t                                                                                           <  2:2786932/91‑1 (MQ=255)

  ttATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttg‑tt                                                                                          >  2:1146306/1‑91 (MQ=255)

    aTTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATttg‑tttga                                                                                       <  1:1534154/91‑4 (MQ=255)

     tttAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGaa                                                                                       <  1:1670934/91‑1 (MQ=255)

     tttAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGaa                                                                                       >  1:1886298/1‑91 (MQ=255)

      ttAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGaaa                                                                                      >  1:301757/1‑91 (MQ=255)

      ttAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGaaa                                                                                      >  1:2778025/1‑91 (MQ=255)

       tAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAAc                                                                                     >  1:1347349/1‑91 (MQ=255)

       tAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAAc                                                                                     <  2:298271/91‑1 (MQ=255)

       tAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAAc                                                                                     >  2:2897989/1‑91 (MQ=255)

       tAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAAc                                                                                     >  1:874732/1‑91 (MQ=255)

       tAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAAc                                                                                     >  1:2216655/1‑91 (MQ=255)

        aaTTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACa                                                                                    >  2:2683429/1‑91 (MQ=255)

         ttgttaGATCAGCTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACaa                                                                                   <  2:3092390/87‑1 (MQ=255)

         aTTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACat                                                                                   <  1:354695/91‑2 (MQ=38)

         aTTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACat                                                                                   <  2:2380947/91‑2 (MQ=38)

          ttATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACaaa                                                                                  >  2:1164198/1‑91 (MQ=255)

          ttATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACatt                                                                                  <  2:1347059/91‑3 (MQ=38)

          ttATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACatt                                                                                  <  2:1589212/91‑3 (MQ=38)

          ttATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACatt                                                                                  <  2:2797799/91‑3 (MQ=38)

          ttATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACatt                                                                                  <  1:1485637/91‑3 (MQ=38)

          ttATAGATCAACTCGAATTTCATAAAAATCAAGATGGAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACatt                                                                                  <  2:1856298/91‑3 (MQ=37)

           tataGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTC‑AACAAATTTGTTTGAAACatta                                                                                 <  1:230039/91‑4 (MQ=38)

             taGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAatg                                                                               <  1:3628947/91‑1 (MQ=255)

              aGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAatga                                                                              >  1:3561941/1‑91 (MQ=255)

              aGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAatga                                                                              >  1:1800388/1‑91 (MQ=255)

              aGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAatga                                                                              >  2:2498976/1‑91 (MQ=255)

                     cTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTaa                                                                       <  2:874732/91‑1 (MQ=255)

                      tCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAat                                                                      <  1:2134204/91‑1 (MQ=255)

                        gAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAatat                                                                    <  2:2368438/91‑1 (MQ=255)

                        gAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAatat                                                                    >  2:1059502/1‑91 (MQ=255)

                         aaTTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAAtata                                                                   >  2:2546699/1‑91 (MQ=255)

                         aaTTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAAtata                                                                   <  2:3488556/91‑1 (MQ=255)

                            ttCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAAtg                                                                >  1:3529286/1‑91 (MQ=255)

                             tCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAAtgt                                                               >  1:2711701/1‑91 (MQ=255)

                              cATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAAtgtg                                                              >  2:2466558/1‑91 (MQ=255)

                               aTAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATgtgt                                                             >  1:2557077/1‑91 (MQ=255)

                               aTAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATgtgt                                                             >  1:1758664/1‑91 (MQ=255)

                               aTAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATgtgt                                                             >  1:683461/1‑91 (MQ=255)

                                tAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTa                                                            >  1:810519/1‑91 (MQ=255)

                                   aCATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTAttt                                                         <  2:1351962/91‑1 (MQ=255)

                                      tCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACa                                                      >  2:2631920/1‑91 (MQ=255)

                                      tCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACa                                                      <  1:1920410/91‑1 (MQ=255)

                                      tCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACa                                                      <  2:2905789/91‑1 (MQ=255)

                                       cAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACAt                                                     <  1:982153/91‑1 (MQ=255)

                                        aaGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTCG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATa                                                    >  2:3408423/1‑91 (MQ=255)

                                        aaGATGAAGATCGATAAGATTAACTTTAAAATAATGTTTCAAAGAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATa                                                    <  1:3873940/91‑1 (MQ=255)

                                         aGATGAAGATCGCTAAGATTAACTTTTAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATAc                                                   <  1:33446/91‑1 (MQ=255)

                                          gATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACt                                                  <  1:3494905/91‑1 (MQ=255)

                                          gATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACt                                                  >  1:3386191/1‑91 (MQ=255)

                                           aTGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTa                                                 >  2:183178/1‑91 (MQ=255)

                                              aaGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaaaa                                              >  1:3790760/1‑91 (MQ=255)

                                                 aTCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATa                                           >  1:1036979/1‑91 (MQ=255)

                                                 aTCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATa                                           <  1:1998539/91‑1 (MQ=255)

                                                  tCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATaa                                          <  2:516331/91‑1 (MQ=255)

                                                  tCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATaa                                          >  1:3107947/1‑91 (MQ=255)

                                                   cGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATgtgt                                                             <  2:432830/71‑1 (MQ=255)

                                                   cGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAAc                                         >  2:1696523/1‑91 (MQ=255)

                                                        aGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGAt                                    >  1:1588994/1‑91 (MQ=255)

                                                         gATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATa                                   >  2:650765/1‑91 (MQ=255)

                                                         gATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATa                                   <  1:1849734/91‑1 (MQ=255)

                                                          aTTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATaa                                  <  1:3380369/91‑1 (MQ=255)

                                                           ttAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATaaa                                 <  2:3287973/91‑1 (MQ=255)

                                                            tAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAAAGTGTATTTACATACTAAAAATAACAAGATAAAc                                <  2:3430418/91‑1 (MQ=255)

                                                             aaCTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACg                               >  2:2334706/1‑91 (MQ=255)

                                                              aCTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACga                              >  1:2478908/1‑91 (MQ=255)

                                                              aCTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACga                              >  1:1901610/1‑91 (MQ=255)

                                                              aCTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACga                              >  1:1904132/1‑91 (MQ=255)

                                                              aCTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACga                              >  1:2976122/1‑91 (MQ=255)

                                                              aCTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACga                              <  2:3774347/91‑1 (MQ=255)

                                                              aCTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACga                              <  2:1214101/91‑1 (MQ=255)

                                                               cTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgat                             >  2:2702376/1‑91 (MQ=255)

                                                                 ttAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgattt                           >  1:1093622/1‑91 (MQ=255)

                                                                  tAAAATAATGTTTCAATCAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttg                          <  2:3714804/91‑1 (MQ=255)

                                                                  tAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttg                          >  2:1157440/1‑91 (MQ=255)

                                                                  tAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttg                          <  2:3117057/91‑1 (MQ=255)

                                                                    aaaTAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttgat                        <  2:587630/91‑1 (MQ=255)

                                                                       taatGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTa                     <  2:2245153/91‑1 (MQ=255)

                                                                         atGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAg                   >  1:3582018/1‑91 (MQ=255)

                                                                          tGTTTC‑AACAAATTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAgg                  >  2:2135361/1‑91 (MQ=39)

                                                                           gTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGc                 <  2:762292/91‑1 (MQ=255)

                                                                           gTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGc                 >  1:806928/1‑91 (MQ=255)

                                                                           gTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGc                 <  1:650765/91‑1 (MQ=255)

                                                                           gTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGc                 >  1:2072346/1‑91 (MQ=255)

                                                                            tttCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCa                <  2:3518524/91‑1 (MQ=255)

                                                                             ttCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCaa               >  1:3332381/1‑91 (MQ=255)

                                                                              tCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAg              <  1:1146306/91‑1 (MQ=255)

                                                                               tcaacaaaTTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAg              >  1:121995/3‑91 (MQ=255)

                                                                               caaacaaaTTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGc             >  2:3032935/1‑91 (MQ=255)

                                                                               caaacaaaTTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGc             >  1:47017/1‑91 (MQ=255)

                                                                                 aacaaaTTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCa            >  1:196366/1‑91 (MQ=255)

                                                                                    aaaTTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGt        >  1:2854446/1‑91 (MQ=255)

                                                                                    aaaTTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGt        >  1:3702454/1‑91 (MQ=255)

                                                                                     aaTTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGtt       >  2:728000/1‑91 (MQ=255)

                                                                                        ttg‑ttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGc    <  1:2631920/91‑1 (MQ=255)

                                                                                         tg‑ttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGca   >  1:436896/1‑91 (MQ=255)

                                                                                          g‑ttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGcac  >  1:3012443/1‑91 (MQ=255)

                                                                                          g‑ttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGcac  >  2:1968429/1‑91 (MQ=255)

                                                                                          g‑ttgAAACAAAATGATGATCAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGcac  >  2:3394018/1‑91 (MQ=255)

GGTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTG‑TTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGCAC  >  USA300TCH1516_ALE/837450‑837630

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