BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,925	T→A	46.7%	intergenic (‑1008/‑160)	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,925	0	T	A	46.7%	44.9 / 202.5	156	intergenic (‑1008/‑160)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (45/38); new base A (40/33); total (85/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

CTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACT  >  CP000730/1929786‑1930061

 ttgttgatTTGCCGATGAGATAAGAGTAGAGTCTCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                                                                                                <  2:954764/135‑1 (MQ=255)

 ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCATGCGTCAATT‑TAGCCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                                                                                                <  2:866725/141‑1 (MQ=255)

  tGTATATTTGCCGAGGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTTAATT‑AAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                                                                               <  1:933634/141‑1 (MQ=255)

  tGTATATTCGCCGATGATATAATATTAGAGTCGCCTATCCCTCAGGCGTAAATC‑TAGACACAGAGAGAAGATTAA‑AAA‑G‑G‑C‑GA‑TGCTGAGTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCCATTGCGttttt                                                                                                                                               <  2:640740/141‑1 (MQ=255)

  tGTATATTCACAGACGGTATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                                                                              <  1:214668/141‑3 (MQ=255)

   gcatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                                                                             <  2:684907/130‑4 (MQ=255)

   gTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt                                                                                                                                              <  1:157319/141‑1 (MQ=255)

      tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                                                                                          >  2:172044/9‑135 (MQ=255)

       attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                                                                                         >  2:513247/8‑140 (MQ=255)

         tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTAtt                                                                                                                                       >  2:920493/6‑141 (MQ=255)

           caagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                     <  1:863092/138‑1 (MQ=255)

           acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATTGAAGAGGCATGGACG‑GGCAT‑‑TAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                     <  1:823026/138‑1 (MQ=255)

           acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑‑TCGACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTGTACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                     <  2:994198/138‑1 (MQ=255)

             cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                                    >  1:1050144/1‑141 (MQ=255)

             agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGCGGAGCTTCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                   <  1:1040905/140‑1 (MQ=255)

              gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                  >  2:326619/1‑140 (MQ=255)

              gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                                                                                                   >  1:542007/1‑141 (MQ=255)

               atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                                                                                                                  <  2:931290/141‑1 (MQ=255)

               atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TACACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATTGCACCAGTCAGCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                                                                                                                  <  2:399891/141‑1 (MQ=255)

                tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                                >  2:930199/4‑138 (MQ=255)

                tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑G‑T‑T‑GC‑TAATCATTTTCGTTCACATTATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                >  2:583068/4‑141 (MQ=255)

                 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                               <  2:186077/139‑1 (MQ=255)

                 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                               <  2:156225/139‑1 (MQ=255)

                 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                               >  1:1049589/3‑141 (MQ=255)

                 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                               <  1:655840/139‑1 (MQ=255)

                 ggtataAGAATCGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                               <  2:33520/139‑1 (MQ=255)

                 ggtataAAAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                               <  2:614834/139‑1 (MQ=255)

                  gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACTTAGTa                                                                                                                              >  1:59171/2‑141 (MQ=255)

                  gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                              >  2:18345/2‑141 (MQ=255)

                  gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                              <  2:539936/140‑1 (MQ=255)

                  atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCCGTGAGGGGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                               <  1:767544/141‑1 (MQ=255)

                  agataataGTAGAGTCGCCTATCTCGCAGGCGACAATT‑TAGACGCAGAGAGGAGGGGTA‑TTA‑G‑G‑T‑GA‑TGCTTATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                               <  1:395164/139‑1 (MQ=255)

                   tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                             >  1:337062/1‑140 (MQ=255)

                    ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                            >  2:144909/1‑139 (MQ=255)

                     taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                           <  2:810339/141‑4 (MQ=255)

                     aaataatAGAGCCGCCTATCGCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                           <  2:873598/140‑4 (MQ=255)

                       ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                         <  2:944018/141‑1 (MQ=255)

                       ataagagagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                         <  2:178958/136‑1 (MQ=255)

                           tagaggcctcTATCTCTCAGGTGTCAATT‑‑‑GACGAAGTGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTCTTGGCTTAGTAAACGCa‑‑aca                                                                                                                   <  1:746053/132‑3 (MQ=255)

                           tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aca                                                                                                                   >  2:737965/1‑139 (MQ=255)

                           tagAGTCGCCTATCTCTCAGGCGTCAA‑‑‑TGGACGAAGAGAGGAGGTGCA‑TTA‑G‑T‑TGGC‑TAATCATTTTCGTGCACGTCATAGCACCAGTCACCAGTGGCCGTGCTGTGGCGTATTATACTTGTTGGCATCGTAGACCCa‑‑aca                                                                                                                   >  2:548803/1‑139 (MQ=255)

                            agagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                  >  2:810943/1‑138 (MQ=255)

                             gagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                 <  1:203723/141‑5 (MQ=255)

                             gagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                 >  1:218811/1‑137 (MQ=255)

                               gTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                               >  1:40418/1‑135 (MQ=255)

                               gTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGGGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCTGTTGCGTATGATACTTATTGGCTTAGTAAAAGCa‑‑acaaata                                                                                                               >  2:970926/1‑135 (MQ=255)

                                tCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataa                                                                                                              <  1:901024/141‑8 (MQ=255)

                                 cGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTCAGTAGACGCA‑‑ATACaa‑‑aa                                                                                                            >  2:974809/1‑141 (MQ=255)

                                  gCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AAt                                                                                                           >  1:598935/1‑141 (MQ=255)

                                  gCCTATCTCTCAGGCGTCAAAT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAGTCATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATCATACTTATTGGCTTAGTACACGCa‑‑acaaataatt                                                                                                            >  1:645268/1‑132 (MQ=255)

                                   ccTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa‑‑aaga                                                                                                          >  1:830701/1‑139 (MQ=255)

                                   ccTATCTCTCAGGCGTCAATT‑‑‑GACGAAGATAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                           <  2:42778/141‑11 (MQ=255)

                                    cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGCGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATCGCACCAGTCCTCAGTGGCTGTGCTGTTTCGTATTAGACTTATTGGCTTAGTAAACGCacgaagtaattat                                                                                                            >  2:298574/1‑129 (MQ=255)

                                    cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                          >  1:699709/1‑130 (MQ=255)

                                    cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                          >  1:975657/1‑130 (MQ=255)

                                     tATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑TACTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGGATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaaa                                                                                                         >  1:438285/1‑129 (MQ=255)

                                      aTCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGt                                                                                                       >  2:846315/1‑141 (MQ=255)

                                      aTCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGt                                                                                                       >  1:1023345/1‑141 (MQ=255)

                                      aTCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                        >  1:189835/1‑128 (MQ=255)

                                        ctttcAGGCGTCAATT‑TAGACGCAGAGAGTAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGa                                                                                                     <  1:286633/138‑1 (MQ=255)

                                            cAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑a‑ta                                                                                               <  1:399780/141‑1 (MQ=255)

                                             aGGCGTCAATT‑TAGACGCGGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTACAA‑‑AATAGGTGAA‑a‑tat                                                                                              >  1:210056/1‑141 (MQ=255)

                                             aGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAAGCA‑‑ATACAA‑‑AATAGGTGAC‑a‑tat                                                                                              >  1:135758/1‑141 (MQ=255)

                                                    aaTT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTACAA‑‑AATAGGTGCC‑A‑TAT‑AGCC‑Gcc                                                                                     >  2:708616/1‑140 (MQ=255)

                                                    aaTT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCa                                                                                     >  1:608925/1‑141 (MQ=255)

                                                      tt‑tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTAGAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCAcc                                                                                   >  1:487394/1‑141 (MQ=255)

                                                      tt‑tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCCTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCAcc                                                                                   <  2:836007/141‑1 (MQ=255)

                                                          aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAAt                                                                                <  1:794324/141‑1 (MQ=255)

                                                            aCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATGGATTTTGTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATaa                                                                              >  1:771971/1‑141 (MQ=255)

                                                                   gaggagGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATccc                                                                       <  1:699398/141‑1 (MQ=255)

                                                                    aggaggTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GTACCAATAAAAATcccc                                                                      <  1:528252/141‑1 (MQ=255)

                                                                     ggaggGGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                     <  2:292634/141‑1 (MQ=255)

                                                                        ggTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                                  <  1:23388/141‑1 (MQ=255)

                                                                                aa‑G‑G‑T‑GG‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                           >  2:350452/1‑141 (MQ=255)

                                                                                 a‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTGTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                          >  2:908138/1‑141 (MQ=255)

                                                                                        ggA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATa                                                       <  2:346207/140‑1 (MQ=255)

                                                                                          gatgttgtttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAg                                                      <  1:496544/133‑1 (MQ=255)

                                                                                              cgtATGTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCCTTTTTTTCTTATTGGCCAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCCCTCTACCGCAAATAGTGAg                                                  <  2:724028/139‑1 (MQ=255)

                                                                                              cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                  >  2:468686/1‑141 (MQ=255)

                                                                                               ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTATTTTCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATACAAATCCCCTCACTACCGCAAATAGTGAgg                                                 >  2:11435/1‑141 (MQ=255)

                                                                                               ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                 <  2:188102/141‑1 (MQ=255)

                                                                                               ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTGTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                 >  1:1007127/1‑141 (MQ=255)

                                                                                                ttttttCGTTCCCATCATAGCACCAGTCATCAGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                <  1:175705/139‑1 (MQ=255)

                                                                                                tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                >  1:55913/1‑141 (MQ=255)

                                                                                                 aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                               >  2:764792/1‑141 (MQ=255)

                                                                                                  ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                              <  2:44583/141‑1 (MQ=255)

                                                                                                     tCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AGACAA‑‑GATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAGTg                                           >  2:528525/1‑141 (MQ=255)

                                                                                                           acaTCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtat                                     <  2:301165/141‑1 (MQ=255)

                                                                                                           acaTCATAGCACCAGTCATTAGTGCCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtat                                     <  1:858724/141‑1 (MQ=255)

                                                                                                           acaTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTT‑TATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGAGGGGATTgg                                          >  2:17004/1‑141 (MQ=255)

                                                                                                                atAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGAC‑TTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcat                                     >  1:415073/1‑141 (MQ=255)

                                                                                                                atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGTGGTATAAGTa                                >  2:243478/1‑141 (MQ=255)

                                                                                                                         gTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATAGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATCGTGAGGGGATCGGTGTAAAAGTAAAAACtttt                       >  2:969388/1‑139 (MQ=255)

                                                                                                                           catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                     >  2:726838/1‑141 (MQ=255)

                                                                                                                            atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCGACTAGTGAAGGGGGTGGTGGAGAGGTAAGTACtgcttat                    >  1:325314/1‑135 (MQ=255)

                                                                                                                             tcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTc                   >  1:3721/1‑141 (MQ=255)

                                                                                                                                 tGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTg               <  2:1050144/141‑1 (MQ=255)

                                                                                                                                  ggCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGATGTAGTAAATACTTATTTTCGTTGt              <  2:1023864/141‑1 (MQ=255)

                                                                                                                                         ccATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAtt       >  1:957166/1‑141 (MQ=255)

                                                                                                                                              gCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACt  >  1:462097/1‑141 (MQ=255)

CTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACT  >  CP000730/1929786‑1930061

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