BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,804	A→G	56.6%	intergenic (‑887/‑281)	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,804	0	A	G	56.6%	65.0 / 123.6	119	intergenic (‑887/‑281)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (29/22); new base G (26/40); total (56/63)
Fisher's exact test for biased strand distribution p-value = 6.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 147 total aligned reads displayed.

TAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAG  >  CP000730/1929667‑1929943

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

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

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

           tcTTTCTTAAATTAATAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGAGTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGt                                                                                                                                     <  2:469758/140‑1 (MQ=255)

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

                      ttAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtctct                                                                                                                          >  2:781378/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)

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

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

                                             aTTTGTGCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagag                                                                                                   <  2:34544/141‑1 (MQ=255)

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

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

                                                 gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga                                                                                               >  2:829688/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‑TCATTTTCGTTCACATCATAGGACCAGTCATCAGTGGc                                  >  1:470594/15‑141 (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)

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

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

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

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

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

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

                                                                                                         cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                <  1:273132/129‑1 (MQ=255)

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

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

                                                                                                           gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCtgtg                              <  1:615758/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‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                             <  2:111537/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)

                                                                                                             aGTGAATATGCTTGTATATTCATAGATGGTATAATAATAGTGTCGTCTATCTCTCAGGCGTCAATT‑‑GACGAAGTGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                            <  2:213878/141‑2 (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)

                                                                                                               tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                          <  2:743504/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)

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

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

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

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

                                                                                                                 aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                         >  1:185008/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‑‑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‑TCATTTTCGTTCACATCATAGCACCAGTCATTAGTGGCGGTGCTGTTGCGt                     >  2:623953/2‑141 (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)

                                                                                                                    atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                     >  1:567940/2‑141 (MQ=255)

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

                                                                                                                       cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                  >  2:887091/1‑141 (MQ=255)

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

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

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

                                                                                                                             tcttcacagatggtATAATAATAGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt            <  1:867040/133‑1 (MQ=255)

                                                                                                                               ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAAGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt          >  1:604855/7‑141 (MQ=255)

                                                                                                                                tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt         >  2:728504/6‑141 (MQ=255)

                                                                                                                                tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt         >  2:193956/6‑141 (MQ=255)

                                                                                                                                  acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg       <  1:105422/138‑1 (MQ=255)

                                                                                                                                   cagatggtATAATACTAGAGTCGCCTATCTCTCGGGCCTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc      <  1:113392/139‑1 (MQ=255)

                                                                                                                                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc      <  1:296647/139‑1 (MQ=255)

                                                                                                                                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc      >  1:269137/3‑141 (MQ=255)

                                                                                                                                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc      <  2:353140/139‑1 (MQ=255)

                                                                                                                                    cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc      <  1:777770/141‑1 (MQ=255)

                                                                                                                                    agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt     >  1:839811/2‑141 (MQ=255)

                                                                                                                                    agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt     <  2:872901/140‑1 (MQ=255)

                                                                                                                                     gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt    >  1:234033/1‑140 (MQ=255)

                                                                                                                                     gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt    >  1:393125/1‑140 (MQ=255)

                                                                                                                                     gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt    <  2:872467/141‑2 (MQ=255)

                                                                                                                                      atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa    <  2:69720/141‑1 (MQ=255)

                                                                                                                                       tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTCCGTTTTTTTCTTATTGGCCaa   >  1:169448/1‑141 (MQ=255)

                                                                                                                                        gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg  <  1:649515/141‑1 (MQ=255)

                                                                                                                                        gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg  >  2:494724/1‑141 (MQ=255)

TAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAG  >  CP000730/1929667‑1929943

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