BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,744,362	G→A	57.0%	H104Y (CAC→TAC)	apt ←	Adenine phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,744,362	0	G	A	57.0%	47.1 / 127.0	100	H104Y (CAC→TAC)	apt	Adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand): ref base G (19/24); new base A (24/33); total (43/57)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAA  >  USA300TCH1516_ALE/1744272‑1744451

cTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTtgtg                                                                                           >  1:3517069/1‑91 (MQ=255)

cTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTtgtg                                                                                           <  1:1424833/91‑1 (MQ=255)

    ttttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACAtt                                                                                       >  1:1760362/1‑91 (MQ=255)

      ttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGt                                                                                     <  1:3036029/91‑1 (MQ=255)

      ttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGt                                                                                     <  1:3724199/91‑1 (MQ=255)

        aTTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTa                                                                                   >  1:2908551/1‑91 (MQ=255)

         ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGttca                                                                                  >  1:186299/1‑89 (MQ=255)

         ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaa                                                                                  <  1:1153465/91‑1 (MQ=255)

         ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTaa                                                                                  <  1:464836/91‑1 (MQ=255)

          tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTACCACACGTTGACCTGGTTTAATTGCATCTTCGTACATTGTTaaa                                                                                 >  1:4152364/1‑91 (MQ=255)

          tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaaa                                                                                 <  1:3471855/91‑1 (MQ=255)

          tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTaaa                                                                                 >  1:77569/1‑91 (MQ=255)

            ctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAAc                                                                               >  1:2122244/1‑91 (MQ=255)

             tgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACa                                                                              <  1:1628963/91‑1 (MQ=255)

                ttCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAttt                                                                           >  1:2589559/1‑91 (MQ=255)

                ttCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAttt                                                                           >  1:857569/1‑91 (MQ=255)

                     tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTAcc                                                                      >  1:2529458/1‑91 (MQ=255)

                     tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTAcc                                                                      >  1:254293/1‑91 (MQ=255)

                     tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTAcc                                                                      >  1:2820175/1‑91 (MQ=255)

                     tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAcc                                                                      <  1:595973/91‑1 (MQ=255)

                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCa                                                                     <  1:786007/91‑1 (MQ=255)

                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCa                                                                     >  1:1392089/1‑91 (MQ=255)

                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCa                                                                     >  1:3070709/1‑91 (MQ=255)

                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCa                                                                     <  1:3655928/91‑1 (MQ=255)

                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGAACCa                                                                     >  1:3460411/1‑91 (MQ=255)

                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCa                                                                     <  1:750296/91‑1 (MQ=255)

                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCa                                                                     <  1:1021407/91‑1 (MQ=255)

                         accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCatat                                                                  >  1:2770553/1‑91 (MQ=255)

                         accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                  <  1:2594552/91‑1 (MQ=255)

                         accaccAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                  <  1:1049007/91‑1 (MQ=255)

                          ccaccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATAtt                                                                 >  1:1772286/1‑91 (MQ=255)

                            accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCt                                                               >  1:1355657/1‑91 (MQ=255)

                            accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCt                                                               >  1:2125863/1‑91 (MQ=255)

                            accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCt                                                               >  1:2125934/1‑91 (MQ=255)

                             ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTa                                                              >  1:3629497/1‑91 (MQ=255)

                               aGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAgg                                                            <  1:1395503/91‑1 (MQ=255)

                               aGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAgg                                                            <  1:1812260/91‑1 (MQ=255)

                                gTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGt                                                           >  1:1133084/1‑91 (MQ=255)

                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCa                                                         >  1:207907/1‑91 (MQ=255)

                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCa                                                         <  1:124989/91‑1 (MQ=255)

                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCa                                                         <  1:3259707/91‑1 (MQ=255)

                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCa                                                         >  1:207941/1‑91 (MQ=255)

                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGATCa                                                         >  1:3211907/1‑91 (MQ=255)

                                   gCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCAt                                                        <  1:2664348/91‑1 (MQ=255)

                                    cTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATa                                                       >  1:1249971/1‑91 (MQ=255)

                                     taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATAc                                                      <  1:2202066/91‑1 (MQ=255)

                                     taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATAc                                                      <  1:1665446/91‑1 (MQ=255)

                                      aataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACt                                                     >  1:1797883/1‑91 (MQ=255)

                                      aataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACt                                                     <  1:3837425/91‑1 (MQ=255)

                                       ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTc                                                    <  1:2581324/91‑1 (MQ=255)

                                        taaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATGTGTACCATATTCTAGGTCATACTCa                                                   >  1:3853396/1‑91 (MQ=255)

                                         aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCAt                                                  >  1:2638189/1‑91 (MQ=255)

                                             catcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACg                                              >  1:2163766/1‑91 (MQ=255)

                                               tcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGaa                                            >  1:272460/1‑91 (MQ=255)

                                               tcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGaa                                            <  1:1660238/91‑1 (MQ=255)

                                               tcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGaa                                            <  1:783150/91‑1 (MQ=255)

                                               tcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGaa                                            <  1:1931130/91‑1 (MQ=255)

                                                 tGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATg                                          >  1:519756/1‑91 (MQ=255)

                                                   tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGAc                                        <  1:2250236/91‑1 (MQ=255)

                                                   tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGAc                                        >  1:1549058/1‑91 (MQ=255)

                                                       tAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCa                                    >  1:2175850/1‑91 (MQ=255)

                                                        aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCAc                                   <  1:1502433/91‑1 (MQ=255)

                                                         acacacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACg                                  >  1:3487277/1‑91 (MQ=255)

                                                          cacacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGa                                 <  1:404587/91‑1 (MQ=255)

                                                           acacGTTGACCTGGTTTAGTTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAg                                <  1:3807571/91‑1 (MQ=255)

                                                           acacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAg                                <  1:3589853/91‑1 (MQ=255)

                                                           acacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAg                                <  1:3016877/91‑1 (MQ=255)

                                                             acGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGt                              <  1:1532659/91‑1 (MQ=255)

                                                               gTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTaa                            <  1:1174627/91‑1 (MQ=255)

                                                               gTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTaa                            <  1:465660/91‑1 (MQ=255)

                                                               gTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTaa                            <  1:588521/91‑1 (MQ=255)

                                                                 tGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAGACATTTGTACCATATTCTAGGTCATACTCATAGCGAATGACTTCACGAGGTGAtt                          <  1:1769329/91‑1 (MQ=255)

                                                                 tGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAtt                          <  1:1171314/91‑1 (MQ=255)

                                                                 tGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAtt                          <  1:1841449/91‑1 (MQ=255)

                                                                 tGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAtt                          <  1:3781806/91‑1 (MQ=255)

                                                                   aCCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTc                        <  1:733155/91‑1 (MQ=255)

                                                                    ccTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                       >  1:3586076/1‑91 (MQ=255)

                                                                    ccTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                       >  1:556444/1‑91 (MQ=255)

                                                                    ccTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                       >  1:3586030/1‑91 (MQ=255)

                                                                    ccTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                       >  1:4019772/1‑91 (MQ=255)

                                                                    ccTGGTTTAATTGCAGCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                       <  1:507464/91‑1 (MQ=255)

                                                                     cTGGTTTAATTGCATCTTTGTATATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTccc                      <  1:3900208/91‑1 (MQ=255)

                                                                     cTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTccc                      <  1:320391/91‑1 (MQ=255)

                                                                      tGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTGCCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCct                     >  1:2438285/1‑91 (MQ=255)

                                                                       ggTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCctt                    <  1:1455826/91‑1 (MQ=255)

                                                                       ggTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCctt                    <  1:1455379/91‑1 (MQ=255)

                                                                        gTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttc                   >  1:722550/1‑91 (MQ=255)

                                                                          ttAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttctt                 <  1:2755419/91‑1 (MQ=255)

                                                                            aaTTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt               <  1:3217738/91‑1 (MQ=255)

                                                                            aaTTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt               <  1:2292116/91‑1 (MQ=255)

                                                                            aaTTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt               >  1:2625839/1‑91 (MQ=255)

                                                                                  aTCTTTGTGCATTGTTACAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACa         <  1:3552083/91‑1 (MQ=255)

                                                                                  aTCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACa         >  1:1090825/1‑91 (MQ=255)

                                                                                  aTCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACa         <  1:467783/91‑1 (MQ=255)

                                                                                  aTCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACa         <  1:1799134/91‑1 (MQ=255)

                                                                                   tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg        <  1:2309909/91‑1 (MQ=255)

                                                                                      ttGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTg     <  1:388744/91‑1 (MQ=255)

                                                                                        gtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCa   <  1:3076461/91‑1 (MQ=255)

                                                                                         tgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaa  >  1:4120255/1‑91 (MQ=255)

                                                                                         tgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaa  >  1:1483420/1‑91 (MQ=255)

CTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAA  >  USA300TCH1516_ALE/1744272‑1744451

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