BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	110,660	T→A	9.6%	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	9.6%	1513.0 / 119.7	429	I173I (ATT→ATA)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (163/222); new base A (24/17); total (188/241)
Fisher's exact test for biased strand distribution p-value = 6.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 443 total aligned reads displayed.

TGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAT  >  USA300TCH1516_ALE/110570‑110750

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

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

      aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCatat                                                                                      >  1:2359410/1‑91 (MQ=255)

       aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCAtata                                                                                     <  1:2074772/91‑1 (MQ=255)

        tAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATaa                                                                                    >  1:817808/1‑91 (MQ=255)

          aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                  >  2:292463/1‑91 (MQ=255)

          aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                  >  2:421528/1‑91 (MQ=255)

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

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

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

            aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCa                                                                                >  2:503809/1‑91 (MQ=255)

            aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCa                                                                                >  2:446779/1‑91 (MQ=255)

             aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                               >  1:2756178/1‑91 (MQ=255)

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

                 aaaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGt                                                                           <  1:1068492/91‑1 (MQ=255)

                 aaaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGt                                                                           <  2:5154899/91‑1 (MQ=255)

                    tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                        <  2:3570070/91‑1 (MQ=255)

                    tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                        >  1:1323854/1‑91 (MQ=255)

                     cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                       <  2:2071636/91‑1 (MQ=255)

                     cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                       >  1:3743479/1‑91 (MQ=255)

                       aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                     <  1:3204494/91‑1 (MQ=255)

                        aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                    >  1:2648486/1‑91 (MQ=255)

                         aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGaata                                                                   <  2:1503052/91‑4 (MQ=255)

                         aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGtt                                                                   >  2:2443107/1‑91 (MQ=255)

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

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

                              tAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCg                                                              <  1:3852073/91‑1 (MQ=255)

                              tAAGTTCTTTGCTCAATATGGCGACTTTAAAAATTTGAAAAGTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCg                                                              <  2:5642277/91‑1 (MQ=255)

                               aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                             <  1:446063/91‑1 (MQ=255)

                               aaGTTCTATGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCGTATAATCCAGAGGTGCCGAGTTATTCgg                                                             <  2:5397798/91‑1 (MQ=255)

                                aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                            >  2:857964/1‑91 (MQ=255)

                                aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                            <  1:1691099/91‑1 (MQ=255)

                                aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                            >  1:1773189/1‑91 (MQ=255)

                                  ttCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                          <  2:4320198/91‑1 (MQ=255)

                                   tCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                         <  2:4383335/91‑1 (MQ=255)

                                     tttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                       >  1:4911857/1‑91 (MQ=255)

                                     tttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                       >  1:4632352/1‑91 (MQ=255)

                                     tttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                       >  2:1137376/1‑91 (MQ=255)

                                     tttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                       <  1:3564477/91‑1 (MQ=255)

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

                                       tGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                     >  1:3440999/1‑91 (MQ=255)

                                            aaTATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                <  1:2265827/91‑1 (MQ=255)

                                              tatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                              <  1:5579907/91‑1 (MQ=255)

                                               atGGCGACTTTAGAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                             <  2:3622312/91‑1 (MQ=255)

                                               atGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTAATATAATCCGGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                             <  2:5446284/91‑1 (MQ=255)

                                                tGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                            <  1:1013101/91‑1 (MQ=255)

                                                tGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                            >  2:4917786/1‑91 (MQ=255)

                                                tGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                            <  1:2171809/91‑1 (MQ=255)

                                                 ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGGTAAATATCAATTAACTa                                           <  1:4272586/91‑1 (MQ=255)

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

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

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

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

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

                                                         tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAGTCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                   >  2:5224566/1‑91 (MQ=255)

                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                  <  1:4855390/91‑1 (MQ=255)

                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                  <  2:3201300/91‑1 (MQ=255)

                                                          aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                  <  1:197845/91‑1 (MQ=255)

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

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

                                                             aaTTTGAAAAATTATAACGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                               >  1:2691368/1‑91 (MQ=255)

                                                             aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATACAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                               >  2:3484368/1‑91 (MQ=255)

                                                                ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                            >  2:2746182/1‑91 (MQ=255)

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

                                                                  gAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATCATGTa                          >  2:5306048/1‑91 (MQ=255)

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

                                                                   aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                         <  1:2534846/91‑1 (MQ=255)

                                                                    aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                        >  2:3320743/1‑91 (MQ=39)

                                                                    aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                        <  2:3756351/91‑1 (MQ=40)

                                                                    aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                        >  1:5226872/1‑91 (MQ=40)

                                                                     aaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                       <  2:2366144/91‑1 (MQ=12)

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

                                                                         tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                   <  1:922513/91‑1 (MQ=39)

                                                                           taAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtga                 >  2:772535/1‑89 (MQ=37)

                                                                            aaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTat                >  1:2905939/1‑90 (MQ=38)

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

                                                                              aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc              <  1:1487824/91‑1 (MQ=38)

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

                                                                               gatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa             >  2:1956926/4‑91 (MQ=12)

                                                                               gatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa             >  1:2449935/4‑91 (MQ=255)

                                                                               gatggACATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa             <  1:4716341/88‑1 (MQ=255)

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

                                                                               gACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa             <  2:2546342/91‑1 (MQ=39)

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

                                                                                  ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGACTATAATGTAAAACAATTACGAaaaa          <  2:643411/91‑1 (MQ=255)

                                                                                  ggAGATATCTCATACGATCTAGAGATGCCGAGTTATTCGACTAAATATCAATTAACTAATGATGATTATAATGTATATCAATTACGCaaaa          <  2:1817876/91‑1 (MQ=11)

                                                                                  ggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa          >  2:2832170/1‑91 (MQ=255)

                                                                                  ggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa          <  2:1337561/91‑1 (MQ=255)

                                                                                   gagaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                     >  2:3401207/1‑79 (MQ=40)

                                                                                   gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAg         >  1:566428/1‑91 (MQ=255)

                                                                                     gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat       <  2:2775215/91‑1 (MQ=255)

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

                                                                                     gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat       <  2:2247916/91‑1 (MQ=255)

                                                                                     gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat       <  2:1465348/91‑1 (MQ=255)

                                                                                     gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat       >  2:1417051/1‑91 (MQ=25)

                                                                                      atatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgata      <  2:483616/91‑1 (MQ=255)

                                                                                      atatTTCATATAATCCAGAGATGCCGAGTTATTCAGCTAAATATCAATTAACTAATGATGATTATCATGTAAAACAATTACGCAAAAgata      <  2:4732446/91‑1 (MQ=14)

                                                                                       tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatat     >  2:2173909/1‑91 (MQ=255)

                                                                                       tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatat     >  1:959745/1‑91 (MQ=255)

                                                                                         tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatga   <  1:1875348/91‑1 (MQ=255)

TGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAT  >  USA300TCH1516_ALE/110570‑110750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`