BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,864	A→T	21.8%	intergenic (‑947/‑221)	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,864	0	A	T	21.8%	317.2 / 76.8	156	intergenic (‑947/‑221)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (50/72); new base T (20/14); total (70/86)
Fisher's exact test for biased strand distribution p-value = 7.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

Only 100 of 160 total aligned reads displayed.

AAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929724‑1930004

    gtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑t                                                                                                                                           >  1:658926/1‑141 (MQ=255)

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

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

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

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

                              aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTAGTTCACATCATAg                                                                                                                >  2:603952/1‑141 (MQ=255)

                                 atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCAc                                                                                                             >  2:327764/1‑141 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAGTT‑‑GCCGACGATAGACGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTAGTTCGCATCATAGCCCCCGTCAGCAGTGGCTgcgctgct                                                                                      >  2:219394/6‑137 (MQ=255)

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

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

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

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

                                                          tatgCTTGTATATTCACAGAAGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                   <  1:467640/139‑1 (MQ=255)

                                                           gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                   <  1:35556/141‑1 (MQ=255)

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

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

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

                                                            tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCACTTC‑A‑T‑T‑GCAAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta                                                                                 >  1:822430/1‑140 (MQ=255)

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

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

                                                                    tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                                         >  1:618788/9‑141 (MQ=255)

                                                                    tatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt                                                                          <  2:255223/141‑1 (MQ=255)

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

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

                                                                           aggtggtataATAAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                  <  2:514985/135‑1 (MQ=255)

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

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

                                                                            gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                 >  1:48936/1‑140 (MQ=255)

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

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

                                                                              gggtataATAATAGAGTCCCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                               <  1:711134/138‑4 (MQ=255)

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

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

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

                                                                               ggtacaataatAGATTCGCCTATCTCTCAGGCGTCCATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                              <  2:555807/136‑1 (MQ=255)

                                                                                gtataATAGGGGGGGCGCCTATCTCTCAGGCGTGAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                              <  2:191874/140‑1 (MQ=255)

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

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

                                                                                 tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                            >  2:115692/1‑140 (MQ=255)

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

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

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

                                                                                    aataatAGAGTGGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                         <  2:239390/141‑1 (MQ=255)

                                                                                     ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                        <  1:810716/141‑1 (MQ=255)

                                                                                     ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCGAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATTCTTATTGGCTTAGTAAACGc                                                        <  1:542625/141‑1 (MQ=255)

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

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

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

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

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

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

                                                                                           gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                  <  1:376056/141‑5 (MQ=255)

                                                                                           gagTCGCCGATGTGTCAGGCTGCAATCGAGGCGCAGAGGGGAGGTGTA‑TAA‑G‑G‑G‑GATGCTTATTTTCGTTCACATCATAGCTCCAGTCATTAGTGGCTGTGCCAGTGCGTTTTTGTCTGATGGGCTAAGTAGACGCAATAc                                                   <  2:171833/141‑1 (MQ=255)

                                                                                            agTCGCCTATCTGTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                 >  1:592453/1‑136 (MQ=255)

                                                                                            agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                 >  2:198806/1‑136 (MQ=255)

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

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

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

                                                                                              tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaa                                                >  2:314793/1‑141 (MQ=255)

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

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

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

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

                                                                                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                            >  1:311927/1‑131 (MQ=255)

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

                                                                                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCCTAGTAAACGCaccaaataatta                                            >  1:191735/1‑130 (MQ=255)

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

                                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                                          >  2:284760/1‑129 (MQ=255)

                                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCTTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                          >  1:781601/1‑129 (MQ=255)

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

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

                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTCATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                         >  2:131690/1‑128 (MQ=255)

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

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

                                                                                                      ctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑A‑G‑T‑GATGCTTATCTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATCCAAAATACGTgg                                        >  2:139871/1‑140 (MQ=255)

                                                                                                          cAGGCTTGAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                    <  2:663794/141‑1 (MQ=255)

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

                                                                                                            ggCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGAGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAtata                                  >  1:788125/1‑141 (MQ=255)

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

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

                                                                                                                      tggACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCaa                        <  1:423127/139‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                     ggTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc         <  1:347823/141‑1 (MQ=255)

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

AAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929724‑1930004

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