BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,875	T→C	61.0%	intergenic (‑958/‑210)	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,875	0	T	C	61.0%	175.7 / 176.3	186	intergenic (‑958/‑210)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (36/37); new base C (59/53); total (96/90)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

AATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGG  >  CP000730/1929739‑1930015

aaTCTATCAGCTAGAATCACATAACTTATCAGAAGCAAGTGACTGTGCTTGTAACTTCGCTCTGCATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTAttt                                                                                                                                           <  1:365724/141‑1 (MQ=255)

                          tATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGtcatc                                                                                                                <  1:8121/141‑1 (MQ=255)

                               gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCCAATCATAGCACCAGTCATCAGTgg                                                                                                           >  2:660837/1‑141 (MQ=255)

                               atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                          >  2:288542/15‑141 (MQ=255)

                               atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                          >  2:463539/15‑141 (MQ=255)

                                tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                                         <  1:715230/128‑1 (MQ=255)

                                gcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                                         <  1:87935/128‑1 (MQ=255)

                                aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                          <  1:652195/141‑1 (MQ=255)

                                aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATATCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGtcagcctttga                                                                                                          >  2:19820/1‑133 (MQ=255)

                                 cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                        >  2:400108/13‑141 (MQ=255)

                                 aGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                                         >  2:568180/1‑141 (MQ=255)

                                  ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                       <  1:63586/130‑1 (MQ=255)

                                  ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                       <  1:62742/130‑1 (MQ=255)

                                  gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                        >  2:683064/1‑141 (MQ=255)

                                   gaagtagatatgCTTGTTTATTCACAGATGATATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGATGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                      <  2:725921/131‑1 (MQ=255)

                                   gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                      <  2:379441/131‑1 (MQ=255)

                                   gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                      <  2:562541/131‑1 (MQ=255)

                                    gagtcgatatgCTTGTATATTCACAGATGGTATCAGAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                     <  1:245784/132‑1 (MQ=255)

                                    aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                     >  1:434977/10‑141 (MQ=255)

                                     agtcgatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TTGTGAATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                    <  2:498252/133‑2 (MQ=255)

                                     agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGGGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCCCCAGTCATCAGTGGCTGTGct                                                                                                    >  2:74521/9‑140 (MQ=255)

                                       tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                                  <  1:288542/135‑4 (MQ=255)

                                       tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                                                                                                   >  2:756786/1‑141 (MQ=255)

                                        gACTGCGTTTGTATACTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTTGTCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                                  <  1:1077/141‑1 (MQ=255)

                                        gACCGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATGTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCCTCAGTGGCTGTGCCAt                                                                                                  <  2:702695/141‑1 (MQ=255)

                                         gatttgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                                <  1:48652/137‑6 (MQ=255)

                                         gatatgCTTTTAAATTCACAGATGGTATAATAATAGAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                                <  1:272892/137‑6 (MQ=255)

                                         gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                                >  1:257151/5‑136 (MQ=255)

                                         gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                                <  1:396443/137‑6 (MQ=255)

                                         gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACAAGTCATCAGTGGCTGTGctgttg                                                                                                >  1:11691/5‑136 (MQ=255)

                                         aatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                                <  1:727175/137‑6 (MQ=255)

                                           tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTGC‑GTTCCCATCATAGCACCAGTCATCTGTGGCTGTGCCAtttc                                                                                               >  2:331535/1‑139 (MQ=255)

                                           tgtgCTTGTATATTTGCCGATGATATAATAGTAAAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                               >  1:350511/1‑141 (MQ=255)

                                              gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                           <  1:704492/141‑3 (MQ=255)

                                               cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                          <  1:416547/141‑4 (MQ=255)

                                                ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                                          <  2:236979/141‑1 (MQ=255)

                                                 tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                        <  2:685505/141‑6 (MQ=255)

                                                  gTATATTCCCAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                       <  1:414022/141‑4 (MQ=255)

                                                   tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                      <  1:384396/141‑8 (MQ=255)

                                                       ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                  >  2:482053/7‑141 (MQ=255)

                                                       ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                  >  2:172752/7‑141 (MQ=255)

                                                       ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCTTAt                                                                                  >  2:116941/7‑141 (MQ=255)

                                                         cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACTAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                >  1:26139/5‑141 (MQ=255)

                                                         cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACTAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                >  1:46857/5‑141 (MQ=255)

                                                         cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                >  2:672549/5‑141 (MQ=255)

                                                         cacagatggtATAATAATAAAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGGGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTG‑GCTCACAGCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAttt                                                                                >  2:229468/5‑140 (MQ=255)

                                                           cagatggtATAATAATATAGTGGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                              <  2:183178/139‑1 (MQ=255)

                                                           cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                              >  1:692460/3‑141 (MQ=255)

                                                              atggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                           <  1:622819/141‑3 (MQ=255)

                                                               tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGGGGCGGTCCTGTTGCGTATTATACTTATTGGCTTAg                                                                          >  2:61134/4‑141 (MQ=255)

                                                               aggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                          <  2:690362/138‑4 (MQ=255)

                                                                ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGcatagt                                                                         >  1:586731/3‑136 (MQ=255)

                                                                 gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGTGAGTTGGGGCA‑TTA‑C‑TTTGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttatta                                                                        >  1:720412/2‑136 (MQ=255)

                                                                 gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttcgta                                                                        >  2:210700/2‑136 (MQ=255)

                                                                  tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGGGGAGGTGCATTCC‑T‑T‑TGCTA‑ATCATTTTCACATCGCCGGA‑AGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                       >  1:177596/1‑140 (MQ=255)

                                                                  tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTggtctagtaa                                                                       >  2:15738/1‑133 (MQ=255)

                                                                  tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                       <  1:309208/141‑2 (MQ=255)

                                                                  tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                       <  1:642279/141‑2 (MQ=255)

                                                                    taataaTGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                     <  2:251347/141‑4 (MQ=255)

                                                                    taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                     <  2:657295/141‑4 (MQ=255)

                                                                    taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                     >  2:300211/1‑141 (MQ=255)

                                                                      ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                   <  1:608995/137‑1 (MQ=255)

                                                                        aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                 <  1:702817/139‑1 (MQ=255)

                                                                        aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                 <  1:323271/139‑1 (MQ=255)

                                                                         gtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCaa                                                                 <  1:103884/141‑1 (MQ=255)

                                                                         atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                <  1:230485/140‑2 (MQ=255)

                                                                         atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTA‑ATAATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                >  2:195353/2‑140 (MQ=255)

                                                                         atagAGTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTACACTTATTGGCTTAGTAAACGCaac                                                                <  1:492174/140‑2 (MQ=255)

                                                                           agagTTCCCTATCGCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGCCGCAATa                                                               <  2:515252/141‑1 (MQ=255)

                                                                           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATGCTTATTGGCTTAGTAAACGCaacaa                                                              <  1:98392/141‑4 (MQ=255)

                                                                           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAAGCaacaa                                                              >  2:283390/1‑138 (MQ=255)

                                                                           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                              >  1:13729/1‑138 (MQ=255)

                                                                           agagTCGCCTATCTCTCAGGAGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                              >  1:466623/1‑138 (MQ=255)

                                                                            gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                             <  2:521906/141‑5 (MQ=255)

                                                                             agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                                             >  1:132538/1‑141 (MQ=255)

                                                                             agTCGCCTATCTCTCAAGCGTCAATTTAGACGCAGAGAGGAGGTGAA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAGTGCGTTTTTTTCTTAGTGGCCAAGTAGACGCAATACa                                                             >  1:193572/1‑141 (MQ=255)

                                                                              gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                            <  1:193340/141‑1 (MQ=255)

                                                                               tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                          >  1:653965/1‑134 (MQ=255)

                                                                               gcGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaa                                                           <  1:79128/140‑1 (MQ=255)

                                                                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                        <  1:317320/141‑10 (MQ=255)

                                                                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                        >  2:412835/1‑132 (MQ=255)

                                                                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                        >  2:596056/1‑132 (MQ=255)

                                                                                  ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                        >  2:263200/1‑141 (MQ=255)

                                                                                   cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTTAGTAGATGCAATACAAAATAg                                                       >  2:712681/1‑141 (MQ=255)

                                                                                       ctatctggggTCCATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGa                                                   <  2:269466/132‑1 (MQ=255)

                                                                                          tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                                <  2:36661/141‑1 (MQ=255)

                                                                                                  cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                                        <  1:668149/141‑1 (MQ=255)

                                                                                                   aaTTTAGACGCAGATAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGTCGCa                                       <  1:547667/141‑1 (MQ=255)

                                                                                                      ttAGACGCAGAGGGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                    <  1:146770/141‑1 (MQ=255)

                                                                                                      ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                    <  1:742979/141‑1 (MQ=255)

                                                                                                        aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACGTATAGCCGCACCAAt                                  <  1:568180/141‑1 (MQ=255)

                                                                                                        aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                  >  1:292570/1‑141 (MQ=255)

                                                                                                              agagagGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAAt                            >  1:503593/1‑141 (MQ=255)

                                                                                                              agagagGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACAGCATAGCAACAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAAt                            >  1:369635/1‑141 (MQ=255)

                                                                                                               gagagGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAACAGGTGACATATAGCCGCACCAATAAAAATc                           <  2:582040/141‑1 (MQ=255)

                                                                                                                      ggTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                    >  1:511171/1‑141 (MQ=255)

                                                                                                                      ggTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                    >  2:321042/1‑141 (MQ=255)

                                                                                                                              aa‑G‑G‑GGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa             <  1:233930/141‑1 (MQ=255)

                                                                                                                                           cTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg    >  2:255398/1‑141 (MQ=255)

                                                                                                                                             tATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg  <  1:756786/141‑1 (MQ=255)

AATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGG  >  CP000730/1929739‑1930015

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