BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,929,872	2 bp→AA	intergenic (‑955/‑212)	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,872	0	G	A	83.7%	351.0 / 57.9	159	intergenic (‑955/‑213)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (11/15); new base A (70/63); total (81/78)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01
*	CP000730	1,929,873	0	C	A	84.1%	374.1 / 62.9	158	intergenic (‑956/‑212)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (10/15); new base A (69/63); total (80/78)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

Only 100 of 160 total aligned reads displayed.

ATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAG  >  CP000730/1929735‑1930013

||

aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCtt                                                                                                                                             <  1:501029/141‑1 (MQ=255)

   cAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGCTGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAGGattatt                                                                                                                                          >  1:31649/1‑141 (MQ=255)

      tCTATCAGCTAGAATCATATAACTCATCAGAAGCAAGTGCCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTc                                                                                                                                       <  1:156611/141‑1 (MQ=255)

               tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCAcatc                                                                                                                              <  1:2878/141‑1 (MQ=255)

                                   atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                         >  1:784066/15‑141 (MQ=255)

                                   atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                         >  1:308116/15‑141 (MQ=255)

                                    tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                                        >  1:8360/14‑141 (MQ=255)

                                    tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                                        >  1:223499/14‑141 (MQ=255)

                                    tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGAATTAC‑T‑T‑TGCTAATCATTTTCGTTTACATCATAGAACCAGTAATCAGTAGCt                                                                                                        >  1:611557/14‑141 (MQ=255)

                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGTGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                      <  1:127338/130‑1 (MQ=255)

                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                      <  1:689727/130‑1 (MQ=255)

                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                      <  1:752439/130‑1 (MQ=255)

                                       gaagtagatattcTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                                     <  1:858467/129‑1 (MQ=255)

                                        aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                    <  1:533139/132‑1 (MQ=255)

                                        aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                    <  1:234519/132‑1 (MQ=255)

                                         agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                   <  1:781268/133‑2 (MQ=255)

                                         agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                   <  1:196085/133‑2 (MQ=255)

                                         agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                   <  1:330582/133‑2 (MQ=255)

                                          gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTTATCAGTGGCTGTGctg                                                                                                  <  1:461429/134‑3 (MQ=255)

                                          gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                  <  1:319392/134‑3 (MQ=255)

                                          gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                  <  1:116348/134‑3 (MQ=255)

                                           tagatatgCTTGTATATTCACCGATGGTATAATAATCGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                                 <  1:4528/135‑4 (MQ=255)

                                           tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                                 <  1:868248/135‑4 (MQ=255)

                                           tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                                 <  1:789671/135‑4 (MQ=255)

                                           tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCTTCAGTGGCTGTGctgt                                                                                                 <  1:460774/135‑4 (MQ=255)

                                             gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                               <  1:747086/137‑6 (MQ=255)

                                             gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                               >  1:775013/5‑136 (MQ=255)

                                               tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGc                                                                                             <  1:342258/141‑1 (MQ=255)

                                               tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGc                                                                                             <  1:367868/141‑1 (MQ=255)

                                               tatgCTTGTATATTCACAGATGGTATGATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCGTCAGTGGCTGTGctgttgcg                                                                                             <  1:561547/139‑8 (MQ=255)

                                               tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCg                                                                                            <  1:259539/139‑1 (MQ=255)

                                               tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg                                                                                             <  1:453797/139‑8 (MQ=255)

                                                 tgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGt                                                                                           >  1:593053/1‑141 (MQ=255)

                                                 tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGta                                                                                          >  1:362093/1‑140 (MQ=255)

                                                   cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGtatt                                                                                        >  1:486822/1‑138 (MQ=255)

                                                     tgtatattcacagatggtATCATAATAGAGTCGCCTATCTCTCAAGCGTCAATT‑‑GACGTAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGCGCTGT‑TGCGtattat                                                                                      >  1:562686/13‑139 (MQ=255)

                                                     tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGtattat                                                                                      <  1:927463/141‑6 (MQ=255)

                                                     tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGtattat                                                                                      <  1:49769/141‑6 (MQ=255)

                                                      gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGtattata                                                                                     >  1:808288/12‑138 (MQ=255)

                                                      gTATATTTGCCGATGCTATAATAGGAGAGTCGCCTATCTCGCAGGCGGCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTGCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGtttttt                                                                                      <  1:28385/141‑1 (MQ=255)

                                                       tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGtattatac                                                                                    <  1:390988/141‑8 (MQ=255)

                                                           ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTAt                                                                                >  1:381859/7‑141 (MQ=255)

                                                            tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTAtt                                                                               >  1:312953/6‑141 (MQ=255)

                                                            tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGC‑TGCGTATTATACTTAtt                                                                               >  1:118684/6‑141 (MQ=255)

                                                             cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTg                                                                              >  1:340588/5‑141 (MQ=255)

                                                              acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTgg                                                                             <  1:272184/138‑1 (MQ=255)

                                                               cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGc                                                                            >  1:851500/3‑141 (MQ=255)

                                                                 gatgatATAAGAGTAGCGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCt                                                                           <  1:710638/141‑1 (MQ=255)

                                                                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCttag                                                                        <  1:652955/138‑4 (MQ=255)

                                                                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAg                                                                        >  1:890406/4‑141 (MQ=255)

                                                                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAg                                                                        >  1:434614/4‑141 (MQ=255)

                                                                   tggtataATAATAGAGTCGCCTATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCttag                                                                        <  1:20763/138‑4 (MQ=255)

                                                                    ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGt                                                                       <  1:9032/139‑1 (MQ=255)

                                                                     gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTAGTTAAGAACTTAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTa                                                                      >  1:674037/2‑141 (MQ=255)

                                                                     gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTa                                                                      >  1:522067/2‑141 (MQ=255)

                                                                      tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaa                                                                     >  1:745403/1‑140 (MQ=255)

                                                                       ataataataGAGTTGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaa                                                                    <  1:734922/141‑3 (MQ=255)

                                                                       ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaa                                                                    <  1:232504/141‑3 (MQ=255)

                                                                       ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaa                                                                    >  1:218836/1‑139 (MQ=255)

                                                                       ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaa                                                                    <  1:726146/141‑3 (MQ=255)

                                                                       ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaa                                                                    >  1:409396/1‑139 (MQ=255)

                                                                       agaataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaa                                                                    <  1:322338/139‑3 (MQ=255)

                                                                        taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaac                                                                   <  1:542973/141‑4 (MQ=255)

                                                                        taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaac                                                                   <  1:389641/141‑4 (MQ=255)

                                                                        taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaac                                                                   <  1:92418/141‑4 (MQ=255)

                                                                        taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaac                                                                   <  1:904626/141‑4 (MQ=255)

                                                                        aaataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTaaac                                                                   <  1:543048/140‑4 (MQ=255)

                                                                         aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATTCTTATTGGCTTAGTAAACg                                                                  <  1:509398/141‑1 (MQ=255)

                                                                         aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACg                                                                  <  1:212055/141‑1 (MQ=255)

                                                                          ataatagAGTCGCCTATCTCTCAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGc                                                                 <  1:568654/137‑1 (MQ=255)

                                                                               agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTGATCGGTGGCTGTGCTGT‑TGCGTAATATACTTATTGGCTTAGTAAACGGaacca                                                            >  1:102543/1‑138 (MQ=255)

                                                                               agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                            >  1:615824/1‑138 (MQ=255)

                                                                               agagTCGCCTATCTCTCAGCCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                            >  1:604522/1‑138 (MQ=255)

                                                                                gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATACGAGGATCT‑CGATGTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                           >  1:628115/1‑137 (MQ=255)

                                                                                 agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                          <  1:921946/141‑6 (MQ=255)

                                                                                 agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                          <  1:259167/141‑6 (MQ=255)

                                                                                 agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAA‑T‑T‑TGCTAATCCTTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                          >  1:183951/1‑136 (MQ=255)

                                                                                 agTCGCCTATCTCTCAGGCGACAATTTACACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGCTGCTTATTTTCGTTCACATCATAGAACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                                           >  1:599856/1‑141 (MQ=255)

                                                                                   tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                        >  1:306349/1‑134 (MQ=255)

                                                                                    cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                       >  1:796012/1‑133 (MQ=255)

                                                                                     gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                      >  1:908762/1‑132 (MQ=255)

                                                                                     gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                      >  1:750229/1‑132 (MQ=255)

                                                                                     gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                      <  1:849069/141‑10 (MQ=255)

                                                                                        tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATCCTTATTGGCTTAGTAAACGCaacaaataattaat                                                   >  1:70653/1‑129 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  >  1:245740/1‑128 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  <  1:667681/141‑14 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  >  1:606846/1‑128 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  >  1:767529/1‑128 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  >  1:740294/1‑128 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  >  1:747941/1‑128 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  >  1:854406/1‑128 (MQ=255)

                                                                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                  >  1:408851/1‑128 (MQ=255)

                                                                                          tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGT‑TGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                 >  1:794522/1‑127 (MQ=255)

                                                                                            tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCGTAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTCTTGGCTAAGGAGACGCAATACAAAATAgaggac                                                >  1:33226/1‑136 (MQ=255)

                                                                                                  gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAg                                          <  1:849353/141‑1 (MQ=255)

                                                                                                        aTTTAGACGCGGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAc                                    >  1:228735/1‑141 (MQ=255)

                                                                                                             gACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTCTCTTATTGGCTAAGTAGACGCAATACAAAAGAGGTGACATATAGCCGCACCAATa                               <  1:677712/141‑1 (MQ=255)

                                                                                                                         aggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                   >  1:104774/1‑141 (MQ=255)

                                                                                                                         aggTGTA‑TAA‑A‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                   <  1:250547/141‑1 (MQ=255)

                                                                                                                                              cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg  <  1:841347/141‑1 (MQ=255)

||

ATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAT‑TGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAG  >  CP000730/1929735‑1930013

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