BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	110,660	T→A	7.0%	I173I (ATT→ATA)	USA300HOU_RS00515 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	110,660	0	T	A	7.0%	844.4 / 34.9	251	I173I (ATT→ATA)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (103/130); new base A (7/11); total (110/141)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

TGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAT  >  USA300TCH1516_ALE/110570‑110750

tGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAtt                                                                                            >  2:2710588/1‑91 (MQ=255)

 gAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                           <  2:2299454/91‑1 (MQ=255)

 gAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                           <  2:1657026/91‑1 (MQ=255)

    aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                        >  1:1331158/1‑91 (MQ=255)

    aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                        <  1:3864850/91‑1 (MQ=255)

     aaaTAAAAAAAGAAATCGAAAACTTGAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                       <  1:3602063/91‑1 (MQ=255)

      aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCatat                                                                                      <  2:2732301/91‑1 (MQ=255)

        tAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATaa                                                                                    <  2:237805/91‑1 (MQ=255)

         aaaaaaaGAAATCGTAAACTTTATGTTCTTTGTTCAATATGGCTACTTTATAAATTTGAAAAATTATAAAGACGGAGATATTTGATATAAt                                                                                   <  2:3611917/91‑1 (MQ=18)

           aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATcc                                                                                 >  2:2892990/1‑91 (MQ=255)

             aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                               <  2:1402768/91‑1 (MQ=255)

              aaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCaga                                                                              <  1:3721503/91‑1 (MQ=255)

               aGAAATCGAAAACTTTAAGTTCTCTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCCTATAATCCagag                                                                             >  2:3227652/1‑91 (MQ=255)

                gAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAgg                                                                            <  2:1997356/91‑1 (MQ=255)

                  aaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTG                                                                          <  1:2717042‑M1/91‑1 (MQ=255)

                   gtcGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGC                                                                         <  2:506136‑M1/90‑1 (MQ=255)

                   aTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGC                                                                         >  1:2752790‑M1/1‑91 (MQ=255)

                    tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCC                                                                        >  1:3647507‑M1/1‑91 (MQ=255)

                      gAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGA                                                                      >  1:31146‑M1/1‑91 (MQ=255)

                       aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                     >  2:397691/1‑91 (MQ=255)

                        aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                    <  2:1751601/91‑1 (MQ=255)

                         aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAata                                                                   <  1:3563096‑M1/91‑4 (MQ=255)

                         aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAata                                                                   <  1:2234178‑M1/91‑4 (MQ=255)

                         aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAata                                                                   <  2:3622905‑M1/91‑4 (MQ=255)

                         aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAAGTATAAAGACGGAGATATTTCATATAATCCCGAGGTGCCGAGtt                                                                   >  1:3877111/1‑91 (MQ=255)

                          aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAGAGACGGAGATATTTCATATAATCCAGAGGTGCCGAataa                                                                  <  2:2424606‑M1/91‑5 (MQ=255)

                          aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAataa                                                                  <  1:394171‑M1/91‑5 (MQ=255)

                          aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTa                                                                  <  1:3464042/91‑1 (MQ=255)

                           cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAataac                                                                 <  1:1795058‑M1/91‑6 (MQ=255)

                           cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                 <  1:2865928/91‑1 (MQ=255)

                            tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAataact                                                                <  1:3478935‑M1/91‑7 (MQ=255)

                            tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAataact                                                                <  1:1013113‑M1/91‑7 (MQ=255)

                            tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATt                                                                >  1:3296761‑M2/1‑90 (MQ=255)

                              tAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATtcg                                                              >  1:1756207‑M2/1‑88 (MQ=255)

                               aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGTAGATATTTCATATAATCCAGAGGTGCCGAGTTATtcgg                                                             >  2:3819147‑M2/1‑87 (MQ=255)

                                    cTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                        >  2:1495845/1‑91 (MQ=255)

                                      ttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                      <  1:2710588/91‑1 (MQ=255)

                                      ttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                      >  2:1219926/1‑91 (MQ=255)

                                       tGTTCAATATGGCGACTTTAAAAATTTGAAAAATTAGAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                     <  2:3062415/91‑1 (MQ=255)

                                         ttCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAATAACtcggcacctctgga                                                   >  1:2116498‑M2/1‑77 (MQ=255)

                                          tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                  <  2:2255800/91‑1 (MQ=255)

                                          tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                  <  2:2923282/91‑1 (MQ=255)

                                          tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAATAACtcggcacctctggat                                                  >  1:2093230‑M2/1‑76 (MQ=255)

                                           cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                 >  1:1470504/1‑91 (MQ=255)

                                           cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                 >  2:2895952/1‑91 (MQ=255)

                                               atGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                             <  1:1391940/91‑1 (MQ=255)

                                                  gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                          <  1:2154421/91‑1 (MQ=255)

                                                  gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                          <  1:1252325/91‑1 (MQ=255)

                                                    gACTTTATAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                        <  1:1932211/91‑1 (MQ=255)

                                                    gACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                        <  1:2371757/91‑1 (MQ=255)

                                                    gACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                        >  1:1612658/1‑91 (MQ=255)

                                                     aCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                       >  2:1062541/1‑91 (MQ=255)

                                                     aCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                       >  2:755149/1‑91 (MQ=255)

                                                      cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                      >  1:559072/1‑91 (MQ=255)

                                                      cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                      <  2:611114/91‑1 (MQ=255)

                                                      cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                      >  2:1886332/1‑91 (MQ=255)

                                                       tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                     <  2:3688810/91‑1 (MQ=255)

                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                  >  1:1466139/1‑91 (MQ=255)

                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                  >  2:1281798/1‑91 (MQ=255)

                                                           aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                 >  2:2526771/1‑91 (MQ=255)

                                                            aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                <  2:3563778/91‑1 (MQ=255)

                                                            aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                >  1:1075052/1‑91 (MQ=255)

                                                             aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                               >  1:313521/1‑91 (MQ=255)

                                                             aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                               >  2:3489384/1‑91 (MQ=255)

                                                                ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                            <  1:3753720/91‑1 (MQ=255)

                                                                 tGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGt                           >  1:2555599/1‑91 (MQ=255)

                                                                  gAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                          >  2:3425550/1‑91 (MQ=255)

                                                                  gAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                          >  1:1617131/1‑91 (MQ=255)

                                                                   aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                         >  1:1534533/1‑91 (MQ=255)

                                                                   aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                         >  2:877736/1‑91 (MQ=255)

                                                                   aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                         >  2:878402/1‑91 (MQ=255)

                                                                    aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAATAACtcggcacctctggattatatgaaatatctccgtctttataa                        >  2:1797127‑M2/1‑50 (MQ=255)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                       <  1:117832/91‑1 (MQ=40)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                       <  1:2475072/91‑1 (MQ=40)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                       >  2:2068869/1‑91 (MQ=40)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                       <  2:3277323/91‑1 (MQ=40)

                                                                     aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAATAACtcggcacctctggattatatgaaatatctccgtctttataat                       >  2:1119699‑M2/1‑49 (MQ=255)

                                                                      aaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAAc                      <  2:2588017/91‑1 (MQ=40)

                                                                       aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATtcggcacctctggattatatgatatatctccatctttataattc                     <  2:2336079‑M2/91‑45 (MQ=255)

                                                                       aTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                     <  2:2872737/91‑1 (MQ=40)

                                                                          ataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                  >  2:634390/1‑91 (MQ=39)

                                                                          ataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                  >  2:861643/1‑91 (MQ=39)

                                                                             aaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg               <  1:1666864/91‑1 (MQ=39)

                                                                              agatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc              <  2:3119307/87‑1 (MQ=9)

                                                                              aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc              <  2:2468225/91‑1 (MQ=12)

                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc              <  1:1508542/91‑1 (MQ=40)

                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc              >  1:816095/1‑91 (MQ=40)

                                                                               gatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa             <  2:2653400/88‑1 (MQ=16)

                                                                               gACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa             <  2:776328/91‑1 (MQ=255)

                                                                               gACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa             >  1:2396516/1‑91 (MQ=255)

                                                                                aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa            <  2:376074/91‑1 (MQ=255)

                                                                                  ggacatatatCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa          >  2:2127277/5‑91 (MQ=255)

                                                                                  ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                <  1:761933/85‑1 (MQ=40)

                                                                                  ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                >  2:761947/1‑85 (MQ=40)

                                                                                  ggAGATATTTCATATAATCCAGAGGTGCCGAATAACtcggcacctctggattatatgaaatatctccgtctttataatttttcaaattttt          >  2:928772‑M2/1‑36 (MQ=255)

                                                                                    agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga        <  1:96399/91‑1 (MQ=255)

                                                                                     gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat       >  1:602666/1‑91 (MQ=255)

                                                                                      atatatCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgata      >  1:2512323/1‑91 (MQ=255)

                                                                                        atTTCATATAATCCAGAGGTGCCGAGTTATtcggcacctctggattataagaaatatctccgtcattatacttccacaactctttaaagtt    <  2:565260‑M2/91‑62 (MQ=255)

                                                                                         tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatga   >  1:2032255/1‑91 (MQ=255)

TGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAT  >  USA300TCH1516_ALE/110570‑110750

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