BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,848	C→A	63.3%	intergenic (‑931/‑237)	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,848	0	C	A	63.3%	136.3 / 192.8	183	intergenic (‑931/‑237)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (36/31); new base A (61/52); total (99/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 184 total aligned reads displayed.

AGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAA  >  CP000730/1929709‑1929988

aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCGAGAATCATATAACTTATCAGAGGCAAGTGACTGTGCTTGTATATTTGCCGATGAGATAATATTGGAGTCGCCTATCTCTAAGGCG‑TCAATTGAGCCGCa                                                                                                                                                   >  2:396959/1‑141 (MQ=255)

 gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCag                                                                                                                                                  <  2:692173/141‑1 (MQ=255)

      tGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCACATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAg‑agg                                                                                                                                            <  1:66082/141‑1 (MQ=255)

         ccAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑aggagg                                                                                                                                         >  1:539066/1‑141 (MQ=255)

          cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                                                                        <  2:175449/141‑1 (MQ=255)

          cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                                                                        <  2:417259/141‑1 (MQ=255)

           aaTTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgtg                                                                                                                                       <  2:156881/141‑1 (MQ=255)

                aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCGGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑CTAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑aa                                                                                                                                >  1:437386/1‑141 (MQ=255)

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

                 gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGCGTA‑T‑CA‑g                                                                                                                              >  2:100176/1‑141 (MQ=255)

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

                  tgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑g‑g                                                                                                                            <  1:601325/141‑1 (MQ=255)

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

                        cGATTCAATCTATCGGCTAAAATCAAATAACTTATCAGAAGCAACTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGc                                                                                                                     >  1:423686/1‑141 (MQ=255)

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

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

                              aaTCTTTCAGCTGAGATCATATCACTGATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTAttt                                                                                                               <  2:638542/141‑1 (MQ=255)

                                 cTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCg                                                                                                            <  1:534176/141‑1 (MQ=255)

                                  tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGt                                                                                                           >  1:32723/1‑141 (MQ=255)

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

                                          aGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCCcatca                                                                                                   >  2:455407/1‑141 (MQ=255)

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

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

                                                      cTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGtca                                                                                       >  2:643536/1‑141 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                    gtagatatgCCTGTATACTCACAGATGGTACAATAATAGCGTCGCCCACCTCGCGTCCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                        <  1:21780/134‑3 (MQ=255)

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

                                                                      agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                      >  1:550046/6‑137 (MQ=255)

                                                                       ccggtgCTTGGAGATTTGCCGATGATATAATAGTAGAGTAGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGCGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCCTCAGTGGCTGTGCCAtt                                                                      <  2:209532/138‑1 (MQ=255)

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

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

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

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

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

                                                                         tatgCTTGTATATGCACAGATGGTAAAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                   <  1:698596/139‑1 (MQ=255)

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

                                                                          gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGGGGCGGA‑C‑AA‑G‑C‑AGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATGGCg                                                                   >  1:177912/1‑141 (MQ=255)

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

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

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

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

                                                                                   tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                         >  1:1491/9‑141 (MQ=255)

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

                                                                                      tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                      >  1:736519/6‑141 (MQ=255)

                                                                                        acagatgggataataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                    <  1:71719/132‑1 (MQ=255)

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

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

                                                                                           gatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                 <  1:165398/141‑2 (MQ=255)

                                                                                             tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCC‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                               >  1:297551/4‑141 (MQ=255)

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

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

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

                                                                                                tataATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                             >  2:697637/1‑141 (MQ=255)

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

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

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

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

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

                                                                                                    ataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                        <  1:293933/141‑1 (MQ=255)

                                                                                                    ataGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                         <  1:247138/141‑1 (MQ=255)

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

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

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

                                                                                                       atagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                     >  2:565761/2‑140 (MQ=255)

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

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

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

                                                                                                          gagTCGCCTATCCCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                   <  2:41598/141‑1 (MQ=255)

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

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

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

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

                                                                                                             tCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                               >  2:113229/1‑134 (MQ=255)

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

                                                                                                               gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                             <  2:154858/141‑10 (MQ=255)

                                                                                                               gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                             <  2:231325/141‑10 (MQ=255)

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

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

                                                                                                                ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCC‑TCAC‑G‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataataa                            >  1:730734/1‑131 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                                        tcAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGAGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTAGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAACATGGGTGGCat                     >  2:256390/1‑141 (MQ=255)

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

                                                                                                                          aGGCG‑TCAATTTAGACGCAGAG‑AGGAGGGGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGGGGCatat                   >  2:587540/1‑141 (MQ=255)

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

                                                                                                                                            cAGAG‑CGGAACTGTG‑T‑AA‑G‑C‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa  <  1:756816/141‑1 (MQ=255)

AGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAA  >  CP000730/1929709‑1929988

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