BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,172,078	T→A	14.7%	F43I (TTT→ATT)	leuB →	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,172,078	0	T	A	14.7%	326.7 / 40.7	108	F43I (TTT→ATT)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (39/53); new base A (8/8); total (47/61)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 115 total aligned reads displayed.

ATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTG  >  USA300TCH1516_ALE/2171988‑2172168

aTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAt                                                                                             >  1:1463826/1‑91 (MQ=255)

aTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAt                                                                                             <  1:572963/91‑1 (MQ=255)

 tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAtt                                                                                            <  1:3103461/91‑1 (MQ=255)

 tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAtt                                                                                            <  1:3242057/91‑1 (MQ=255)

   ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATCATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAatt                                                                                           <  1:1040654/91‑1 (MQ=255)

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

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

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

        aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTggtggt                                                                                     <  1:3243723/91‑1 (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)

           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)

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

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

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

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

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

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

                           ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCg                                                                 >  1:3038080/1‑91 (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:1677327/91‑1 (MQ=255)

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

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

                                   gCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCt                                                         <  1:2180825/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:2663713/91‑1 (MQ=255)

                                        aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                    <  1:3661968/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:3780467/1‑91 (MQ=255)

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

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

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

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

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

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

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

                                                 gTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAAc                                           <  1:2942860/91‑1 (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:623413/1‑91 (MQ=255)

                                                     taaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTa                                       <  1:2614547/91‑1 (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:886659/91‑1 (MQ=255)

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

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

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

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

                                                                ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                            >  1:3778428/1‑91 (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:3141185/1‑91 (MQ=255)

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

                                                                            tAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                <  1:2850637/91‑1 (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: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: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)

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

                                                                                  caccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt           <  1:3907323/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)

                                                                                         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)

ATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTG  >  USA300TCH1516_ALE/2171988‑2172168

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