BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,172,088	Δ1 bp	10.9%	coding (137/1047 nt)	leuB →	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,172,088	0	C	.	10.9%	356.3 / 37.7	110	coding (137/1047 nt)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand): ref base C (39/59); new base . (9/3); total (48/62)
Fisher's exact test for biased strand distribution p-value = 2.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 115 total aligned reads displayed.

AAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTG  >  USA300TCH1516_ALE/2171998‑2172177

aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGc                                                                                            >  1:3560500/1‑91 (MQ=255)

aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGc                                                                                            <  1:2164325/91‑1 (MQ=255)

 aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGc                                                                                            <  1:1462916/91‑1 (MQ=37)

 aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGcc                                                                                           <  1:549606/91‑1 (MQ=255)

    tttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGctct                                                                                         <  1:106875/91‑4 (MQ=25)

       gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑tt                                                                                    >  1:2505624/1‑91 (MQ=255)

       gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑tt                                                                                    >  1:459321/1‑91 (MQ=255)

       gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATC‑AATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑ttt                                                                                   <  1:1150262/91‑2 (MQ=18)

        aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAAGTTTAATTATC‑AATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑tttg                                                                                  <  1:3290348/91‑3 (MQ=12)

        aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑ttt                                                                                   <  1:377702/91‑2 (MQ=21)

        aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTg                                                                                   <  1:3339601/91‑1 (MQ=255)

         aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGa                                                                                  >  1:3820329/1‑91 (MQ=255)

         aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGa                                                                                  <  1:2555488/91‑1 (MQ=255)

          cGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGAt                                                                                 <  1:608925/91‑1 (MQ=255)

            gATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATAc                                                                               <  1:1073467/91‑1 (MQ=255)

                 ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCg                                                                          >  1:3038080/1‑91 (MQ=255)

                  tATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCgg                                                                         <  1:2067208/91‑1 (MQ=255)

                    tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCg                                                                       <  1:2489723/91‑1 (MQ=255)

                      aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAg                                                                     <  1:3091691/91‑1 (MQ=255)

                       ttGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGCGCCTCTA‑TTGATACATTCGGCGAGc                                                                    >  1:3874071/1‑91 (MQ=255)

                        tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGcc                                                                   <  1:928159/91‑1 (MQ=255)

                        tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGcc                                                                   <  1:1677327/91‑1 (MQ=255)

                         gCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCt                                                                  <  1:2180709/91‑1 (MQ=255)

                            tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTCCCTCTA‑TTGATACATTCGGCGAGCCTTTa                                                               >  1:3323587/1‑91 (MQ=255)

                              aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                             <  1:3661968/91‑1 (MQ=255)

                              aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                             <  1:2663713/91‑1 (MQ=255)

                               aaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACt                                                            <  1:403733/91‑1 (MQ=255)

                                  tataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgag                                                         <  1:393516/91‑1 (MQ=255)

                                   ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                                        <  1:3220523/91‑1 (MQ=255)

                                   ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                                        >  1:341749/1‑91 (MQ=255)

                                   ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                                        >  1:3780467/1‑91 (MQ=255)

                                   ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                                        >  1:491827/1‑91 (MQ=255)

                                   ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                                        >  1:1066213/1‑91 (MQ=255)

                                   ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                                        >  1:1065216/1‑91 (MQ=255)

                                    taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                                       <  1:2701040/91‑1 (MQ=255)

                                    taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                                       <  1:2700864/91‑1 (MQ=255)

                                    taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                                       >  1:2616744/1‑91 (MQ=255)

                                    taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                                       >  1:1275318/1‑91 (MQ=255)

                                     aaGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaaa                                                      >  1:162105/1‑91 (MQ=255)

                                       gTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAAc                                                    <  1:1773807/91‑1 (MQ=255)

                                       gTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAAc                                                    <  1:1774288/91‑1 (MQ=255)

                                        taataaATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAAcc                                                   <  1:2609495/91‑1 (MQ=255)

                                          ataaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCtt                                                 >  1:622479/1‑91 (MQ=255)

                                          ataaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCtt                                                 >  1:623413/1‑91 (MQ=255)

                                            aaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCttca                                               >  1:3999599/1‑89 (MQ=255)

                                             aataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTaaa                                              <  1:900614/91‑1 (MQ=255)

                                              atatAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAAt                                             >  1:1101165/1‑91 (MQ=255)

                                                ataACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGc                                           <  1:2559663/91‑1 (MQ=255)

                                                 taACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                          <  1:1438783/91‑1 (MQ=255)

                                                 taACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                          <  1:886659/91‑1 (MQ=255)

                                                     tttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTa                                      >  1:2872003/1‑91 (MQ=255)

                                                      ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                                     >  1:3778428/1‑91 (MQ=255)

                                                      ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                                     <  1:1560519/91‑1 (MQ=255)

                                                      ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                                     <  1:1370715/91‑1 (MQ=255)

                                                          ttATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAaga                                 <  1:1436316/91‑1 (MQ=255)

                                                             tCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa                              >  1:7166/1‑91 (MQ=255)

                                                             atcaaTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc                               >  1:2032531/4‑91 (MQ=17)

                                                              cAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                             >  1:3141185/1‑91 (MQ=255)

                                                               caaTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                             >  1:238784/2‑91 (MQ=18)

                                                               aaaTAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTACAAGAGCAGa                            >  1:3973077/1‑91 (MQ=255)

                                                                 aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt                           >  1:3138144/1‑91 (MQ=18)

                                                                 aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt                           >  1:3138656/1‑91 (MQ=18)

                                                                 aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt                           >  1:1944061/1‑91 (MQ=18)

                                                                  tAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                         <  1:2850637/91‑1 (MQ=255)

                                                                  tAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                         >  1:1724508/1‑91 (MQ=255)

                                                                   aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCt                        <  1:1949898/91‑1 (MQ=255)

                                                                   aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCt                        <  1:3667850/91‑1 (MQ=255)

                                                                   aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCt                        <  1:3668319/91‑1 (MQ=255)

                                                                   aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAAAAAACCTTAAATGCGTGTAAAATAGCATCTGCt                        <  1:1877348/91‑1 (MQ=255)

                                                                   aga‑gCACCACGAATTTGGTGGTGCCCCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGATAACCTTAAATGCGTGTAAAATAGCATCTGCt                        <  1:3653883/91‑1 (MQ=255)

                                                                    ga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa                       <  1:1918169/91‑1 (MQ=255)

                                                                     aggCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa                       >  1:2116223/3‑91 (MQ=21)

                                                                     aggCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa                       >  1:2115604/3‑91 (MQ=21)

                                                                     a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGTGTGTAAAAGAGCAGATGCTAt                      <  1:1689832/91‑1 (MQ=255)

                                                                     a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt                      >  1:136700/1‑91 (MQ=255)

                                                                     a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt                      >  1:136574/1‑91 (MQ=255)

                                                                     a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt                      <  1:1257914/91‑1 (MQ=255)

                                                                      ggCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt                      >  1:592848/2‑91 (MQ=25)

                                                                       gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                     >  1:3836259/1‑91 (MQ=255)

                                                                       gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                     >  1:3836762/1‑91 (MQ=255)

                                                                        caccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt                    <  1:3707147/91‑1 (MQ=255)

                                                                        caccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt                    <  1:1533814/91‑1 (MQ=255)

                                                                            acGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCATCTGCTCTTTTAca                <  1:48418/91‑2 (MQ=38)

                                                                            acGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACt                <  1:2404351/91‑1 (MQ=255)

                                                                             cGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTg               >  1:298515/1‑91 (MQ=255)

                                                                               aaTTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTggg             >  1:1643235/1‑91 (MQ=255)

                                                                               aaTTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTggg             >  1:1807093/1‑91 (MQ=255)

                                                                               aaTTTGGAGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAATCCTTAAATGCGTGTAAAAGAGCAGATGCTAGTTTACTggg             <  1:2087089/91‑1 (MQ=255)

                                                                                 tttGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg           <  1:1553480/91‑1 (MQ=255)

                                                                                   tggtggTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCa         <  1:145573/91‑1 (MQ=255)

                                                                                    ggtggtGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa        <  1:512425/91‑1 (MQ=255)

                                                                                    ggtggtGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa        <  1:1150378/91‑1 (MQ=255)

                                                                                    cgtggtg‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa        >  1:281488/2‑91 (MQ=255)

                                                                                      tggtgCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATc      >  1:4087233/1‑91 (MQ=255)

                                                                                      tggtgCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATc      >  1:4087201/1‑91 (MQ=255)

                                                                                      tggtgCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATc      <  1:3807299/91‑1 (MQ=255)

                                                                                       ggtgCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCg     >  1:197347/1‑91 (MQ=255)

                                                                                         tgCCTCTA‑TTGATACATTCGGCGAGCCTTTAACGGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGt   <  1:353587/91‑1 (MQ=255)

                                                                                          gCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTg  <  1:4028384/91‑1 (MQ=255)

                                                                                          gCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTg  <  1:1943983/91‑1 (MQ=255)

AAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTG  >  USA300TCH1516_ALE/2171998‑2172177

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