BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,331,236	T→A	8.9%	S189T (TCT→ACT)	USA300HOU_RS06555 →	Monoacylglycerol lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,331,236	0	T	A	8.9%	711.3 / 46.1	225	S189T (TCT→ACT)	USA300HOU_RS06555	Monoacylglycerol lipase
Reads supporting (aligned to +/- strand): ref base T (117/88); new base A (9/11); total (126/99)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 233 total aligned reads displayed.

GCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATT  >  USA300TCH1516_ALE/1331146‑1331326

gCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGAttt                                                                                            <  2:642950/91‑1 (MQ=255)

     gaAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGt                                                                                       <  2:2161781/91‑1 (MQ=255)

     gaAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGt                                                                                       <  2:1124925/91‑1 (MQ=255)

      aaaCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAgaa                                                                                      <  2:3069208/91‑3 (MQ=255)

         cGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATc                                                                                   >  2:1224698/1‑91 (MQ=255)

         cGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATc                                                                                   <  1:3645087/91‑1 (MQ=255)

          gAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATcc                                                                                  >  2:3630402/1‑91 (MQ=255)

          gAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATcc                                                                                  >  1:146566/1‑91 (MQ=255)

          gAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATcc                                                                                  <  2:133004/91‑1 (MQ=255)

          gAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATcc                                                                                  <  1:2724682/91‑1 (MQ=255)

          gAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATcc                                                                                  >  1:2385634/1‑91 (MQ=255)

          gAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAGAAATcc                                                                                  <  1:1005262/91‑1 (MQ=255)

           aGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCa                                                                                 >  2:3026917/1‑91 (MQ=255)

             ttGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAAt                                                                               >  2:3318131/1‑91 (MQ=255)

             ttGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGAATTCTAGTAATCCAAt                                                                               >  2:83832/1‑91 (MQ=255)

              tGAATGGGTTAACCAGTTAGTATCAAATAGCTTTAATAAAAGAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAAtt                                                                              >  2:732809/1‑91 (MQ=255)

                aaTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGa                                                                            >  2:3547048/1‑91 (MQ=255)

                aaTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGa                                                                            >  2:2705323/1‑91 (MQ=255)

                 aTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGaa                                                                           >  1:3062822/1‑91 (MQ=255)

                 aTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGaa                                                                           >  1:2139405/1‑91 (MQ=255)

                 aTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGaa                                                                           >  1:3700700/1‑91 (MQ=255)

                    ggTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTa                                                                        <  1:3397413/91‑1 (MQ=255)

                    ggTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTa                                                                        >  1:3747386/1‑91 (MQ=255)

                    ggTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTa                                                                        <  2:975982/91‑1 (MQ=255)

                      ttAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGa                                                                      >  2:1849550/1‑91 (MQ=255)

                       tAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGAt                                                                     <  1:2679835/91‑1 (MQ=255)

                        aaCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATa                                                                    >  1:3326108/1‑91 (MQ=255)

                         aCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATaa                                                                   >  2:3193850/1‑91 (MQ=255)

                         aCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATaa                                                                   <  1:1259027/91‑1 (MQ=255)

                         aCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATaa                                                                   >  2:32183/1‑91 (MQ=255)

                         aCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATaa                                                                   >  2:276470/1‑91 (MQ=255)

                          ccAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAAc                                                                  <  2:445776/91‑1 (MQ=255)

                          ccAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAAc                                                                  <  1:231339/91‑1 (MQ=255)

                           cAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACt                                                                 >  2:1624598/1‑91 (MQ=255)

                            aGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACtt                                                                >  2:1254426/1‑91 (MQ=255)

                              ttAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTtta                                                              <  2:86794/91‑1 (MQ=255)

                              ttAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTtta                                                              >  1:1080050/1‑91 (MQ=255)

                              ttAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTtta                                                              <  2:249462/91‑1 (MQ=255)

                                aGTATCAAATAGATTTAATAAATAAATACGTCCATTACGTACACAAAGTTATTGGATTTCTAGTAATACAATTGAAGTAGATAACTttatt                                                            <  2:2947447/91‑1 (MQ=255)

                                    tCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAg                                                        <  1:2938605/91‑1 (MQ=255)

                                     cAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGa                                                       <  2:487770/91‑1 (MQ=255)

                                     cAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGa                                                       >  2:3704355/1‑91 (MQ=255)

                                      aaaTAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGAt                                                      >  2:1054136/1‑91 (MQ=255)

                                      aaaTAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGAt                                                      <  1:261221/91‑1 (MQ=255)

                                      aaaTAGTATTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGAt                                                      <  1:1942636/91‑1 (MQ=255)

                                       aaTAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATc                                                     <  2:1657643/91‑1 (MQ=255)

                                       aaTAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATc                                                     >  1:209443/1‑91 (MQ=255)

                                         tAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCa                                                   >  1:1739580/1‑91 (MQ=255)

                                            ttttAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCatat                                                >  2:1919949/1‑91 (MQ=255)

                                             tttAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCAtata                                               <  2:1533113/91‑1 (MQ=255)

                                             tttAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCAtata                                               >  1:15417/1‑91 (MQ=255)

                                              ttAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAg                                              <  2:3340350/91‑1 (MQ=255)

                                               taataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGt                                             <  1:1224698/91‑1 (MQ=255)

                                                   aaaaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGAt                                         <  2:2228717/91‑1 (MQ=255)

                                                   aaaaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGAt                                         <  1:3421150/91‑1 (MQ=255)

                                                    aaaaaaTACGTACATTACGTACACAAAGTGATTGGATTTCTAGTAAGCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGAtt                                        <  2:2589102/91‑1 (MQ=255)

                                                     aaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGAttt                                       >  1:125633/1‑91 (MQ=255)

                                                         aTACGTCCATTACGTAAACAAAGTGTTTGAATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtg                                   >  2:1005921/1‑91 (MQ=255)

                                                          tACGCCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtgt                                  >  1:3068358/1‑91 (MQ=255)

                                                           aCGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtgtg                                 >  1:3041918/1‑91 (MQ=255)

                                                            cGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATgtgt                                <  1:351272/91‑1 (MQ=255)

                                                            cGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATgtgt                                <  1:1761368/91‑1 (MQ=255)

                                                            cGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATgtgt                                >  1:2965576/1‑91 (MQ=255)

                                                            cGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATgtgt                                >  2:1640967/1‑91 (MQ=255)

                                                             gTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTc                               >  1:369907/1‑91 (MQ=255)

                                                                cATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCaaa                            <  2:2173525/91‑1 (MQ=255)

                                                                 tggaCGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAAt                           >  2:760851/4‑91 (MQ=255)

                                                                  ttACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATc                          <  2:146566/91‑1 (MQ=255)

                                                                  ttACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATc                          >  1:66613/1‑91 (MQ=255)

                                                                   tACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCa                         <  1:1982090/91‑1 (MQ=255)

                                                                      gTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAatt                      <  1:2173019/91‑1 (MQ=255)

                                                                         cacaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAattatt                   >  2:3800214/1‑91 (MQ=255)

                                                                         cacaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAattatt                   <  1:412313/91‑1 (MQ=255)

                                                                          acaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAAttatta                  <  1:3137483/91‑1 (MQ=255)

                                                                          acaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAAttatta                  <  1:1045669/91‑1 (MQ=255)

                                                                            aaaGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATtata                <  1:774933/91‑1 (MQ=255)

                                                                             aaGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTatat               <  2:791409/91‑1 (MQ=255)

                                                                              aGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATc              >  1:2835786/1‑91 (MQ=255)

                                                                              aGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATc              >  1:242072/1‑91 (MQ=255)

                                                                              aGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATc              >  1:2401765/1‑91 (MQ=255)

                                                                               gTGATTGGATTTCTAGTAATCCAATTTAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCa             <  1:1502537/91‑1 (MQ=255)

                                                                               gTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCa             <  1:160138/91‑1 (MQ=255)

                                                                               gTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCa             >  2:1059227/1‑91 (MQ=255)

                                                                               gTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCa             <  1:491434/91‑1 (MQ=255)

                                                                                 gATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCaaa           <  1:2076376/91‑1 (MQ=255)

                                                                                 gATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCaaa           >  2:1866246/1‑91 (MQ=255)

                                                                                 gATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCaaa           >  1:2987449/1‑91 (MQ=255)

                                                                                  aTTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAGGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAAc          <  1:3870677/91‑1 (MQ=255)

                                                                                   ttGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACa         >  2:1947258/1‑91 (MQ=255)

                                                                                   ttGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACa         <  2:1594883/91‑1 (MQ=255)

                                                                                   ttGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACa         >  2:1433636/1‑91 (MQ=255)

                                                                                    tGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAg        >  1:1627608/1‑91 (MQ=255)

                                                                                    tGGATTACTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAg        >  2:3263659/1‑91 (MQ=255)

                                                                                      gATTACTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGcc      >  1:287741/1‑91 (MQ=255)

                                                                                       aTTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCt     >  1:1038580/1‑91 (MQ=255)

                                                                                        tttCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCta    <  2:2688854/91‑1 (MQ=255)

                                                                                        tttCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCta    <  1:760851/91‑1 (MQ=255)

                                                                                        tttCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCta    <  1:3551826/91‑1 (MQ=255)

                                                                                         ttCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtat   <  2:2492138/91‑1 (MQ=255)

                                                                                          tCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtatt  <  1:1380068/91‑1 (MQ=255)

GCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATT  >  USA300TCH1516_ALE/1331146‑1331326

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