BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,331,241	T→A	7.9%	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	7.9%	897.2 / 50.6	266	S190R (AGT→AGA)	USA300HOU_RS06555	Monoacylglycerol lipase
Reads supporting (aligned to +/- strand): ref base T (126/117); new base A (17/4); total (144/122)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 277 total aligned reads displayed.

GAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATG  >  USA300TCH1516_ALE/1331151‑1331331

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

   aCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAAt                                                                                         <  2:1125602/91‑1 (MQ=255)

   aCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAAt                                                                                         <  2:124128/91‑1 (MQ=255)

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

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

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

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

       gTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCaa                                                                                     <  1:3355379/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

            aTGGGTTAACCAGTTAGTGTCAGATAGTTTCAATAAAAAACTACGTCCTTTACGGACAGCAAGTGATTGGATATCTAGTACACCAATCGag                                                                                >  2:2046312/1‑90 (MQ=255)

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

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

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

               ggTTAACCAGTTAGTAGCAAATAGATTTAATAAAAAAATACTTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTa                                                                             <  2:3475697/91‑1 (MQ=255)

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

                     ccAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACAGAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAAc                                                                       <  1:3423533/91‑1 (MQ=255)

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

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

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

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

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

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

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

                                       ttttAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCatat                                                     <  1:2516668/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                            ataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTgg                                                <  1:1480596/91‑1 (MQ=255)

                                            ataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTgg                                                <  2:710261/91‑1 (MQ=255)

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

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

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

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

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

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

                                                      acttCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtgtg                                      >  1:696453/4‑91 (MQ=255)

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

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

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

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

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

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

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

                                                             ttACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTGGATAACTTCATTAAAGATCCATATAGTGGATTTAATGTGTCAAATc                               <  1:3158056/91‑1 (MQ=255)

                                                              tACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCa                              >  2:3237812/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)

                                                                   acacAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAattat                         <  2:3124554/91‑1 (MQ=255)

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

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

                                                                     acaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAAttatta                       >  1:1807229/1‑91 (MQ=255)

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

                                                                     acaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAAttatta                       >  2:3647507/1‑91 (MQ=255)

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

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

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

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

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

                                                                             aTTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAAc               <  1:1909009/91‑1 (MQ=255)

                                                                             aTTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAAc               <  2:3722676/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                ggATTACTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATCTAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGc            >  2:577578/1‑91 (MQ=255)

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

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

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

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

                                                                                  attacTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCt          >  2:2277348/5‑91 (MQ=255)

                                                                                  aTTACTAGAAATCCAATTGAAGTCGATAACTTTATTAAAGATCCATATAGTGGATTTAATGCGTCAAATCAATCATTATATCAAACAGCCt          >  2:25926/1‑91 (MQ=255)

                                                                                   tttCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAAGAGCCta         <  1:3078791/91‑1 (MQ=255)

                                                                                    tacTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtat        >  2:1903825/3‑91 (MQ=255)

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

                                                                                      cTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtatta      >  1:3787696/1‑91 (MQ=255)

                                                                                        aGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATtata    <  1:3570827/91‑1 (MQ=255)

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

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

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

                                                                                          tAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATg  <  2:1568284/91‑1 (MQ=255)

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

GAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATG  >  USA300TCH1516_ALE/1331151‑1331331

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