BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,799,069	G→T	8.0%	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	8.0%	735.7 / 5.9	232	Q516K (CAA→AAA)	USA300HOU_RS14380	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (97/109); new base T (24/2); total (121/111)
Fisher's exact test for biased strand distribution p-value = 5.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

GTCTAATAAATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCTATA  >  USA300TCH1516_ALE/2798929‑2799208

gTCTAATAAATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTttg                                                                                                                                             >  1:391530/1‑141 (MQ=255)

  cTAATAAATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTttgat                                                                                                                                           >  4:45035/1‑141 (MQ=255)

     ataaATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTttgattga                                                                                                                                        <  3:682263/141‑1 (MQ=255)

         aTCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGtct                                                                                                                                    <  2:1381515/141‑1 (MQ=255)

          tCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGtctc                                                                                                                                   <  2:1634860/141‑1 (MQ=255)

            tGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCtt                                                                                                                                 >  3:338976/1‑141 (MQ=255)

            tGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCtt                                                                                                                                 <  3:261705/141‑1 (MQ=255)

             gATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCttt                                                                                                                                >  2:804745/1‑141 (MQ=255)

                   tGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTGTTTTTTTGAGTCTCttttttgaa                                                                                                                          >  2:1554957/1‑136 (MQ=255)

                      ttCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTc                                                                                                                       >  1:1584257/1‑141 (MQ=255)

                       tCATTAAATTTGCTAACTTATTTAAATCTGTGGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTCTTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGCGTCTCTTTTGTAAAttta                                                                                                                      >  1:366638/1‑139 (MQ=255)

                        cATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCaa                                                                                                                     <  2:550777/141‑1 (MQ=255)

                          ttAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATa                                                                                                                   >  1:1215001/1‑141 (MQ=255)

                          ttAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATa                                                                                                                   >  4:600027/1‑141 (MQ=255)

                           tAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATaa                                                                                                                  >  1:1099564/1‑141 (MQ=255)

                           tAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATaa                                                                                                                  >  4:61723/1‑141 (MQ=255)

                             aaTTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTCGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTACATCCAATAAat                                                                                                                >  1:1039515/1‑141 (MQ=255)

                               tttGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTTTTTTTGTCTCTTTTGTAAATTCAATAAatca                                                                                                              >  4:237135/1‑141 (MQ=255)

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

                                   cTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCg                                                                                                          >  1:178352/1‑141 (MQ=255)

                                    tAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGc                                                                                                         <  1:1438932/141‑1 (MQ=255)

                                    tAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGc                                                                                                         <  2:1520375/141‑1 (MQ=255)

                                     aaCTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCt                                                                                                        <  4:582434/141‑1 (MQ=255)

                                      aCTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTGTTTTTTTTAGTCTCTTTTGTAAATTCAATAAATCATCCGCtt                                                                                                       >  4:442054/1‑141 (MQ=255)

                                             tAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCAtt                                                                                                >  1:620041/1‑141 (MQ=255)

                                                 tCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGTTTTGTGCTCTTTTGTAAATTCAATTAATCATTCGCTTTACCttttatt                                                                                            >  2:1252889/1‑137 (MQ=255)

                                                 tCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAAt                                                                                            >  1:630247/1‑141 (MQ=255)

                                                  cTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATg                                                                                           <  1:521790/141‑1 (MQ=255)

                                                  cTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATg                                                                                           <  1:1330225/141‑1 (MQ=255)

                                                   tGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGc                                                                                          <  2:102879/141‑1 (MQ=255)

                                                       gTTTTCCCATTTTTATGCAATCGATCTAGTAAACTTGGACGATTTACTATTTGTGATAAAATATAATCTATATCTTTTGTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGAt                                                                                      <  2:73045/141‑1 (MQ=255)

                                                        ttttCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGAtt                                                                                     >  3:392406/1‑141 (MQ=255)

                                                            cccATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTTTGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAAc                                                                                 >  4:539842/1‑141 (MQ=255)

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

                                                             ccATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACt                                                                                <  1:1604635/141‑1 (MQ=255)

                                                              cATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTTAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTa                                                                               >  4:635391/1‑141 (MQ=255)

                                                               aTTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTaa                                                                              <  3:61723/141‑1 (MQ=255)

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

                                                                 ttttATTCAATCGATCTAGTAAACTTGGACGCTTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCCATTTTGTAAATTTACTAAATCATACGCTTTACCATTAAATGCAGCCTTAACTAAAt                                                                            >  1:73045/1‑141 (MQ=255)

                                                                  tttATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAAtt                                                                           <  2:451555/141‑1 (MQ=255)

                                                                     aTTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTaaa                                                                        >  3:139878/1‑141 (MQ=255)

                                                                       tCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTaaaaa                                                                      >  1:1501291/1‑141 (MQ=255)

                                                                        cAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAat                                                                     >  2:1457670/1‑141 (MQ=255)

                                                                          aTCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAattt                                                                   <  1:893005/141‑1 (MQ=255)

                                                                           tCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTTTTTTTTTTTTGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCGTTCAATGTCGATTTAACTACATTAAAAatttt                                                                  >  2:164796/1‑141 (MQ=255)

                                                                              aTCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAattttatt                                                               <  1:477793/141‑1 (MQ=255)

                                                                                cttgtcaaCTTGGACGATTTACTAGTGGTGATAAAATATAATCTATATCGTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCGATGCCGATTTAACTAAATTAGAAattttatttt                                                             <  2:1034672/135‑1 (MQ=255)

                                                                                cTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAattttatttt                                                             <  2:307163/141‑1 (MQ=255)

                                                                                cTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAattttatttt                                                             <  2:573582/141‑1 (MQ=255)

                                                                                  aGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATGCAATGCCGATTTAACTAAATTAAAAACTTTATTTTGa                                                           >  3:410034/1‑141 (MQ=255)

                                                                                   gTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGcc                                                                                         <  4:656518/110‑1 (MQ=255)

                                                                                     aaaCTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCt                                                        >  1:1278030/1‑141 (MQ=255)

                                                                                       aCTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCCTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGt                                                      >  1:1428103/1‑141 (MQ=255)

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

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

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

                                                                                          tGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTc                                                   >  1:917409/1‑141 (MQ=255)

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

                                                                                                  ttACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACAAAAATAAAAATTTTATTTTGATCTGTTTCTATTTTAg                                           <  2:395256/141‑1 (MQ=255)

                                                                                                   tACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGt                                          <  2:436196/141‑1 (MQ=255)

                                                                                                    aCTATCGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTa                                         <  4:705904/141‑1 (MQ=255)

                                                                                                       aTTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAAtt                                      <  3:64624/141‑1 (MQ=255)

                                                                                                        ttGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTa                                     <  1:1022480/141‑1 (MQ=255)

                                                                                                        ttGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTa                                     <  1:1215350/141‑1 (MQ=255)

                                                                                                        ttGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTa                                     <  3:834/141‑1 (MQ=255)

                                                                                                        ttGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTa                                     >  4:113415/1‑141 (MQ=255)

                                                                                                           gTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTATCTTGTAAATTCAATAAATCATCCGCTTTACCATTCAAGGTCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAAt                                  <  2:1473432/141‑1 (MQ=255)

                                                                                                             gATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCt                                >  1:1049837/1‑141 (MQ=255)

                                                                                                               tAAAATATAATCTATATCTTTTTTCTTTTTTTTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATATGATTTTATTTTAGTACTTAAATCTGc                              >  2:1327404/1‑141 (MQ=255)

                                                                                                                 aaaTATAATCTATATCTTTTTTCGTTTGATTGAGTCTCttt                                                                                                                                >  1:1471961/1‑41 (MQ=255)

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

                                                                                                                  aaTATAATCTATATCTTTTTTTTTTTTTTTTAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATTTTCCAg                           >  2:1116286/1‑141 (MQ=255)

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

                                                                                                                    tataATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATGATCCGCTTTACCATTCAATGCCGATTTAACTAACTTAAAAATTTTAATTTGATCTGTATCTATTTTAGTAATTAAATCTGCTAGTa                         >  2:1535922/1‑141 (MQ=255)

                                                                                                                    tataATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTa                         >  3:486993/1‑141 (MQ=255)

                                                                                                                      taATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAAt                       <  4:315765/141‑1 (MQ=255)

                                                                                                                       aaTCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAAtt                      <  3:271226/141‑1 (MQ=255)

                                                                                                                       aaTCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAAtt                      <  2:938684/141‑1 (MQ=255)

                                                                                                                        aTCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGtgacca                     >  3:534645/1‑136 (MQ=255)

                                                                                                                          ccagatCTTTTTTCGTTTTATTGAGTCTCTTTTGTACATTCAATAAATCATCCGCTTTACCATTCTTTGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTg                   <  4:306402/137‑1 (MQ=255)

                                                                                                                          cTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGA‑TTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGc                  <  2:201955/141‑1 (MQ=255)

                                                                                                                           tataTCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGc                  <  2:472487/141‑1 (MQ=255)

                                                                                                                             tatCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCt                >  3:110199/1‑141 (MQ=255)

                                                                                                                             tatCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCt                <  1:227459/141‑1 (MQ=255)

                                                                                                                             tatCTTTTGTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAAAAAAAAAAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCt                <  4:41372/141‑1 (MQ=255)

                                                                                                                                cTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCtttt             >  1:1297306/1‑141 (MQ=255)

                                                                                                                                 ttttttCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTg            <  3:406868/141‑1 (MQ=255)

                                                                                                                                  tttttCGTTTTTTTTAGTCTCTTTGGTAATTTCAATAAATCATCCGCTTTACCATTCGATGCTGATTTAACTAATTTAAAAATTTTATATTGATTTGTTTCTGTTTTAGTAATTAATTCTGCCAGTAATTTTGCCTTTTGt           >  2:490195/1‑141 (MQ=255)

                                                                                                                                  tttttCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGt           <  1:1516128/141‑1 (MQ=255)

                                                                                                                                   ttttCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTc          <  2:222635/141‑1 (MQ=255)

                                                                                                                                      tCGTTTGNTTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCttt       <  2:1239918/141‑1 (MQ=255)

                                                                                                                                      tCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCt                >  3:35986/1‑132 (MQ=255)

                                                                                                                                      tCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCt                <  4:35986/132‑1 (MQ=255)

                                                                                                                                       cGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAg                                           >  2:1157404/1‑104 (MQ=255)

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

                                                                                                                                        gTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCt     <  2:445236/141‑1 (MQ=255)

                                                                                                                                         tttGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCta    <  2:181388/141‑1 (MQ=255)

                                                                                                                                          ttgattgaGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCtat   >  1:423277/1‑141 (MQ=255)

                                                                                                                                           tgattgaGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCtata  >  1:263093/1‑141 (MQ=255)

                                                                                                                                           tgattgaGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCtata  >  1:1638283/1‑141 (MQ=255)

GTCTAATAAATCTGATCCTTGATTCATTAAATTTGCTAACTTATTTAAATCTGTCGTTTTCCCATTTTTATTCAATCGATCTAGTAAACTTGGACGATTTACTATTGGTGATAAAATATAATCTATATCTTTTTTCGTTTGATTGAGTCTCTTTTGTAAATTCAATAAATCATCCGCTTTACCATTCAATGCCGATTTAACTAAATTAAAAATTTTATTTTGATCTGTTTCTATTTTAGTAATTAAATCTGCCAGTAATTTTGCCTTTTGTCTTTCTATA  >  USA300TCH1516_ALE/2798929‑2799208

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