BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,172,081	G→C	12.7%	G44R (GGT→CGT)	leuB →	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,172,081	0	G	C	12.7%	322.5 / 33.7	102	G44R (GGT→CGT)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (36/53); new base C (8/5); total (44/58)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 113 total aligned reads displayed.

GTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAA  >  USA300TCH1516_ALE/2171992‑2172171

gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTtgg                                                                                            >  1:118502/1‑91 (MQ=255)

  ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATCAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTtggtg                                                                                          >  1:3589313/1‑91 (MQ=255)

  ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTtggtg                                                                                          >  1:2937679/1‑91 (MQ=255)

     gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                        <  1:3363351/91‑1 (MQ=38)

     gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                        <  1:1881879/91‑1 (MQ=38)

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

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

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

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

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

             gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑tt                                                                              >  1:2505624/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:2555488/91‑1 (MQ=255)

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

               aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGa                                                                            <  1:2053230/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:928214/91‑1 (MQ=255)

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

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

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

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

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

                                    aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                       <  1:3661968/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: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:1275318/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)

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

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

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

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

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

                                                ataaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCtt                                           >  1:622479/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:3778499/1‑91 (MQ=255)

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

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

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

                                                                 tATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag                          >  1:4111202/1‑91 (MQ=255)

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

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

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

                                                                    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:3138656/1‑91 (MQ=18)

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

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

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

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

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

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

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

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

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

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

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

                                                                           a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt                >  1:136700/1‑91 (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:3907323/91‑1 (MQ=255)

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

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

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

                                                                               accacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt             <  1:3611866/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)

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

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

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

                                                                                         tggtggTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCa   <  1:144516/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)

GTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAA  >  USA300TCH1516_ALE/2171992‑2172171

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