BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,872	G→A	70.0%	intergenic (‑955/‑213)	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	70.0%	238.9 / 155.0	196	intergenic (‑955/‑213)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (34/25); new base A (67/70); total (101/95)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 196 total aligned reads displayed.

TCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAT  >  CP000730/1929741‑1930007

tCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGGG‑TA‑TAA‑G‑G‑TGATGCTTATTTTc                                                                                                                                >  1:89949/1‑141 (MQ=255)

         tAGAATCATATAACCTATCAGAAGCATCTGATTGGGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTACTTTCGTTCAcatc                                                                                                                       <  2:510327/141‑1 (MQ=255)

             aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAg                                                                                                                   >  1:120961/1‑141 (MQ=255)

                          tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                                                                                                      >  2:702850/1‑141 (MQ=255)

                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                  >  2:216025/15‑141 (MQ=255)

                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                  <  1:243169/127‑1 (MQ=255)

                                ggaagtggatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                               <  1:303739/130‑1 (MQ=255)

                                ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAAGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCagt                                                                                               >  2:295305/12‑141 (MQ=255)

                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                              >  2:446819/11‑141 (MQ=255)

                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                              <  1:260132/131‑1 (MQ=255)

                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                              <  1:79981/131‑1 (MQ=255)

                                 gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                              <  2:202120/131‑1 (MQ=255)

                                 cAAGTGACTGTGCTTGTATTGTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                               <  2:712038/141‑1 (MQ=255)

                                  aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                             <  1:194823/132‑1 (MQ=255)

                                  aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                             <  2:707931/132‑1 (MQ=255)

                                   agttgatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                            <  2:340707/133‑2 (MQ=255)

                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                           <  2:261174/134‑3 (MQ=255)

                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAAGCATCAGTGGCTGTGctt                                                                                           >  1:648556/8‑139 (MQ=255)

                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                           >  1:665767/8‑139 (MQ=255)

                                    gtagatatgCTTGTACATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                           <  1:378046/134‑3 (MQ=255)

                                     tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                          <  1:22001/135‑4 (MQ=255)

                                       gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGgtgttctgttt                                                                                        >  1:297441/5‑134 (MQ=255)

                                        ttatgCTTGTATATTCACAGTTGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                       <  2:119996/138‑7 (MQ=255)

                                        agatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                       <  2:612720/138‑7 (MQ=255)

                                         tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGGGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                      >  2:111990/3‑141 (MQ=255)

                                         tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                      >  1:564133/3‑141 (MQ=255)

                                         gatgcgggTATATTCCCAGATGGTATAAAAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTGCACATCATAGCACCAGCCATCAGTGGCTGTGctgttgcg                                                                                      <  2:72792/134‑8 (MQ=255)

                                          ttgCTTGTATATTCACAGATGGTATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATGGGACCATTCATCAGTGGCTGTGCTGTTGCGt                                                                                     <  2:398364/140‑1 (MQ=255)

                                          atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGt                                                                                     >  2:124454/2‑141 (MQ=255)

                                          atgCTTGTATATTCAAAGATGGTATAATAATAGAGACGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGt                                                                                     <  2:479971/140‑1 (MQ=255)

                                              ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 <  2:631932/141‑5 (MQ=255)

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

                                                gTATATTTGCCGATGATATACTAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGCCGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTGTCGTTCACATCATAGCACCGGTCATCCGTGGCTGTGCCATTTCGattttt                                                                                >  2:242328/1‑141 (MQ=255)

                                                 tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                              <  2:58166/141‑8 (MQ=255)

                                                 tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                              <  1:333845/141‑5 (MQ=255)

                                                 tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                              <  1:474462/141‑8 (MQ=255)

                                                  atattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                             >  1:69526/10‑136 (MQ=255)

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

                                                     ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCGT‑TAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTCTTATACTTAt                                                                          >  1:633238/7‑141 (MQ=255)

                                                     ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                          >  2:77152/7‑141 (MQ=255)

                                                     ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                          >  2:589971/7‑141 (MQ=255)

                                                      ttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAt                                                                          >  1:671676/1‑141 (MQ=255)

                                                      tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                         >  2:780498/6‑141 (MQ=255)

                                                       cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGGAG‑CATTAC‑T‑T‑TGCTAATCATTTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTTCGTACTATACTTATTg                                                                        >  2:229152/5‑141 (MQ=255)

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

                                                        acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATACTACTTATTgg                                                                       >  2:373499/4‑141 (MQ=255)

                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      <  2:646429/139‑1 (MQ=255)

                                                           gatagtataATAAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    <  2:40769/136‑2 (MQ=255)

                                                            atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                   >  2:764145/5‑139 (MQ=255)

                                                            actatatAATATTGGCGTGACCTTTCTTTAAGGGCTCAACTTACACGCAGTGGGGCGGTG‑TT‑GCA‑G‑G‑TGATGCTTACTTTGGTTCAGCTCATGGCTCCGGTTATTAGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTa                                                                    <  2:17241/138‑1 (MQ=255)

                                                             tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                  >  1:588513/4‑138 (MQ=255)

                                                             tgctataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGATAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                  <  1:146154/138‑4 (MQ=255)

                                                             tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa                                                                   >  1:223955/1‑141 (MQ=255)

                                                             tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTA‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTaa                                                                   >  2:19668/1‑141 (MQ=255)

                                                                tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                               >  1:47314/1‑140 (MQ=255)

                                                                tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                               >  1:240191/1‑140 (MQ=255)

                                                                 ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTGTTGGCTTAGTaaa                                                              >  2:599572/1‑139 (MQ=255)

                                                                  taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATTCTTATTGGCTTAGTaaac                                                             <  2:736297/141‑4 (MQ=255)

                                                                  taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                             >  2:292369/1‑141 (MQ=255)

                                                                   aataGTTGGCTCGCCTGTCTCCTAGTCGTCAATTTAGACGCAGAGAGTAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAc                                                             <  2:743024/141‑1 (MQ=255)

                                                                    ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                           <  1:109307/137‑1 (MQ=255)

                                                                    ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGtgagcgc                                                           >  1:356719/5‑135 (MQ=255)

                                                                    ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                            >  2:605459/1‑141 (MQ=255)

                                                                         agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                      >  1:504381/1‑138 (MQ=255)

                                                                          gattCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                     <  2:62385/138‑5 (MQ=255)

                                                                          gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                     >  1:775920/1‑137 (MQ=255)

                                                                          gagTCGCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                     <  2:46574/141‑5 (MQ=255)

                                                                            gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                   >  1:400650/1‑135 (MQ=255)

                                                                            gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                   <  1:527599/141‑7 (MQ=255)

                                                                              cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                 >  1:768522/1‑133 (MQ=255)

                                                                              cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                 >  1:314514/1‑133 (MQ=255)

                                                                               gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                <  1:342295/141‑10 (MQ=255)

                                                                               gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                >  2:185601/1‑132 (MQ=255)

                                                                                gatATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                               >  2:7948/3‑131 (MQ=255)

                                                                                ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                <  1:55901/141‑1 (MQ=255)

                                                                                ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                               >  2:65493/1‑131 (MQ=255)

                                                                                 cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaccaaatgattaa                                              >  1:173919/1‑129 (MQ=255)

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

                                                                                 cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                              <  1:71526/141‑12 (MQ=255)

                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGATGTG‑CA‑TTA‑C‑TTTGCTAATCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                            <  2:422526/141‑14 (MQ=255)

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

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

                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTACTGGCTTAGTGAACGCaacaaataattaata                                            >  2:500437/1‑128 (MQ=255)

                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatacttaata                                            >  1:539569/1‑128 (MQ=255)

                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                            >  2:217934/1‑128 (MQ=255)

                                                                                    tttctcAGGCATCAATTTAGACGAAGAGAGGATGTC‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTg                                            <  1:46378/139‑1 (MQ=255)

                                                                                    tctctcAGGCGTCAATTT‑‑ACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                           <  2:510007/141‑15 (MQ=255)

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

                                                                                        tcAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                        <  1:231264/141‑1 (MQ=255)

                                                                                          aGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGATTGCGTTTTTTTCTTATTGGCTCAGTAGACGCAATACAAAATAGGTGACatat                                      >  2:582931/1‑141 (MQ=255)

                                                                                          aGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                      <  1:361699/141‑1 (MQ=255)

                                                                                          aGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTAACTGGAGATGTGAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACACAATAGGTGACatat                                      <  1:413639/141‑1 (MQ=255)

                                                                                               tCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGCCATATAGccc                                 >  2:424931/1‑140 (MQ=255)

                                                                                               tCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCg                                 >  1:543142/1‑141 (MQ=255)

                                                                                                    ttAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATCGGTGACATATAGCCGCACCa                            <  2:89949/141‑1 (MQ=255)

                                                                                                      aGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                          >  2:419533/1‑141 (MQ=255)

                                                                                                       gACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCAGCACCAATa                         >  2:233968/1‑141 (MQ=255)

                                                                                                              agagGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcc                  >  1:25230/1‑141 (MQ=255)

                                                                                                                             aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATACCCTCACTACCGCa     >  2:20756/1‑141 (MQ=255)

                                                                                                                                  g‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACCAAATCGATGACCTATAGCAGCACACAGCAACATCCCCTCACTGCCGccata  >  2:541548/1‑137 (MQ=255)

TCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAT  >  CP000730/1929741‑1930007

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