BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,918	A→G	59.9%	intergenic (‑1001/‑167)	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,918	0	A	G	59.9%	124.6 / 170.9	173	intergenic (‑1001/‑167)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (40/29); new base G (60/43); total (100/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 174 total aligned reads displayed.

GTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTAT  >  CP000730/1929783‑1930058

ttgCTTGGATATTTTCCTATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                                                                                <  2:892052/140‑1 (MQ=255)

gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGTCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                                                                                <  1:718776/141‑1 (MQ=255)

atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                               <  2:745872/140‑1 (MQ=255)

atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAAGCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGt                                                                                                                                               >  1:478695/2‑141 (MQ=255)

  gcgtGTATATTAACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                                                                             <  1:245101/138‑3 (MQ=255)

  gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                                                                             <  1:112256/141‑3 (MQ=255)

  gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                                                                             <  1:176831/141‑3 (MQ=255)

   cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                                            >  2:712100/1‑138 (MQ=255)

     tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                                                                          >  2:84588/1‑139 (MQ=255)

     tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                                                                          <  2:446697/141‑6 (MQ=255)

     tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCTGTTGCGtattat                                                                                                                                          <  2:395724/141‑6 (MQ=255)

     tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                                                                          <  2:761835/141‑6 (MQ=255)

      gTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt                                                                                                                                          <  2:932464/141‑1 (MQ=255)

        atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttttt                                                                                                                                        >  2:145239/1‑140 (MQ=255)

        atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc                                                                                                                                        >  2:1101970/1‑141 (MQ=255)

         tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                                                                                      <  2:455949/133‑7 (MQ=255)

             cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                  >  1:724216/5‑141 (MQ=255)

             cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                  >  1:27522/5‑141 (MQ=255)

             cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAACTT‑‑ACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                  >  1:485173/5‑141 (MQ=255)

              acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                 <  1:65185/138‑1 (MQ=255)

              acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                 <  1:15864/138‑1 (MQ=255)

                agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                               >  1:452312/2‑141 (MQ=255)

                agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                               >  1:919763/2‑141 (MQ=255)

                 gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                              >  1:1074319/1‑140 (MQ=255)

                 gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                              >  1:123878/1‑140 (MQ=255)

                  atgatCTCAGAGGAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                                                                                                              <  2:1054650/141‑1 (MQ=255)

                   tgttataATAAGAGACGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAATAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                            <  1:222035/138‑4 (MQ=255)

                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                            >  1:334817/4‑141 (MQ=255)

                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTTTTGGCttat                                                                                                                            >  2:500914/4‑138 (MQ=255)

                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                            >  2:19058/4‑141 (MQ=255)

                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                            >  2:115399/4‑141 (MQ=255)

                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAACACCGATCAACAGTGGCTGTGCTGTTGCGTCTTATCCTTATTGGCTTAg                                                                                                                            >  2:827289/4‑141 (MQ=255)

                   tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                            <  1:734483/138‑4 (MQ=255)

                   tggtataAGAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                            <  1:965861/138‑4 (MQ=255)

                   tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTATACGAAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGATGTTTTCTTGtttggtca                                                                                                                             >  2:582386/1‑135 (MQ=255)

                    gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                                                            >  1:857951/1‑141 (MQ=255)

                     gtataATAATAGAGTCGGCTATCTCTCAGGGGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCGCCAGTCATCAGTGGCTGTGCTGTTGCGTGTTATACTTATTGGTTTAGTa                                                                                                                          >  2:256610/2‑141 (MQ=255)

                     gtataATAATAGAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGCTG‑CATTAC‑T‑T‑TGCTAATAATTTTCGTTCACATCATAGCACCAATCATCAGTGGCTGTGCTGTTGCGTATTATATGTATTGGCTTAGTa                                                                                                                          >  2:348131/2‑141 (MQ=255)

                     gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                          >  1:800885/2‑141 (MQ=255)

                     gtataATAATAGAGTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCCTACGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                          >  1:571971/2‑141 (MQ=255)

                      tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCTTATTGGCTTAGtca                                                                                                                         >  2:518844/1‑139 (MQ=255)

                      tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGACGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                         >  2:714008/1‑140 (MQ=255)

                        taataaTAGATTCTCCTATCTCCCGGGCGTCAATT‑‑GACGAAGAGAGGACGTGCCT‑TAC‑T‑T‑TGCTTCTCATTTTCGTTCACATCATAGCACCATTCATCAGTTGCTGTGCTGTTGCGTATTATTCTTATTGGCTTAGTAAAc                                                                                                                       <  2:115887/141‑1 (MQ=255)

                         aagaatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                      <  2:227011/136‑1 (MQ=255)

                          ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                     <  1:966583/141‑1 (MQ=255)

                          ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                     <  2:339298/141‑1 (MQ=255)

                           tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                     <  1:1071881/141‑1 (MQ=255)

                           tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                     <  1:1004251/141‑1 (MQ=255)

                            aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                     >  1:646883/3‑139 (MQ=255)

                             atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                  >  2:1066255/2‑140 (MQ=255)

                              tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                 <  1:571420/141‑3 (MQ=255)

                              tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                 <  1:280341/141‑3 (MQ=255)

                               agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                >  2:595827/1‑138 (MQ=255)

                               agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                <  2:298789/141‑4 (MQ=255)

                                 agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTCTTGGCTTAGTAAACGccacaaat                                                                                                              >  2:359396/1‑134 (MQ=255)

                                   tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAGGTCGCCGCAAT‑‑ACaa‑‑a                                                                                                         >  1:315166/1‑141 (MQ=255)

                                    cGCCTATCTCTCAGTCGTCAACT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCAGCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTGGTGAACGCaacaaataga                                                                                                           >  2:758655/1‑133 (MQ=255)

                                    cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aa                                                                                                        <  2:201777/141‑1 (MQ=255)

                                     gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑AA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aaa                                                                                                       >  2:551636/1‑140 (MQ=255)

                                     gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                          <  1:741885/141‑10 (MQ=255)

                                     gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                          <  2:526243/141‑10 (MQ=255)

                                      ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCATT‑‑ACAA‑‑AATa                                                                                                      >  1:610897/1‑141 (MQ=255)

                                      ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                         >  2:748463/1‑131 (MQ=255)

                                       cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAg                                                                                                     <  1:285229/141‑1 (MQ=255)

                                       cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAg                                                                                                     <  1:603150/141‑1 (MQ=255)

                                       cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                        >  2:318013/1‑130 (MQ=255)

                                       cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                        >  2:870350/1‑130 (MQ=255)

                                       cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataagtaa                                                                                                        >  2:742861/1‑130 (MQ=255)

                                        tATCTCTTAGCCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                       <  2:178232/141‑13 (MQ=255)

                                        tATCTCTCAGGCGTCAATT‑‑GACGAAGGGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                       <  2:510179/141‑13 (MQ=255)

                                         aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                      >  1:361686/1‑128 (MQ=255)

                                         aTCTCTCACGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGTATTATACTTAGTGGCTTAGTAAACGCaacaaataattaata                                                                                                      >  2:815114/1‑128 (MQ=255)

                                          tctctcAGGCTTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                     <  2:681465/141‑15 (MQ=255)

                                          tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTAAGTAGACGCAAT‑‑ACAC‑‑AATAGGta                                                                                                  >  2:897083/1‑140 (MQ=255)

                                           ctctcAGGCGTCAATATAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCAT‑‑ACAA‑‑AATAGGTGa                                                                                                 >  1:820861/1‑141 (MQ=255)

                                                    gTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGTTTATCTTCGTTCAAATCCTATCTTCCGTGATCAGTGTCTGTGCTATCGCGTTTTTTACTTATTGTCTTACTAACCGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑GGcc                                                                                      <  2:458261/141‑1 (MQ=255)

                                                     tCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATCGGTGAC‑ATAT‑AGCC‑g                                                                                    >  1:1002530/1‑141 (MQ=255)

                                                                     gaggagGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAACAAAccc                                                                    >  2:1071356/1‑141 (MQ=255)

                                                                      aggaggTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATTATTGGTTTTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCCATAAAAATcccc                                                                   >  1:526248/1‑141 (MQ=255)

                                                                       ggaggTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                  <  2:857951/141‑1 (MQ=255)

                                                                       ggaggTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATCGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                  <  2:729976/141‑1 (MQ=255)

                                                                         aggTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCa                                                                >  2:903794/1‑141 (MQ=255)

                                                                          ggTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                               <  1:1035769/141‑1 (MQ=255)

                                                                           cgataa‑taA‑G‑G‑AAATGCTGGTCTGCGTTCACATCAATGCGCCCCTCATCAGCGGCTGCGCCTTTTCGTTTTCTTCTTACCGGCAAACTATACGCAAT‑‑AGAA‑‑AATATGTGAC‑ATAT‑AGCT‑TCACCAATAAAAATCCCCTCAc                                                               <  2:857658/136‑1 (MQ=255)

                                                                             g‑tA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTTGCTAAGTAGACTCAAT‑‑ACAG‑‑AATAGGTGAC‑ATAT‑AGAC‑GAACCAATAAAAATCCCCTGAcgaa                                                            >  1:194561/1‑138 (MQ=255)

                                                                                  taA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGc                                                         >  2:632019/1‑141 (MQ=255)

                                                                                  taA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGAAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGc                                                         >  1:636235/1‑141 (MQ=255)

                                                                                    a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                       >  1:861554/1‑141 (MQ=255)

                                                                                    a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                       <  1:951174/141‑1 (MQ=255)

                                                                                           gATGCTGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAg                                                   <  1:390148/141‑1 (MQ=255)

                                                                                                ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                              >  2:988190/1‑141 (MQ=255)

                                                                                                 tactttCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTGTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGCGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                             <  2:279408/138‑1 (MQ=255)

                                                                                                  aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                            >  1:840814/1‑141 (MQ=255)

                                                                                                   tcttCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                           <  1:897083/139‑1 (MQ=255)

                                                                                                       cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                       >  2:778228/1‑141 (MQ=255)

                                                                                                       cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑AGAT‑AGCC‑GCACCAATAAAAATCCCCTCCCTACCGCAAAAAGTGAGGGGATTgg                                       >  1:21586/1‑141 (MQ=255)

                                                                                                                  tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCCCACTACTGC‑AATAGTGTGGGGATTGGTGcatt                                 >  1:216538/1‑140 (MQ=255)

                                                                                                                    gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAt                          <  1:342894/141‑1 (MQ=255)

                                                                                                                       ccAGTCATACATGGCTGTTCCATTTCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTGAAGACt                       >  2:1072782/1‑141 (MQ=255)

                                                                                                                                            aTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATtat  <  1:874452/141‑1 (MQ=255)

GTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTAT  >  CP000730/1929783‑1930058

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