BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,331,241	T→A	9.5%	S190R (AGT→AGA)	USA300HOU_RS06555 →	Monoacylglycerol lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,331,241	0	T	A	9.5%	772.4 / 56.2	241	S190R (AGT→AGA)	USA300HOU_RS06555	Monoacylglycerol lipase
Reads supporting (aligned to +/- strand): ref base T (120/98); new base A (16/7); total (136/105)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 248 total aligned reads displayed.

GAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATG  >  USA300TCH1516_ALE/1331151‑1331331

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

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

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

   aCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAGAAAt                                                                                         <  1:495745/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

        ttGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAGAAATCCAAt                                                                                    <  1:1257028/91‑1 (MQ=255)

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

         tGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAAtt                                                                                   >  2:732837/1‑91 (MQ=255)

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

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

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

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

             tGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAg                                                                               >  2:3284819/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

                         ttAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGAATTCTAGTAATCCAATTGAAGTAGATAACTtta                                                                   <  2:3808913/91‑1 (MQ=255)

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

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

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

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

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

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

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

                                      gTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCata                                                      >  2:563402/1‑91 (MQ=255)

                                         ttAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAg                                                   >  2:3698471/1‑91 (MQ=255)

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

                                          taataaAAAAATACGTCCATTACGTACACAATGTGATTGGAATTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGt                                                  >  2:1563611/1‑91 (MQ=255)

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

                                          taataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATATCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGt                                                  >  2:2762604/1‑91 (MQ=255)

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

                                               aaaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGAtt                                             >  1:3698530/1‑91 (MQ=255)

                                               aaaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGAtt                                             <  2:1639118/91‑1 (MQ=255)

                                                  aaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTaa                                          <  2:2505502/91‑1 (MQ=255)

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

                                                   aaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAt                                         >  1:2944473/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)

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

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

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

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

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

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

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

                                                              tACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCa                              <  1:1982090/91‑1 (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:774713/91‑1 (MQ=255)

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

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

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

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

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

                                                                             aTTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAAc               >  2:1380068/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                  attacTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCt          >  1:2729309/5‑91 (MQ=255)

                                                                                  aTTTCTAGTACTCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCt          <  2:2129356/91‑1 (MQ=255)

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

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

                                                                                   ttacTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCta         >  1:825863/4‑91 (MQ=255)

                                                                                    ttCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGTATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtat        <  2:3196534/91‑1 (MQ=255)

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

                                                                                     acTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtatt       >  1:2491920/2‑91 (MQ=255)

                                                                                     acTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGCGTCAAATCAATTATTATATCAAACAGCCtatt       >  2:3008874/2‑91 (MQ=255)

                                                                                      cTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtatta      >  2:72724/1‑91 (MQ=255)

                                                                                       tAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACCGCCtattat     >  1:3719829/1‑91 (MQ=255)

                                                                                       tAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtattat     <  2:919105/91‑1 (MQ=255)

                                                                                       tAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtattat     <  2:1794007/91‑1 (MQ=255)

                                                                                       tAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtattat     <  1:2164329/91‑1 (MQ=255)

                                                                                       tAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtattat     <  2:2920480/91‑1 (MQ=255)

                                                                                        aGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATtata    >  1:362748/1‑91 (MQ=255)

                                                                                         gTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTatat   <  1:2586819/91‑1 (MQ=255)

                                                                                          tAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATg  >  2:3736158/1‑91 (MQ=255)

GAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATG  >  USA300TCH1516_ALE/1331151‑1331331

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