BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,865	A→C	43.5%	intergenic (‑948/‑220)	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,865	0	A	C	43.5%	29.0 / 255.3	184	intergenic (‑948/‑220)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (61/43); new base C (36/44); total (97/87)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

Only 100 of 187 total aligned reads displayed.

CGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929733‑1930004

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

        tCTATCAGCTAGAATCATATAATTTAGCAGAAGCAAGTGACTGTGCTTTTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTc                                                                                                                              <  2:171522/141‑1 (MQ=255)

               gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTcaca                                                                                                                      >  2:242692/1‑141 (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:463539/15‑141 (MQ=255)

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

                                      tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                              <  1:715230/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:62742/130‑1 (MQ=255)

                                        ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                            <  1:63586/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:562541/131‑1 (MQ=255)

                                         gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                           <  2:379441/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)

                                           aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                          >  2:233191/1‑141 (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)

                                                 tggccTTGTATATTTGCCCATGATATAATAGTAGTGTCGCCTATCGCTCAGGCGTCAATTTAGACGTAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                    <  2:206841/137‑1 (MQ=255)

                                                   tgCTTGTATATTTGGCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGt                                                                                  <  2:286128/141‑1 (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)

                                                       tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                             >  1:762306/1‑139 (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)

                                                           tatgcggcGATGATATAATAGTAGAGCCGCCTATCTCTCAGGTGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt                                                                          <  1:94473/134‑1 (MQ=255)

                                                             ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                       >  2:482053/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:254977/5‑141 (MQ=255)

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

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

                                                                  agatggtATAATAATCGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGTAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACt                                                                  >  2:214743/2‑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)

                                                                        tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGGGGAGGTGCATTCC‑T‑T‑TGCTA‑ATCATTTTCACATCGCCGGA‑AGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                            >  1:177596/1‑140 (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)

                                                                          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)

                                                                           aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                         >  1:81441/1‑141 (MQ=255)

                                                                           aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                         >  2:173057/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)

                                                                               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‑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)

                                                                                      cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                              >  2:149529/1‑133 (MQ=255)

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

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

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

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

                                                                                        ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTA‑ATCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                            >  2:184971/1‑131 (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)

                                                                                                 cAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTA‑CTCATTTTA‑GTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCTGTTCCGTATTATACTTATTGGCTTAGTAAACGCAACACATAAttaataggcgac                                   >  2:111110/1‑130 (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)

CGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATG‑CTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929733‑1930004

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