BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,331,236	T→A	7.1%	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	7.1%	876.9 / 38.2	262	S189T (TCT→ACT)	USA300HOU_RS06555	Monoacylglycerol lipase
Reads supporting (aligned to +/- strand): ref base T (125/118); new base A (10/9); total (135/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

GCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATT  >  USA300TCH1516_ALE/1331146‑1331326

 cTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTc                                                                                           <  2:320864/91‑1 (MQ=255)

 cTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTc                                                                                           >  2:2302446/1‑91 (MQ=255)

     ggaaCGAGTTGAAGGGGTTAACCAGGTAGTATCAAATAGTTTGAATAAAAAAATACGTCCATTACGTACACAAAGTGATGGGATTTCTAGt                                                                                       <  1:1870538/89‑1 (MQ=255)

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

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

     gaAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAga                                                                                       <  2:2186556/91‑2 (MQ=255)

       aaCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTaa                                                                                     >  1:2556590/1‑91 (MQ=255)

        aCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAAt                                                                                    >  1:3028562/1‑91 (MQ=255)

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

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

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

           aGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAGAAATCCa                                                                                 <  2:1084681/91‑1 (MQ=255)

            gTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCaa                                                                                >  1:1163587/1‑91 (MQ=255)

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

              tGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAAtt                                                                              <  1:1853436/91‑1 (MQ=255)

               gAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTg                                                                             >  2:857095/1‑91 (MQ=255)

               gAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTg                                                                             >  1:2742095/1‑91 (MQ=255)

                aaTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGa                                                                            <  1:1119073/91‑1 (MQ=255)

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

                aaTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGa                                                                            <  1:227068/91‑1 (MQ=255)

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

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

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

                   gggTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGt                                                                         >  1:602866/1‑91 (MQ=255)

                   gggTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGt                                                                         <  1:1303333/91‑1 (MQ=255)

                     gTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAg                                                                       <  1:3287559/91‑1 (MQ=255)

                     gTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAg                                                                       >  1:1038599/1‑91 (MQ=255)

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

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

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

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

                             gTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACttt                                                               >  1:208759/1‑91 (MQ=255)

                               tAGTATCAAATGGTTTTAATAAATAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAGTAATCCAATTGAAGTAGATAACTttat                                                             <  2:3022571/91‑1 (MQ=255)

                               tAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTttat                                                             >  2:3121797/1‑91 (MQ=255)

                                  tATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTaa                                                          >  1:491258/1‑91 (MQ=255)

                                   aTCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTaaa                                                         >  1:957037/1‑91 (MQ=255)

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

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

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

                                       aaTAGTTATAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATc                                                     <  2:3192853/91‑1 (MQ=255)

                                          aGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCat                                                  <  2:1215039/91‑1 (MQ=255)

                                           gTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCata                                                 <  2:2226330/91‑1 (MQ=255)

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

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

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

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

                                               taataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTCATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGt                                             >  2:3478664/1‑91 (MQ=255)

                                               taataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGt                                             >  2:1932958/1‑91 (MQ=255)

                                               taataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGt                                             >  2:3844465/1‑91 (MQ=255)

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

                                               aaaaaaaaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGt                                             <  2:1277499/87‑1 (MQ=255)

                                                aataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTg                                            >  2:539648/1‑91 (MQ=255)

                                                aataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTg                                            >  2:2433929/1‑91 (MQ=255)

                                                aataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTg                                            >  2:2681782/1‑91 (MQ=255)

                                                 taaaaaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAGTTGAAGTAGATAACTTTATTAAAGATCCATATAGTgg                                           <  2:1429247/89‑1 (MQ=255)

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

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

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

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

                                                        aaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAt                                    <  2:2175431/91‑1 (MQ=255)

                                                        aaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAt                                    >  1:1366347/1‑91 (MQ=255)

                                                         aTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtg                                   >  1:2855655/1‑91 (MQ=255)

                                                           aCGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtgtg                                 <  2:3062304/91‑1 (MQ=255)

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

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

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

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

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

                                                              tCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCa                              >  2:1525764/1‑91 (MQ=255)

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

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

                                                                    aCGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCaa                        >  1:3019035/1‑91 (MQ=255)

                                                                    aCGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCaa                        >  1:918295/1‑91 (MQ=255)

                                                                     cGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAat                       >  1:1448526/1‑91 (MQ=255)

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

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

                                                                       tACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAatta                     >  1:3099736/1‑91 (MQ=255)

                                                                        acacAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAattat                    >  1:952231/1‑91 (MQ=255)

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

                                                                          acaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATTAAttatta                  <  2:2639542/91‑1 (MQ=255)

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

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

                                                                             aaGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATGTAGTGGATTTAATGTGTCAAATCAATTATTatat               <  2:335948/91‑1 (MQ=255)

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

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

                                                                                tGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCaa            >  2:1460999/1‑91 (MQ=255)

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

                                                                                 gATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATGTAATGTGTCAAATCAATTATTATATCaaa           <  2:3448274/91‑1 (MQ=255)

                                                                                   ttGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATTAAACa         >  2:2551101/1‑91 (MQ=255)

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

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

                                                                                    tGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAg        <  2:2005797/91‑1 (MQ=255)

                                                                                     ggATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGc       >  1:2951126/1‑91 (MQ=255)

                                                                                     ggATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGc       >  2:2837910/1‑91 (MQ=255)

                                                                                     ggATTACTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGc       >  1:3624247/1‑91 (MQ=255)

                                                                                     ggATTACTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGc       >  2:953970/1‑91 (MQ=255)

                                                                                      gATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGcc      <  1:473066/91‑1 (MQ=255)

                                                                                      gATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGcc      >  1:3395039/1‑91 (MQ=255)

                                                                                      gATTTCTAGGAATCAAATTGAAGTAGATAGCTGTATGAAAGATCCATATAGTGGATTGAATGTGTCAAATCAATTATTACATCAAACAGcc      <  2:1381633/91‑1 (MQ=255)

                                                                                       aTTACTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATGCAAATAGTGGATTTAATGTGTCAAATCAAGTATTATATCAAACAGccc     >  2:26084/1‑90 (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: `‑`