BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,799	C→A	33.9%	intergenic (‑882/‑286)	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,799	0	C	A	33.9%	142.1 / 72.9	108	intergenic (‑882/‑286)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (37/34); new base A (14/23); total (51/57)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTA  >  CP000730/1929660‑1929934

gTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCg                                                                                                                                              <  1:33594/141‑1 (MQ=255)

gTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCg                                                                                                                                              >  2:742496/1‑141 (MQ=255)

  cATAATAAAATCGCCTACTTTCTTAAATTAATAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgat                                                                                                                                            <  2:409284/141‑1 (MQ=255)

  cATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgat                                                                                                                                            <  1:143628/141‑1 (MQ=255)

     aataaAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgagtgt                                                                                                                                         >  2:26577/1‑137 (MQ=255)

         aaaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataa                                                                                                                                     >  2:492273/1‑141 (MQ=255)

          aatatccTATTTTCTGTAAATAACATGATCTATTCTCAGACAACTTGTCCAATTAAGTGTAGACGATTCAATCTATCATCTCGAATCATATAACTTATCCGAAGCAAGTGACTGTGCTTGTATATTTGCCGATTATataat                                                                                                                                    <  2:519461/136‑1 (MQ=25)

            tCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtag                                                                                                                                  <  1:346892/141‑1 (MQ=255)

              tccTACTTTCTTAAATTAACAATATCTACTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                <  1:88840/140‑1 (MQ=255)

              cccTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAAGTAAGTGTAGACGATTCCATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagta                                                                                                                                <  1:763590/140‑1 (MQ=255)

                cTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTTTATTTGCCGATGATATAATAGTaga                                                                                                                              >  1:268565/1‑141 (MQ=255)

                      tCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGcc                                                                                                                        >  1:638656/1‑141 (MQ=255)

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                       <  2:764321/141‑1 (MQ=255)

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                       >  2:72687/1‑141 (MQ=255)

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                       <  2:602467/141‑1 (MQ=255)

                       cTTAAATTAACAATATCTATTCTCATAGAATTTGTCAAATTAAGTGTAGACGATTCAATCTATCAGCTAGATTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTTCCGATGATATAATAGTAGAGTCGACt                                                                                                                       >  2:333251/1‑141 (MQ=255)

                         tAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAt                                                                                                                     >  2:534176/1‑141 (MQ=255)

                               aaCAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCa                                                                                                               <  2:769535/141‑1 (MQ=255)

                                      tCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCa                                                                                                       <  1:695758/141‑1 (MQ=255)

                                         aTTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAAtt                                                                                                    >  1:379884/1‑141 (MQ=255)

                                            ctcATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAg                                                                                                 >  1:668581/1‑141 (MQ=255)

                                              cATAGAATTTGTCCAATTAAGTGTAGACGATTCAAACCATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGAc                                                                                               <  1:394381/141‑1 (MQ=255)

                                                 aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATCTAATAGGAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCa                                                                                            >  2:687706/1‑141 (MQ=255)

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

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

                                                  gAATTTGGCCAATTAAGTGGAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCag                                                                                           <  2:216814/141‑1 (MQ=255)

                                                       tGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAg‑agg                                                                                     <  1:602199/141‑1 (MQ=255)

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

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

                                                           cAATTAAGGGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                 <  1:247119/141‑1 (MQ=255)

                                                           cAATAAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgt                                                                                 <  2:25293/141‑1 (MQ=255)

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

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

                                                             aTTAAGGGGAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGgtgt                                                                               <  2:695065/141‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                         cGATTCAATCTATCGGCTAAAATCAAATAACTTATCAGAAGCAACTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGc                                                               >  1:423686/1‑141 (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)

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

                                                                               aaTCTATCAGCTAGAATCATATAACTTATCAGACGCAAGTGACTGTGCTTGTATATTTGCCGATGACATAATAGTAGAGTCGCCTATCTCTCAGGCG‑GCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAttt                                                         >  2:765749/1‑141 (MQ=255)

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

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

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

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

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

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

                                                                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                         >  1:307206/15‑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)

                                                                                                               tcggaagtagatatgCTTGTCTATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                        >  1:471002/14‑141 (MQ=255)

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

                                                                                                                cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCGATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                       >  1:456350/13‑141 (MQ=255)

                                                                                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                      <  1:619722/130‑1 (MQ=255)

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

                                                                                                                  gaagttggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑CGGAGTTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCtgtg                     <  1:255424/131‑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)

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

                                                                                                                     gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGAGGCGGGTTTTCGTTCACATCATGGCACCAGTCATCAGTGGCTGTGcc                   >  1:524609/1‑141 (MQ=255)

                                                                                                                      tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                 <  1:316784/135‑4 (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)

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

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

                                                                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg               >  1:646252/5‑136 (MQ=255)

                                                                                                                        ccggtgCTTGGAGATTTGCCGATGATATAATAGTAGAGTAGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGCGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCCTCAGTGGCTGTGCCAtt                <  2:209532/138‑1 (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)

                                                                                                                          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)

                                                                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGAGCTGTTGCg             <  1:567487/139‑1 (MQ=255)

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

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

                                                                                                                           gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGGGGCGGA‑CAA‑G‑C‑AGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATGGCg             >  1:177912/1‑141 (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)

                                                                                                                                gttattTTGGCCGATGATATAATAGTAGAGTCGCTTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑T‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt        <  2:395127/139‑1 (MQ=255)

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

                                                                                                                                 gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata      <  2:558452/141‑7 (MQ=255)

                                                                                                                                  tataTTCACAGATGGTATAATAATAGATTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACCAAGAG‑AGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac     <  1:365937/141‑8 (MQ=255)

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

                                                                                                                                      tttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCA‑TCAATTTAGACGCAGAG‑ACACGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCACCAGTGACTGTTCCCTTGCGtttttttactt  >  2:174548/1‑137 (MQ=255)

GTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTA  >  CP000730/1929660‑1929934

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