BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,870	A→C	69.8%	intergenic (‑953/‑215)	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,870	0	A	C	69.8%	250.2 / 90.2	157	intergenic (‑953/‑215)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (20/28); new base C (48/58); total (70/87)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

ATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGT  >  CP000730/1929740‑1930010

aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGGGGTGTA‑TAA‑G‑G‑T‑GA‑TTGTTGtttt                                                                                                                                    >  1:720683/1‑141 (MQ=255)

aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTCtttt                                                                                                                                    >  2:819802/1‑141 (MQ=255)

      cAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCGGAGAGGGGGTGTG‑TAA‑G‑G‑T‑GT‑TGCTTTTTTTCGTTca                                                                                                                              >  1:594507/1‑141 (MQ=255)

        gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTcaca                                                                                                                            >  2:35556/1‑141 (MQ=255)

          tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCAcatc                                                                                                                          <  1:808235/141‑1 (MQ=255)

            gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCAcatcat                                                                                                                        >  1:698946/1‑141 (MQ=255)

                  tataACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGAAcc                                                                                                                  >  2:791623/1‑141 (MQ=255)

                          aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                                                          <  1:138099/141‑1 (MQ=255)

                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                     >  2:192086/15‑141 (MQ=255)

                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                     >  2:247314/15‑141 (MQ=255)

                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCCCCAGCCATCAGTGCCt                                                                                                    >  1:42930/14‑141 (MQ=255)

                                cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                   >  2:136253/13‑141 (MQ=255)

                                 ggaagtagatatgCTTGTATTTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCAGTTGCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                  <  2:158240/130‑1 (MQ=255)

                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGGGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                  <  2:154340/130‑1 (MQ=255)

                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                  <  1:803573/130‑1 (MQ=255)

                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                  <  1:488992/130‑1 (MQ=255)

                                 gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                   <  1:167235/141‑1 (MQ=255)

                                  gaagtagataggCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                 <  1:817627/130‑1 (MQ=255)

                                   gagttgatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                <  2:321725/132‑1 (MQ=255)

                                   aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                <  2:558803/132‑1 (MQ=255)

                                   aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TGGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                >  1:337450/10‑141 (MQ=255)

                                   aaGTGACTGTGCTTGGATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                 <  1:59458/141‑1 (MQ=255)

                                    agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                               <  2:673894/133‑2 (MQ=255)

                                    agtagatatgCTTGTATATTCACAGATGGGATAATAATAGAGTCGCCTATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                               <  2:26005/133‑2 (MQ=255)

                                     gtagatttgCTTGTATATTCACAGATGGTATAATAATGGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                              <  2:11155/134‑3 (MQ=255)

                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                              >  1:349963/8‑139 (MQ=255)

                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                              >  1:760732/8‑139 (MQ=255)

                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTACGGCGACTAATT‑‑ACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                              <  1:823123/134‑3 (MQ=255)

                                      tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                             <  1:550427/135‑4 (MQ=255)

                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATAAGTGGCTGTGctgtt                                                                                            >  2:352146/6‑137 (MQ=255)

                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                           >  1:208871/5‑136 (MQ=255)

                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                           <  1:297781/137‑6 (MQ=255)

                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                           >  2:542625/5‑136 (MQ=255)

                                         cttttcTTGTATATTTGGCGATGATATAATAGTCGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                                                                                           <  2:784971/136‑1 (MQ=255)

                                           gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                         <  1:632163/141‑1 (MQ=255)

                                           atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                        <  2:738597/140‑1 (MQ=255)

                                           atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                        >  2:601971/2‑141 (MQ=255)

                                            tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta                                                                                       <  1:479830/141‑2 (MQ=255)

                                             gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                      <  1:550332/141‑3 (MQ=255)

                                             gCTTGGATATTGGCCGATGATAGAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                                                                                       <  2:199097/141‑1 (MQ=255)

                                               ttGTATATTCACAGATGGTGTAATAATAGAGTCGCCTATCTCTCAGGCGTCAACT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                    >  2:511428/1‑140 (MQ=255)

                                                 gTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt                                                                                   <  1:704293/141‑1 (MQ=255)

                                                  tataTTCACAGATGGTATAATAATAGATTCGCCTATCTTTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                 <  2:278859/141‑5 (MQ=255)

                                                    tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCAAATCATAGCACCAGTCATCAGTGCCTGTGCTGTTGCGtattatactt                                                                               >  2:298755/9‑135 (MQ=255)

                                                    tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                                               >  2:55213/9‑141 (MQ=255)

                                                    tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑TACTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                                               >  2:605023/9‑141 (MQ=255)

                                                      ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                             >  2:810716/7‑141 (MQ=255)

                                                       tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                            <  1:361057/136‑1 (MQ=255)

                                                        cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                           >  1:790995/5‑141 (MQ=255)

                                                          cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                         <  1:413275/139‑1 (MQ=255)

                                                           agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                        <  1:786131/140‑1 (MQ=255)

                                                           agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                        <  1:786025/140‑1 (MQ=255)

                                                            gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCATTCTTGGCtt                                                                       >  2:752379/1‑140 (MQ=255)

                                                            gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                                        <  1:619346/141‑1 (MQ=255)

                                                             atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                                                       <  1:490140/141‑1 (MQ=255)

                                                              tgttataAGAATAGAGTCGCCTAACTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTTCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                     <  2:303127/138‑1 (MQ=255)

                                                              tggtataATACTAGATGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                     <  2:756122/138‑4 (MQ=255)

                                                              gggtataAGAGTAGAGTCGCCTATCTTTCGTGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGTACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                     <  2:498392/138‑1 (MQ=255)

                                                               ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                    <  1:361531/139‑1 (MQ=255)

                                                               ggtataATAATAGAGTCGCCTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCTCATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                    <  1:241552/139‑1 (MQ=255)

                                                                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATGATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTGGCttggta                                                                   >  2:197150/2‑136 (MQ=255)

                                                                atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                    >  2:41776/1‑141 (MQ=255)

                                                                 tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                  <  2:233635/141‑2 (MQ=255)

                                                                 tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                  >  2:365784/1‑140 (MQ=255)

                                                                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                 >  1:156826/1‑139 (MQ=255)

                                                                  ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATGTTCGTTCGCAGCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTCAGTAg                                                                  >  2:505018/1‑141 (MQ=255)

                                                                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCAGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                >  2:191068/1‑141 (MQ=255)

                                                                    aataatAGAGTTGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                               >  1:307755/1‑141 (MQ=255)

                                                                      taatagAGTCGCCTATCTCTCAGGCCTCAATT‑‑GAAGAAGAGAGGTGGTGCATTAC‑T‑T‑T‑GT‑TAATCATTTTCCTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                             <  2:379340/138‑1 (MQ=255)

                                                                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                            >  2:650488/3‑141 (MQ=255)

                                                                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                            <  1:111478/139‑1 (MQ=255)

                                                                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCCTAGTAAACGCaa                                                            >  2:39152/3‑141 (MQ=255)

                                                                        atatagTCCCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                           <  1:107199/137‑2 (MQ=255)

                                                                          agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                          >  2:836573/1‑141 (MQ=255)

                                                                             gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GG‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                       >  2:325926/1‑141 (MQ=255)

                                                                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAaaaa                                                      >  1:13125/1‑141 (MQ=255)

                                                                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                      <  2:347035/141‑7 (MQ=255)

                                                                               cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                    >  1:580713/1‑133 (MQ=255)

                                                                                gCCTGTCTCTCAGGCGTCAATT‑‑GGCGAGGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                   <  2:133413/141‑10 (MQ=255)

                                                                                gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                   <  2:568671/141‑10 (MQ=255)

                                                                                gCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                   <  2:641506/141‑10 (MQ=255)

                                                                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                  <  1:814002/141‑11 (MQ=255)

                                                                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAGTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                 <  2:571334/141‑12 (MQ=255)

                                                                                  cTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                 <  1:775842/141‑12 (MQ=255)

                                                                                  atATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTACTAAACGCaccaaataattaa                                                 >  2:514129/2‑129 (MQ=255)

                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                <  2:42152/141‑13 (MQ=255)

                                                                                    ttctctcAGGCGTCAATTTAGACGCAGAGCGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                                <  1:230311/140‑1 (MQ=255)

                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                               >  2:198995/1‑128 (MQ=255)

                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                               >  1:417073/1‑128 (MQ=255)

                                                                                          cAGGCGTTGATTTAGACGAAGAATTTAGTTTTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                          <  2:625527/141‑1 (MQ=255)

                                                                                                tCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCg                                    <  2:122296/141‑1 (MQ=255)

                                                                                                     ttAGCCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                               <  1:331499/141‑1 (MQ=255)

                                                                                                     ttAGCCGCAGAGAGGAGGTGTA‑CAA‑G‑G‑T‑GC‑TGCTTATTTTCGAAAGAGTCATCGTTTATGGCCCCGATCGCTGTGCCATTGCGTTTTGTGCTAGTTTTCCAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                               <  2:804414/141‑1 (MQ=255)

                                                                                                           gCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                         <  2:515932/141‑1 (MQ=255)

                                                                                                           gCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                         <  2:576682/141‑1 (MQ=255)

                                                                                                                  ggaggTGTA‑TAA‑G‑G‑T‑GA‑TGCTGATTTTCGTTCGCATCATAGAACCAGTAATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                  <  2:610102/141‑1 (MQ=255)

                                                                                                                    aggTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTGTTTAGTAGACGCAATACAAAATGGGTGACATATACCCGCGCAAGTACAAATCCCCTcc                >  2:787829/1‑140 (MQ=255)

                                                                                                                     tgtgtA‑TAA‑C‑G‑T‑TA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc               <  1:52103/140‑1 (MQ=255)

                                                                                                                         ta‑taA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc           <  2:83127/141‑1 (MQ=255)

                                                                                                                             aa‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTCCCGCa        >  2:22551/1‑141 (MQ=255)

ATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGT  >  CP000730/1929740‑1930010

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