BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,865:1	+C	10.0%	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	1	.	C	10.0%	650.3 / 77.8	170	intergenic (‑948/‑220)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (82/71); new base C (11/6); total (93/77)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-02

Only 100 of 179 total aligned reads displayed.

AGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAA  >  CP000730/1929725‑1930005

aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑Taa                                                                                                                                                 >  1:156165/1‑141 (MQ=255)

aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGAAA‑ACAAGTTTGACGGCGAG‑AGGAGGTGTa‑aaa                                                                                                                                                 >  1:292766/1‑141 (MQ=255)

 gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACCCAGAG‑CGGCGGTGTA‑AAA‑g                                                                                                                                               >  2:205334/1‑141 (MQ=255)

   gtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGGGTA‑TAA‑g‑gg                                                                                                                                            >  1:246895/1‑140 (MQ=255)

        cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGc                                                                                                                                      >  1:386619/1‑141 (MQ=255)

        cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCA‑TAAATGCAGACGCAGAG‑AGGAGGGGTA‑CAA‑G‑G‑TGATGc                                                                                                                                      >  1:478470/1‑141 (MQ=255)

                  tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGt                                                                                                                            >  2:624848/1‑141 (MQ=255)

                               cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGCTGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGTAGGTGTA‑TAA‑G‑G‑TGCTGCTTATTTTCGTTCACATCATAGCa                                                                                                               >  1:251606/1‑141 (MQ=255)

                                  ataACTTATCAGAAGCAAGAGACTGTGCTTGTATATTTGCCGATGATATAATAGCAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                                                                                            >  1:503233/1‑141 (MQ=255)

                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                <  1:117662/127‑1 (MQ=255)

                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                >  2:561683/15‑141 (MQ=255)

                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                >  1:307206/15‑141 (MQ=255)

                                              tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                               >  1:189291/14‑141 (MQ=255)

                                                ggaagtagatacgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                             <  2:505844/129‑1 (MQ=255)

                                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                            >  2:606465/11‑141 (MQ=255)

                                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                            >  2:606190/11‑141 (MQ=255)

                                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGCCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                            <  2:134865/131‑1 (MQ=255)

                                                  aattagatatgCCTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                           <  2:555817/132‑1 (MQ=255)

                                                  aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                           >  1:29408/10‑141 (MQ=255)

                                                  aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                            <  2:446337/141‑1 (MQ=255)

                                                   agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                          >  1:613140/9‑140 (MQ=255)

                                                   actagatacgCTTGTATATTCACAGACGGTATAATAATAGAGTTGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                          <  1:192797/132‑2 (MQ=255)

                                                   aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                           <  1:742496/141‑1 (MQ=255)

                                                     tagatatgCTTGTATATTCACAGATGGTATAAGAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                        <  2:464364/135‑4 (MQ=255)

                                                     atgatatgattGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                        <  2:156768/132‑4 (MQ=255)

                                                       gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                      >  1:24742/5‑136 (MQ=255)

                                                       aatattcTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                      <  1:706627/135‑6 (MQ=255)

                                                       aCTGTGCTTGTATATTTGCCTATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGGAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                                       <  2:315977/141‑1 (MQ=255)

                                                        atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgg                                                                                     >  1:534432/4‑135 (MQ=255)

                                                        atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                     <  1:315760/138‑7 (MQ=255)

                                                        atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                     <  1:684803/138‑7 (MQ=255)

                                                         tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                     <  1:60674/141‑1 (MQ=255)

                                                         tgtgCTTGTATATTCGCCGATGATATAATAGGAGAGTCGCCTATCTCTCCGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                     <  2:28513/141‑1 (MQ=255)

                                                         tatgCTTGTATACTCACAGAAGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                    <  1:567610/139‑1 (MQ=255)

                                                          gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                    >  2:402339/1‑141 (MQ=255)

                                                          gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCT‑AAAGAAGCAAAAGTGAGCCG‑AGTGAACGTCGAGTA‑TAA‑G‑G‑TGATGCTTGTTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTACCATTGCg                                                                                    >  1:370522/1‑141 (MQ=255)

                                                          atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                   <  1:105831/140‑1 (MQ=255)

                                                            gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                                                                                  >  1:691625/1‑141 (MQ=255)

                                                             cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                >  2:642691/1‑141 (MQ=255)

                                                                gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                             <  1:538346/141‑7 (MQ=255)

                                                                   tatTTACAGATGGTGTAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                          <  1:198709/141‑7 (MQ=255)

                                                                    attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                         >  2:170795/8‑140 (MQ=255)

                                                                    aattaccaatgatATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                          <  1:333251/133‑1 (MQ=255)

                                                                     ttcacagatgatATAAAAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                        <  2:23651/135‑1 (MQ=255)

                                                                      tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                       >  1:734401/6‑141 (MQ=255)

                                                                        acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTgg                                                                     <  2:585235/138‑1 (MQ=255)

                                                                         cagatggtATAATAATAGAGAAGCATAACTCTCTTACG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                    >  1:767767/3‑141 (MQ=255)

                                                                          agttggtataATAATAGAGTCGCCTAGCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                   <  1:639072/135‑1 (MQ=255)

                                                                          agatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                   >  2:211307/2‑141 (MQ=255)

                                                                             tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                >  1:430400/4‑141 (MQ=255)

                                                                             tggtataATAATAGAGGCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGCAGAG‑AGGAGGAGCATTAC‑T‑T‑TGCTAATCATTTTAGTTCACATAATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                >  1:646422/4‑141 (MQ=255)

                                                                             gggtataATAAGAGAGTCGCCTATCTCCCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCGAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                <  1:261843/138‑4 (MQ=255)

                                                                              ggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                               >  1:409270/3‑141 (MQ=255)

                                                                              ggtataAAAATAGAGTCGCCAATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCGCATCATAGCACCAGTCATCAGTGGCTGTGCTGTTACGTATTATACTTATTGGCTTAGt                                                               >  1:488767/3‑141 (MQ=255)

                                                                                tataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                             >  2:668643/1‑140 (MQ=255)

                                                                                 ataataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTACTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                            <  2:576585/141‑3 (MQ=255)

                                                                                 ataataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGATTAGtcaa                                                            >  2:465716/1‑138 (MQ=255)

                                                                                  aaataatAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                           <  2:386510/140‑4 (MQ=255)

                                                                                    ataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTATACGc                                                         <  1:525159/141‑1 (MQ=255)

                                                                                    agtatagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                         <  1:632896/137‑1 (MQ=255)

                                                                                     taatagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                        <  1:416749/138‑1 (MQ=255)

                                                                                      aatagATTCGCCTATCTCTCAGGCG‑TCAA‑‑TGGACTAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                       <  2:770001/139‑1 (MQ=255)

                                                                                      aatagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                       <  2:749139/139‑1 (MQ=255)

                                                                                        tagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaaa                                                     >  1:31180/1‑139 (MQ=255)

                                                                                        tagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                     >  2:758991/1‑139 (MQ=255)

                                                                                        tagAGTCGCCTATCCCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                     <  2:264399/141‑3 (MQ=255)

                                                                                        atgagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                     <  1:544764/139‑3 (MQ=255)

                                                                                         agagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTACACGCaacaa                                                    >  2:416835/1‑138 (MQ=255)

                                                                                          gagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                   >  1:59597/1‑137 (MQ=255)

                                                                                           agTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                  >  2:72282/1‑136 (MQ=255)

                                                                                           agTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                  >  2:480994/1‑136 (MQ=255)

                                                                                            gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                  <  1:72687/141‑1 (MQ=255)

                                                                                            gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                  <  2:32723/141‑1 (MQ=255)

                                                                                            gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                  <  2:698408/141‑1 (MQ=255)

                                                                                               gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTACACGCaacaaatgatt                                              >  1:244536/1‑132 (MQ=255)

                                                                                               gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                              >  2:24855/1‑132 (MQ=255)

                                                                                               gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                              <  2:155011/141‑10 (MQ=255)

                                                                                                ccTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                              >  2:193899/1‑141 (MQ=255)

                                                                                                ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                             >  1:491941/1‑131 (MQ=255)

                                                                                                ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                             >  1:457055/1‑131 (MQ=255)

                                                                                                 cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  2:85138/1‑130 (MQ=255)

                                                                                                 cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  2:101946/1‑130 (MQ=255)

                                                                                                 cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  1:744668/1‑130 (MQ=255)

                                                                                                 cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                            >  1:320741/1‑130 (MQ=255)

                                                                                                  tATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                           <  2:613140/141‑13 (MQ=255)

                                                                                                  tATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                                           >  1:694320/1‑129 (MQ=255)

                                                                                                   aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                          >  2:277348/1‑128 (MQ=255)

                                                                                                   aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                          >  1:325699/1‑128 (MQ=255)

                                                                                                   aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataagtaata                                          >  1:171442/1‑128 (MQ=255)

                                                                                                   aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                          >  2:666719/1‑128 (MQ=255)

                                                                                                    tctctcAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAATACAAAATAGGTg                                          >  2:205758/1‑141 (MQ=255)

                                                                                                        ccAGGCG‑GCAATTTAGACGCAGAG‑AGGAGGGGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAGTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                      <  1:342103/140‑1 (MQ=255)

                                                                                                                tCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATCGCCg                               >  2:644272/1‑141 (MQ=255)

                                                                                                                        gACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATa                       >  1:456067/1‑141 (MQ=255)

                                                                                                                           gctgag‑cgGAGGTGTA‑TAA‑G‑G‑GGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                    <  2:379884/138‑1 (MQ=255)

                                                                                                                                g‑aggagGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc               <  1:624848/141‑1 (MQ=255)

                                                                                                                                  aggaggTGTG‑GAA‑G‑G‑TGATGCTTAGTTTCGTTCACATCATAGCACCCGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATCCAAAATAGGTGACATATAGCCGCACCAATAAAAATcctc              >  1:566864/1‑139 (MQ=255)

                                                                                                                                  aggaggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAGTAAACATcccc              >  1:588547/1‑141 (MQ=255)

                                                                                                                                    gaggTGTA‑TAA‑G‑G‑TGATGATAGGGTGCTGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc            >  2:748512/1‑141 (MQ=255)

                                                                                                                                            ttaA‑G‑G‑TCGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg     <  2:539066/140‑1 (MQ=255)

AGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAA  >  CP000730/1929725‑1930005

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