BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,799,069	G→T	9.2%	Q516K (CAA→AAA)	USA300HOU_RS14380 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,799,069	0	G	T	9.2%	312.3 / 2.1	111	Q516K (CAA→AAA)	USA300HOU_RS14380	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (40/58); new base T (11/2); total (51/60)
Fisher's exact test for biased strand distribution p-value = 5.69e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GTCTAATAAATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCTAT  >  USA300TCH1516_ALE/2798929‑2799207

gTCTAATAAATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTATTATTGGTGATAAAATATAATCTATTTCTTTTTTCGTttg                                                                                                                                            >  1:504171/1‑141 (MQ=255)

   taataaATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTCCTATTGGTTATAAAATATAATCTATATATTTTTTATTttgatt                                                                                                                                         >  1:119537/1‑141 (MQ=255)

   taataaATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTttgatt                                                                                                                                         <  2:115203/141‑1 (MQ=255)

     ataaATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTttgattga                                                                                                                                       <  2:1179665/141‑1 (MQ=255)

     ataaATCTGATCCCTGATTCATTAAATTTGTTAACTTATTTAAATCTGTCCTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTttgattga                                                                                                                                       <  2:781944/141‑1 (MQ=255)

        aaTCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGtc                                                                                                                                    <  1:966512/141‑1 (MQ=255)

          tCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTGCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGtctc                                                                                                                                  >  1:39517/1‑141 (MQ=255)

            tGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCtt                                                                                                                                <  2:917007/141‑1 (MQ=255)

               tCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTg                                                                                                                             <  1:657601/141‑1 (MQ=255)

               tCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTCTTCGTTTGATTGAGTCTCTTTTg                                                                                                                             <  2:855918/141‑1 (MQ=255)

               cccTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTGTATTCAATCGATCTAGTAAACTTGTACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTg                                                                                                                             <  1:1092763/140‑1 (MQ=255)

                ccTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATCAAATATAATCTATATCTTTTTTCGTTTGAGTGAGTCTCTTTTGt                                                                                                                            <  1:654717/141‑1 (MQ=255)

                ccTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGt                                                                                                                            >  1:382049/1‑141 (MQ=255)

                ccTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTATCTTTTGt                                                                                                                            >  2:171021/1‑141 (MQ=255)

                 cttttttCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTa                                                                                                                           <  2:393256/136‑1 (MQ=255)

                 cTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTa                                                                                                                           <  1:825928/141‑1 (MQ=255)

                   tGATTCATTAGATTTACTAACTTATTTAAATCAGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGCTAAAATATAATCAATATCTTTTTTTGTTTTATTGAGTCTCTTTTGTaaa                                                                                                                         >  2:853564/1‑141 (MQ=255)

                   tGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAACATATAATCTATATCTTTTTTCGTTTGATTGAGTTTGTTTTGTaat                                                                                                                         >  1:972673/1‑140 (MQ=255)

                   tGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTTTGTTTCTTTTGTaaa                                                                                                                         >  2:853414/1‑141 (MQ=255)

                   tGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTaat                                                                                                                         >  1:662475/1‑140 (MQ=255)

                   tGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTaaa                                                                                                                         <  1:189700/141‑1 (MQ=255)

                   tGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTaaa                                                                                                                         <  1:189713/141‑1 (MQ=255)

                     aTTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTTAATCTCTTTTTTTAAtt                                                                                                                       >  2:834226/1‑141 (MQ=255)

                     aTTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTGTTTTTTTTATTCTCTATTGTAAAtt                                                                                                                       >  2:1164479/1‑141 (MQ=255)

                     aTTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAAtt                                                                                                                       >  1:13666/1‑141 (MQ=255)

                     aTTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAAtt                                                                                                                       <  2:904522/141‑1 (MQ=255)

                     aTTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGAGTAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTTTAAAtt                                                                                                                       >  1:517699/1‑141 (MQ=255)

                           tAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATGGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATaa                                                                                                                 >  2:5875/1‑141 (MQ=255)

                            aaaTTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATaaa                                                                                                                <  1:196407/141‑1 (MQ=255)

                               tttGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAatca                                                                                                             <  1:599191/141‑1 (MQ=255)

                                   cTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTAAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTTTATTTTGTAAATTTAATAAattatcgg                                                                                                         >  2:117760/1‑139 (MQ=255)

                                     accTTATTTAAATCTGTCCTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCt                                                                                                       <  2:1034823/139‑1 (MQ=255)

                                      aCTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCtt                                                                                                      >  1:547097/1‑141 (MQ=255)

                                       cTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTGTGTCTCTTTTGTAATTTTATTAAATTATTCGCttt                                                                                                     >  2:535535/1‑141 (MQ=255)

                                        ttATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTa                                                                                                    >  2:1017183/1‑141 (MQ=255)

                                            ttAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCAt                                                                                                <  1:969544/141‑1 (MQ=255)

                                             tAAAGCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCATTTGTTTGTGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCAtt                                                                                               >  2:725872/1‑141 (MQ=255)

                                              acaTCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCAGCCGCTTTACCATTc                                                                                              <  1:45923/139‑1 (MQ=255)

                                              aaaTCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTc                                                                                              <  2:409990/141‑1 (MQ=255)

                                              aaaTCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTc                                                                                              <  1:760303/141‑1 (MQ=255)

                                                ctCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCaa                                                                                            <  1:145794/140‑1 (MQ=255)

                                                aTCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATCTCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCaa                                                                                            <  2:897438/141‑1 (MQ=255)

                                                    tttgTTTTCCCATCTTTATTCAATCGTTCTAGTAAACTTGGACGATTTAATCGTAGTGATACATTATAATCTATATCCTTTTTCGTTTGATTGGGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGcc                                                                                        <  1:1121862/138‑1 (MQ=255)

                                                    gTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATGTAATTTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGcc                                                                                        >  2:425293/1‑141 (MQ=255)

                                                      cGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGa                                                                                      <  1:403176/141‑1 (MQ=255)

                                                         tttCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGAttt                                                                                   >  1:1016392/1‑141 (MQ=255)

                                                          ttCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTa                                                                                  <  2:476927/141‑1 (MQ=255)

                                                          ttCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTa                                                                                  <  2:985295/141‑1 (MQ=255)

                                                           tCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTaa                                                                                 >  2:1043480/1‑141 (MQ=255)

                                                            cccATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAAc                                                                                <  2:223999/141‑1 (MQ=255)

                                                            cccATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAAc                                                                                <  2:981170/141‑1 (MQ=255)

                                                            cccATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATGTAAc                                                                                <  1:813567/141‑1 (MQ=255)

                                                                tttttATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTaaa                                                                            >  1:1038448/1‑141 (MQ=255)

                                                                      ttCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTaaaa                                                                      <  1:352169/141‑1 (MQ=255)

                                                                             gagcTAGTAAACTTGGACGATTTACTATTGATGATAAAACATAATCAATCTCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATATATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAattttat                                                               <  2:486004/138‑1 (MQ=255)

                                                                             gATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATATAACTAAATTAAAAattttat                                                               <  1:620687/141‑1 (MQ=255)

                                                                               tCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATCTGATCCATATCTTTTTTCGTTTGATTGAGGCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAACatgttcttt                                                             >  2:53679/1‑141 (MQ=255)

                                                                                cTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTAGTTGTAATTAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAattttatttt                                                            <  1:1177141/141‑1 (MQ=255)

                                                                                 tAGTAAACTTGGACGATTTACTATTGGTGATACATTGTGCTGTATTTCCTTTTTCGGTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATACGCTTTACCATTCAATGCCGATTTAACTAAATTAGAAATTTTATTTTg                                                           >  1:1182961/1‑141 (MQ=255)

                                                                                  aGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTGAATCTCTTTTGTAAATTCAATAAATCATCCGCTTTTCCATTCAATGCCGATTTGTCTTCTTTAAAATTTTTATTTTgt                                                          >  2:619889/1‑140 (MQ=255)

                                                                                  aGTAAAATTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTACAAATTTTATTCTGa                                                          >  2:235245/1‑141 (MQ=255)

                                                                                   gTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTGTTTTTGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGAt                                                         >  2:785578/1‑141 (MQ=255)

                                                                                   gTAAAATTGAACGATGTACGATAGGGGTTAAAATATAGTCGATATCTTTTGTCGTTTGATTGAGTCTCTTTTGTAAAGTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGAt                                                         >  2:880064/1‑141 (MQ=255)

                                                                                        cTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGtt                                                    <  1:788943/141‑1 (MQ=255)

                                                                                        cTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGtt                                                    <  2:39517/141‑1 (MQ=255)

                                                                                         ttGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGttt                                                   <  2:827340/141‑1 (MQ=255)

                                                                                           ggACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCt                                                 <  2:612407/141‑1 (MQ=255)

                                                                                           ggACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCt                                                 <  1:1017183/141‑1 (MQ=255)

                                                                                            gACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATCAAAACTCTTATTTTGATCTCTTTCtt                                                >  2:647006/1‑140 (MQ=255)

                                                                                            gACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCCCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTa                                                <  2:382049/141‑1 (MQ=255)

                                                                                                     cTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTTTTTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACTATTGAATGCCGATTTAATTAAATTATAAATTTTATTTTGTTTTGTTTATATTttattat                                       >  2:1098107/1‑140 (MQ=255)

                                                                                                     cTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTaa                                       <  1:5827/141‑1 (MQ=255)

                                                                                                        ttagTGATAAAATATAATCTATATCTTCTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTa                                    <  2:1176903/138‑1 (MQ=255)

                                                                                                             gATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCCATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCt                               <  2:982385/141‑1 (MQ=255)

                                                                                                                 aaaTATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCa                           >  2:724943/1‑141 (MQ=255)

                                                                                                                  aaTATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAg                          >  2:684672/1‑141 (MQ=255)

                                                                                                                   atatAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGt                         <  2:204031/141‑1 (MQ=255)

                                                                                                                   atatAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGt                         >  2:468052/1‑141 (MQ=255)

                                                                                                                   atatAATCTATATCTATTTTCGTTTGCTTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTCAACTAAATTAAAACTTTTATTTTGATCTGTTTCTATTCTAGTAATTAAATCTGCCAGt                         <  1:725872/141‑1 (MQ=255)

                                                                                                                    tataATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTa                        >  2:136981/1‑141 (MQ=255)

                                                                                                                    tataATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTATATTTTAGTAATTAAATCTGCCAGTa                        >  1:166302/1‑141 (MQ=255)

                                                                                                                     ataATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTaa                       <  1:1149143/141‑1 (MQ=255)

                                                                                                                       aaTCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCCGTTTCTATTTTAGTAATTAAATCTGCCAGTAAtt                     <  2:598290/141‑1 (MQ=255)

                                                                                                                            atatCTTTTTTCCTTTGATTGAGTCTCTTTTGTCAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTATATCTGCCAGTAATTTTGcc                <  2:670531/141‑1 (MQ=255)

                                                                                                                             tatTTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCt               <  2:119537/141‑1 (MQ=255)

                                                                                                                             tatCTTTTGTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCt               >  1:590673/1‑141 (MQ=255)

                                                                                                                              tttttttttCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGGAATTAAATCTGCNAGTAATGTTGCCtt              <  2:602333/138‑1 (MQ=255)

                                                                                                                               tCTTTTTTTTTTTTTTTTTGTCTCTTTTGTAAATTCCATAAATCATCCGCTTTACCATGCAATGCCGATTTACCTTATTTAAAAATTTTATTTTGTTCTGTTTTTTTTTTTGTTATTTATTCTACTAGTATTTTTGCGttt             >  2:1038363/1‑141 (MQ=255)

                                                                                                                               tCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCttt             >  1:106575/1‑141 (MQ=255)

                                                                                                                                cTTTTTTTGTTTGATTGATGCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCtttt            <  2:755056/141‑1 (MQ=255)

                                                                                                                                   ttttCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATGTTATTTTGCTCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGtt         >  1:992432/1‑140 (MQ=255)

                                                                                                                                   ttttCGTTTGATTGAGACTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCCGTTTCTATTTTAGTAGTTAAAGCTGCCAGTAATCTGGCCTTTTGtt         >  2:794169/1‑140 (MQ=255)

                                                                                                                                     ttCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCtt       <  1:65147/141‑1 (MQ=255)

                                                                                                                                      tCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCttt      >  1:843579/1‑141 (MQ=255)

                                                                                                                                      tCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCttt      <  2:547097/141‑1 (MQ=255)

                                                                                                                                       cctttGAGTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACTATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTACTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTc     <  2:950134/139‑1 (MQ=255)

                                                                                                                                       cGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTc     >  1:932766/1‑141 (MQ=255)

                                                                                                                                       cGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAAAATGATTTTGGTCTGTTTCTATTTTAGTAATTAAATTTGCCAGTATTTTTGACTTTTTTTTTTc     >  1:568323/1‑141 (MQ=255)

                                                                                                                                        gTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGAGCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCt    >  1:47880/1‑141 (MQ=255)

                                                                                                                                          ttgattgaGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATACAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCtat  >  1:230509/1‑141 (MQ=255)

GTCTAATAAATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCTAT  >  USA300TCH1516_ALE/2798929‑2799207

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