BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,799	C→A	52.4%	intergenic (‑882/‑286)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,799	0	C	A	52.4%	‑3.8 / 153.6	111	intergenic (‑882/‑286)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (27/25); new base A (31/27); total (58/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 113 total aligned reads displayed.

TGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT  >  CP000730/1929659‑1929936

tGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTACTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCCAGAATCATATACCTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGcc                                                                                                                                                <  2:123669/141‑1 (MQ=255)

tGTCACAATAAAATCGATGACTTTCTGAAATTAACAATATCTATTCTCATAGAATTTGTCCAAGTAAGTGTAGACGATTCAATTTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGcc                                                                                                                                                <  2:198251/141‑1 (MQ=25)

   cATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGCAGCTGAGT‑TCTTTTCTTGTATATTTGCCgat                                                                                                                                             <  2:1068457/141‑1 (MQ=21)

      aataaAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCgatgat                                                                                                                                          <  1:818537/141‑1 (MQ=255)

           aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                     >  1:932464/1‑141 (MQ=255)

           aaTCGCCTACTTTCTTAAAGTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                     >  1:719774/1‑141 (MQ=255)

              cGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                                                  <  2:1019512/141‑1 (MQ=255)

              cGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                                                  <  1:194269/141‑1 (MQ=255)

                   aCTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGGGCTTGTATATTTGCCGAGGATAGAATAGTAGAGt                                                                                                                             >  1:458261/1‑141 (MQ=255)

                     tttCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCg                                                                                                                           >  1:504767/1‑141 (MQ=255)

                         ttAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTa                                                                                                                      <  1:362728/141‑1 (MQ=255)

                          tAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTAt                                                                                                                     >  1:542137/1‑141 (MQ=255)

                             aTTAACAAGATCTATTCTCATAGAATGTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTAtctc                                                                                                                  <  2:353183/141‑1 (MQ=255)

                                 aCAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTGTCAg                                                                                                              >  2:739895/1‑141 (MQ=255)

                                     tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGTTTTTGTTTTTATTTTTGCCGATGATATAATAGTAGAGTCGC‑CTATGTCTCAGGCGt                                                                                                          >  1:857658/1‑141 (MQ=38)

                                     tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGt                                                                                                          >  2:1035769/1‑141 (MQ=255)

                                       tCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCa                                                                                                        >  1:371047/1‑141 (MQ=255)

                                          aTTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGACAAtt                                                                                                     >  2:390148/1‑141 (MQ=255)

                                                       ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                        <  1:245631/141‑1 (MQ=255)

                                                       ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTAAAAGATTAGATAATCGTTTA‑ACGCGCGAGCTCTCAGGCGTCAATTTAGACGCagagag                                                                                        <  1:1093049/141‑1 (MQ=255)

                                                        tGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAgagg                                                                                       >  1:997173/1‑141 (MQ=255)

                                                         gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAgaggg                                                                                      >  1:546631/1‑140 (MQ=255)

                                                           ccAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGaggagg                                                                                    <  1:335820/141‑1 (MQ=255)

                                                                  aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑Taa                                                                            <  2:719774/141‑1 (MQ=255)

                                                                   gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                          >  2:285229/1‑141 (MQ=255)

                                                                   gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCACAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGGGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                          >  2:603150/1‑141 (MQ=255)

                                                                              tCAATCTATCAGCTAGAATCATATAACTTATCAGCCGC‑CAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAt                                                             <  1:748055/141‑1 (MQ=255)

                                                                              tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAt                                                             <  2:138621/141‑1 (MQ=255)

                                                                                            aGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCAcatca                                               <  1:611307/141‑1 (MQ=255)

                                                                                              aaTCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATa                                             >  1:176627/1‑141 (MQ=255)

                                                                                              aaTCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATa                                             >  2:713111/1‑141 (MQ=255)

                                                                                                  atatAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCAc                                         <  1:616798/141‑1 (MQ=255)

                                                                                                   tataACTTATCACAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCAcc                                        <  1:477129/141‑1 (MQ=255)

                                                                                                    ataACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCTGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                       <  2:562550/141‑1 (MQ=255)

                                                                                                    ataACTTATCAGAAGC‑AAGGGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                       <  2:385224/141‑1 (MQ=255)

                                                                                                              aGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATCGTCGATGCGC‑CTATCTCTCAGGCGTCTATTTAGACGCAGAGAGGAGGGGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTGCACATCATAGCACAAGTCATGAGTg                             >  1:213763/1‑141 (MQ=255)

                                                                                                              aGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATAAGTg                             >  2:342894/1‑141 (MQ=255)

                                                                                                               gAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                            <  2:642341/141‑1 (MQ=255)

                                                                                                                atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                           >  1:298789/15‑141 (MQ=255)

                                                                                                                 tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          >  1:694400/14‑141 (MQ=255)

                                                                                                                 tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          >  1:779770/14‑141 (MQ=255)

                                                                                                                 tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          <  1:690224/128‑1 (MQ=255)

                                                                                                                 tcggaagtagatatgCCTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          >  1:654921/14‑141 (MQ=255)

                                                                                                                  cggaagtagatatgCTTGTATATTCACAGATGGTATGATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                         <  2:1075602/129‑1 (MQ=255)

                                                                                                                  cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGTGCCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                         >  1:483337/13‑141 (MQ=255)

                                                                                                                  cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGCGCTCCAGTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                         >  1:753601/13‑141 (MQ=255)

                                                                                                                   tgaaggtgatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:936672/130‑1 (MQ=255)

                                                                                                                   ggaagtagttatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:920351/130‑1 (MQ=255)

                                                                                                                   ggaagtagatatgCTTGTATATTTGCAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:825765/130‑1 (MQ=255)

                                                                                                                   ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCtgg                        >  1:755295/12‑140 (MQ=255)

                                                                                                                   ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:692563/130‑1 (MQ=255)

                                                                                                                   ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:289819/130‑1 (MQ=255)

                                                                                                                   ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        >  1:49264/12‑141 (MQ=255)

                                                                                                                   ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        >  1:755035/12‑141 (MQ=255)

                                                                                                                   c‑AAGTGATTGTGCTTGTATATTTGCAGATGATAGAATAGTAGAGTCGC‑CTACCTCTCAGGCGGCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:451705/141‑1 (MQ=255)

                                                                                                                    gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       >  2:116323/11‑141 (MQ=255)

                                                                                                                    gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       <  1:932688/131‑1 (MQ=255)

                                                                                                                    gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       >  1:341501/11‑141 (MQ=255)

                                                                                                                    gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       <  1:660359/131‑1 (MQ=255)

                                                                                                                     aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                      <  1:1076765/132‑1 (MQ=255)

                                                                                                                     aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                      <  1:175702/132‑1 (MQ=255)

                                                                                                                     aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                      <  1:328843/132‑1 (MQ=255)

                                                                                                                     aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGc                      >  1:561739/10‑141 (MQ=255)

                                                                                                                       gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGGCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                    >  1:173457/8‑139 (MQ=255)

                                                                                                                        tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GGCGAAGAGAGGAGGGGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGctgt                   >  1:563036/7‑138 (MQ=255)

                                                                                                                         agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCCCCAGTCACCAGTGGCGGTGcttgt                  >  2:168351/6‑137 (MQ=255)

                                                                                                                          ggtttgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 <  2:615632/137‑6 (MQ=255)

                                                                                                                          gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 >  1:649318/5‑136 (MQ=255)

                                                                                                                          gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 >  2:713399/5‑136 (MQ=255)

                                                                                                                          gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCggtactgttg                 >  2:1053726/5‑134 (MQ=255)

                                                                                                                          gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 >  2:687062/5‑136 (MQ=255)

                                                                                                                          gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGctgttg                 >  1:96340/5‑136 (MQ=255)

                                                                                                                          gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTAGTTCACATCATATCACCAGTAATCAGTGGCCGTGctgtgg                 >  2:989420/5‑136 (MQ=255)

                                                                                                                           atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                >  2:811502/4‑135 (MQ=255)

                                                                                                                           atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                >  2:889125/4‑135 (MQ=255)

                                                                                                                            tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGCGGCCAACTGC‑G‑T‑TTCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg               >  1:255225/3‑141 (MQ=255)

                                                                                                                            tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATAATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg               >  2:114621/3‑134 (MQ=255)

                                                                                                                             ttgCTTGGATATTTTCCTATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGCg               <  2:892052/140‑1 (MQ=255)

                                                                                                                             gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGT‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg               <  1:718776/141‑1 (MQ=255)

                                                                                                                             gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg               <  2:733135/141‑1 (MQ=255)

                                                                                                                             atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt              <  2:745872/140‑1 (MQ=255)

                                                                                                                             atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAAGCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGt              >  1:478695/2‑141 (MQ=255)

                                                                                                                             atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCCCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt              >  2:397605/2‑141 (MQ=255)

                                                                                                                               gcgtGTATATTAACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat            <  1:245101/138‑3 (MQ=255)

                                                                                                                               gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt             >  2:951174/1‑141 (MQ=255)

                                                                                                                               gCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat            <  1:112256/141‑3 (MQ=255)

                                                                                                                               gCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat            <  1:176831/141‑3 (MQ=255)

                                                                                                                                cTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt           >  2:159904/1‑141 (MQ=255)

                                                                                                                                cTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt           >  2:712100/1‑138 (MQ=255)

                                                                                                                                 ttGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta          <  2:455323/141‑5 (MQ=255)

                                                                                                                                 ttGTATATTCACAGACGGTATAATAATAGAGTCGC‑CTATCTCTCTGGCGTCAATGT‑‑ACGATGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta          <  2:759559/141‑5 (MQ=255)

                                                                                                                                  tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat         <  2:446697/141‑6 (MQ=255)

                                                                                                                                  tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat         >  2:84588/1‑139 (MQ=255)

                                                                                                                                  tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCTGTTGCGtattat         <  2:395724/141‑6 (MQ=255)

                                                                                                                                  tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat         <  2:761835/141‑6 (MQ=255)

                                                                                                                                   gTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt         <  2:932464/141‑1 (MQ=255)

                                                                                                                                    tataTTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt        <  2:504767/141‑1 (MQ=255)

                                                                                                                                     atatTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttttt       >  2:145239/1‑140 (MQ=255)

                                                                                                                                     atatTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc       >  2:1101970/1‑141 (MQ=255)

                                                                                                                                          tGCCGATGATATAATAGTAGAGTCGC‑ATATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt  >  2:687741/1‑141 (MQ=255)

TGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT  >  CP000730/1929659‑1929936

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