BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,843	Δ1 bp	70.5%	intergenic (‑926/‑242)	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,843	0	A	.	70.5%	435.7 / 308.4	237	intergenic (‑926/‑242)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (36/30); new base . (87/80); total (125/112)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

TCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAA  >  CP000730/1929705‑1929982

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

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

 cATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGAc                                                                                                                                                >  2:441428/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)

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

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

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

                                                             aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGtcatca                                                                              <  1:95442/141‑1 (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)

                                                                 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‑TCATTTTCGTTCACATCATAGCACCAGTCATCACTGGCtg                                                                       <  2:521435/129‑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)

                                                                       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)

                                                                        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)

                                                                        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)

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

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

                                                                           gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                               >  1:430087/5‑136 (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)

                                                                           aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTG‑TAA‑G‑G‑T‑GATGC‑TTATTTTCAGATACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                >  1:489748/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)

                                                                              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‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                            >  2:936563/2‑141 (MQ=255)

                                                                              atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                            >  1:567940/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:271449/141‑3 (MQ=255)

                                                                                      atattcaaagatggtATAATAATAGAGTCGCCTCTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                    >  2:344770/10‑136 (MQ=255)

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

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

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

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

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

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

                                                                                              cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTTGc                                             >  2:872463/1‑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)

                                                                                                atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                          >  1:794153/5‑139 (MQ=255)

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

                                                                                                 tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                         <  1:669884/138‑4 (MQ=255)

                                                                                                 tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                         <  1:106099/138‑4 (MQ=255)

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

                                                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                        >  1:612586/3‑141 (MQ=255)

                                                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                        >  1:204431/3‑141 (MQ=255)

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

                                                                                                   gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                       >  1:139950/2‑141 (MQ=255)

                                                                                                     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                     >  1:847028/1‑139 (MQ=255)

                                                                                                     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                     >  1:4748/1‑139 (MQ=255)

                                                                                                        ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                  <  1:38992/141‑1 (MQ=255)

                                                                                                        ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATAATACGGATTGGCTTAGTAAACGc                                  <  1:226877/141‑1 (MQ=255)

                                                                                                           atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCTGTTGCTTTGCTGTTGCGTATTATATTTATTGGCTTATTAAACGCaac                               >  2:147231/2‑140 (MQ=255)

                                                                                                            tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTACTTAGACGCAAt                               >  2:259973/1‑141 (MQ=255)

                                                                                                            tagAGCCGTCTATTTCTCAGTCGGCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTTATCAGCGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                              <  2:852145/141‑3 (MQ=255)

                                                                                                             agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                              <  1:913395/141‑1 (MQ=255)

                                                                                                             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGTACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                             <  2:932192/141‑4 (MQ=255)

                                                                                                             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                             >  2:682593/1‑138 (MQ=255)

                                                                                                              gagTCGCCTATCTTTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                            <  2:378383/141‑5 (MQ=255)

                                                                                                              gagTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                            <  1:733450/141‑5 (MQ=255)

                                                                                                               agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                           <  1:687557/141‑6 (MQ=255)

                                                                                                                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                        >  2:139467/1‑133 (MQ=255)

                                                                                                                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                        >  2:272426/1‑133 (MQ=255)

                                                                                                                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                       >  2:743584/1‑132 (MQ=255)

                                                                                                                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                       <  2:718271/141‑10 (MQ=255)

                                                                                                                   gCCTATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                       <  1:49383/141‑10 (MQ=255)

                                                                                                                    ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                       >  2:891788/1‑141 (MQ=255)

                                                                                                                    ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                      >  1:548759/1‑131 (MQ=255)

                                                                                                                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                     >  1:769679/1‑130 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCACTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                    >  1:937040/1‑141 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGTTTAGTAAACGCaacaaataattgata                   >  1:714068/1‑128 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatg                   >  1:699730/1‑128 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                   >  1:632144/1‑128 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAATGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                   >  1:542668/1‑128 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGCAGGCATTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                   >  1:448853/1‑128 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattcata                   >  2:327059/1‑128 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GATAA‑TCATTATCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                   >  2:649384/1‑128 (MQ=255)

                                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                   >  1:839580/1‑128 (MQ=255)

                                                                                                                        tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTGGACGCAATACAGAAtttgtg                   >  2:657552/1‑136 (MQ=255)

                                                                                                                             cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata              <  2:428250/141‑1 (MQ=255)

                                                                                                                                 cGTTAATTTAGACGCAGGGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc          <  1:309534/141‑1 (MQ=255)

                                                                                                                                     aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa      <  1:101620/141‑1 (MQ=255)

                                                                                                                                       tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑G‑GATGC‑TTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGAAATACACCATAGGAGACATATAGCCGCAcc    <  2:693638/141‑1 (MQ=255)

                                                                                                                                         tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCACTACAAAATAGGTGACATATAGCCGCCCCaa  >  1:377853/1‑141 (MQ=255)

TCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAA  >  CP000730/1929705‑1929982

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