BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,799	C→A	55.3%	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	55.3%	‑0.0 / 149.9	118	intergenic (‑882/‑286)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (29/23); new base A (26/39); total (55/63)
Fisher's exact test for biased strand distribution p-value = 9.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

Only 100 of 121 total aligned reads displayed.

ATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGC  >  CP000730/1929663‑1929939

ataataAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgatg                                                                                                                                                <  1:904323/141‑1 (MQ=255)

    taaAATCGCCTACTTTCTTAAATTAACAATATCTATTCGCATAGAGTGTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCAGATAACTTATCCGAAGCAAGGGACTATGCGTGTATATTTGCCGATGatat                                                                                                                                            >  1:414113/1‑141 (MQ=25)

       aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATACCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                         <  1:182844/141‑1 (MQ=255)

       aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                         <  2:656158/141‑1 (MQ=255)

       aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                         >  1:394821/1‑141 (MQ=255)

          cGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                                                      <  1:147508/141‑1 (MQ=255)

           gCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                     <  2:832994/141‑1 (MQ=255)

               aCTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGt                                                                                                                                 <  2:469942/141‑1 (MQ=255)

                   tCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGcc                                                                                                                             >  2:101620/1‑141 (MQ=255)

                      tAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAATGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAt                                                                                                                          >  1:645611/1‑141 (MQ=255)

                          ttAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtctct                                                                                                                      >  2:781378/1‑141 (MQ=255)

                          ttAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtctct                                                                                                                      >  2:744222/1‑141 (MQ=255)

                                          tcATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGa                                                                                                      >  2:733059/1‑141 (MQ=255)

                                           cATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCCATTTAGAc                                                                                                     >  2:729231/1‑141 (MQ=255)

                                           cATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGAc                                                                                                     >  2:441428/1‑141 (MQ=255)

                                            aTAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGCTATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACg                                                                                                    >  1:22524/1‑141 (MQ=255)

                                              aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCa                                                                                                  >  2:883763/1‑141 (MQ=255)

                                               gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATCTTTGCCGACGATATAATCATAGAGTCGCCTATCTCTCAGGCGGCACTTTCAACGccg                                                                                                 >  2:550462/1‑139 (MQ=255)

                                               gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                                 >  1:520782/1‑141 (MQ=255)

                                                aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaga                                                                                                <  2:234262/141‑1 (MQ=255)

                                                 aTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagag                                                                                               <  2:804382/141‑1 (MQ=255)

                                                     gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga                                                                                           >  2:677448/1‑141 (MQ=255)

                                                               gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                               >  1:747109/1‑141 (MQ=255)

                                                                 gtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGGGGAGGGGTA‑TAA‑G‑G‑t                                                                           >  1:491417/1‑141 (MQ=255)

                                                                  tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTAACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑g                                                                         >  1:623881/1‑141 (MQ=255)

                                                                   aGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATATCTCAGGCGTCAATTTAGTCGCAGTGAGGAGGTGTA‑TAA‑G‑G‑T‑Ga                                                                        <  1:516640/141‑1 (MQ=255)

                                                                       gATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATgat                                                                    >  2:836125/1‑139 (MQ=255)

                                                                           cAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTCCTTGTATATTTGCCGATGATATAATAGTAGCGTCGCCCATCTCTCCGGCGTCCATTTAGACGCAGAGAGGAGGTGTA‑TAG‑G‑G‑A‑GATGC‑TTAtt                                                               >  2:609540/1‑141 (MQ=255)

                                                                              tCTATCAGCTAGAATCATATAACTTCACAGAAGAAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGGGTA‑TAA‑G‑G‑T‑GCTGC‑TTCTTTTc                                                            >  2:839342/1‑141 (MQ=255)

                                                                                 aTCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTTTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGtt                                                         <  1:392400/141‑1 (MQ=255)

                                                                                  tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTc                                                        <  1:945882/141‑1 (MQ=255)

                                                                                    ggCTCGAAGCATATAACTTATCAAAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGGCGCCTATCTCTCTTGCGTCAATTTCGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTcac                                                      <  2:93428/140‑1 (MQ=255)

                                                                                         gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCAcatcat                                                 >  1:136038/1‑141 (MQ=255)

                                                                                                       aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGtcatca                                   <  1:95442/141‑1 (MQ=255)

                                                                                                       aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGtcatca                                   >  1:500526/1‑141 (MQ=255)

                                                                                                           atcggaagtagatatgCTTGTATATTCACCGATGGTATGATAATATTGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                              <  1:405192/127‑1 (MQ=255)

                                                                                                           atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                              >  1:531183/15‑141 (MQ=255)

                                                                                                           atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                              <  1:385809/127‑1 (MQ=255)

                                                                                                           atcggaagtagatatgCTTGGATGTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                              <  1:588986/127‑1 (MQ=255)

                                                                                                            tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                             >  1:800637/14‑141 (MQ=255)

                                                                                                             tggaagtaggtatgCTTGTATATTCACAGATGGTATTATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCGCATCATAGCACCAGTCATCAGTGGCtg                            <  2:220925/129‑1 (MQ=255)

                                                                                                             cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                            >  1:809445/13‑141 (MQ=255)

                                                                                                             cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                            <  2:407656/129‑1 (MQ=255)

                                                                                                             cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCACTGGCtg                            <  2:521435/129‑1 (MQ=255)

                                                                                                             cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                            >  2:552933/13‑141 (MQ=255)

                                                                                                              ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                           <  1:532799/130‑1 (MQ=255)

                                                                                                               gaagtagatatgCTTGTATATTCACAGGTGGTACAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                          <  1:229031/131‑1 (MQ=255)

                                                                                                               gaagtagagtttcTTGGATATTCACAGTTGGTATAATCATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                          <  1:142484/129‑1 (MQ=255)

                                                                                                                aatgagagtagcgtGTATATTCACAGATGGTATAATAATAGGGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                         <  2:409890/128‑1 (MQ=255)

                                                                                                                aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                         <  1:132841/132‑1 (MQ=255)

                                                                                                                aagtagatatacTTGTATATTCACAGGTGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGCAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                         <  2:276159/130‑1 (MQ=255)

                                                                                                                aaggagagatgCTTGTATATTCACAGATTGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                         <  2:762431/132‑1 (MQ=255)

                                                                                                                 aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                         <  1:441428/141‑1 (MQ=255)

                                                                                                                 aGTGACTGTGCTTGTATATTTGCCGATGATATAACAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                         <  2:406851/141‑1 (MQ=255)

                                                                                                                 aGTGAATATGCTTGTATATTCATAGATGGTATAATAATAGTGTCGTCTATCTCTCAGGCGTCAATT‑‑GACGAAGTGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                        <  2:213878/141‑2 (MQ=255)

                                                                                                                  gtaggtatgCTTGTATATTCACAGATAGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑GCATTTTCGTTCACATCATAGCACTAGTCATCAGTGGCTGTGctg                       <  1:573061/134‑3 (MQ=255)

                                                                                                                  gtagatatgattGTATCTTCACCGATGGTATAGGAATAGAGTCGCTTATCGCTCAGGGGGCAATT‑‑GACGAAGAGAGGAGGGGCA‑TTA‑C‑T‑TGGCTAA‑TCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                       <  2:890004/131‑3 (MQ=255)

                                                                                                                  gtagatatgCTTGTATATTCACAGCTGGTAGAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑TGCGAATCATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                       <  2:591319/134‑3 (MQ=255)

                                                                                                                  gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                       <  1:583158/134‑3 (MQ=255)

                                                                                                                  gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGACTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCAGTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTCGCcgttcgt                       >  2:278690/8‑137 (MQ=255)

                                                                                                                  gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                        >  2:588576/1‑141 (MQ=255)

                                                                                                                   tagatctgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                      <  2:23786/135‑4 (MQ=255)

                                                                                                                   tGACTGTGCTTGTATATTTGCGGATGATATAATAGGAGAGTCGCCTATCTCGCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                       <  2:625450/141‑1 (MQ=255)

                                                                                                                    agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATCGCACCAGTCATCACTGGCTgcgctgtt                     >  1:799861/6‑137 (MQ=255)

                                                                                                                    agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                     >  1:745034/6‑137 (MQ=255)

                                                                                                                     gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                    <  1:650162/137‑6 (MQ=255)

                                                                                                                     gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                    >  2:651107/5‑136 (MQ=255)

                                                                                                                     gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                    >  1:430087/5‑136 (MQ=255)

                                                                                                                     gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                    >  1:128562/5‑136 (MQ=255)

                                                                                                                     aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTG‑TAA‑G‑G‑T‑GATGC‑TTATTTTCAGATACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                     >  1:489748/1‑141 (MQ=255)

                                                                                                                     aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                     <  2:529880/141‑1 (MQ=255)

                                                                                                                     aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                     >  1:185008/1‑141 (MQ=255)

                                                                                                                     aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                     >  2:723604/1‑141 (MQ=255)

                                                                                                                      atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                   <  1:510148/138‑7 (MQ=255)

                                                                                                                       tgtgCTTGTATATTTGCCGATGATATAATCGTAGCGTCGACTATCAAAAGGGAGTAAGATGATGCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                   <  2:198155/141‑1 (MQ=255)

                                                                                                                       tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                   <  1:755727/141‑1 (MQ=255)

                                                                                                                       tatgCTTGTATATTCACAGATGGTATAATAATCGAGTCGCCTACCTCTCAGGCGTCAATT‑‑GACGACGAGAGGAGTTGCA‑TTACT‑T‑T‑GCTAT‑GCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGctgttgcg                  <  1:689404/139‑8 (MQ=255)

                                                                                                                       tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg                  <  2:716259/139‑8 (MQ=255)

                                                                                                                       tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                  >  1:481777/3‑141 (MQ=255)

                                                                                                                       tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAAAT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                  >  1:568540/3‑141 (MQ=255)

                                                                                                                        ttgCTTGTATATTCACAGGTGTGATCATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                 <  1:352086/140‑1 (MQ=255)

                                                                                                                        atgCTTGTATATTCACAGATGGTATAATAGTAGAGTCACCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑G‑T‑T‑GCTAA‑TCATTTTGGTTCGCATCATAGAAACCACCATCAGTGGCTGTGctggtgcgt                 >  2:704367/2‑133 (MQ=255)

                                                                                                                        atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtttcgt                 >  2:141103/2‑133 (MQ=255)

                                                                                                                        atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                 >  2:936563/2‑141 (MQ=255)

                                                                                                                        atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                 >  2:469908/2‑141 (MQ=255)

                                                                                                                        atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                 >  1:829767/2‑141 (MQ=255)

                                                                                                                          gCTTGTCTATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat               <  2:801468/141‑3 (MQ=255)

                                                                                                                          gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat               <  1:97220/141‑3 (MQ=255)

                                                                                                                          gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat               <  1:271449/141‑3 (MQ=255)

                                                                                                                          gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat               <  2:589049/141‑3 (MQ=255)

                                                                                                                           cTTGGATATTCACAGATGGTAAAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt              <  1:803171/141‑1 (MQ=255)

                                                                                                                              gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata           <  1:658777/141‑7 (MQ=255)

                                                                                                                              gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata           <  2:144288/141‑7 (MQ=255)

                                                                                                                               tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac          <  1:439643/141‑8 (MQ=255)

                                                                                                                               tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac          <  1:758781/141‑8 (MQ=255)

                                                                                                                               tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac          <  2:436408/141‑8 (MQ=255)

                                                                                                                               tataTTCACAGATGGTATAATAACAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac          <  1:391189/141‑8 (MQ=255)

                                                                                                                                   tttGCCGATGATATAATACTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTa       <  2:922020/141‑1 (MQ=255)

                                                                                                                                        cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTTGc  >  2:872463/1‑141 (MQ=255)

                                                                                                                                        cGATGATATAATAGTAGAGTCGCCTATCACTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATATTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc  >  1:607712/1‑141 (MQ=255)

ATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGC  >  CP000730/1929663‑1929939

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