BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,861	T→C	62.1%	intergenic (‑944/‑224)	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,861	0	T	C	62.1%	112.6 / 110.0	125	intergenic (‑944/‑224)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (25/22); new base C (47/30); total (72/53)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

AAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACC  >  CP000730/1929724‑1930001

aaGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGta‑ta                                                                                                                                               <  1:385609/141‑1 (MQ=255)

aaGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGta‑ta                                                                                                                                               <  2:334283/141‑1 (MQ=255)

  gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                                                                                            <  2:104716/141‑1 (MQ=255)

  gtggAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                                                                                            <  2:469274/141‑1 (MQ=255)

               aaTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTAttt                                                                                                                            <  2:55256/141‑1 (MQ=255)

                aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGGGGAGGTGTA‑TAA‑G‑G‑G‑GGTGCTTGtttt                                                                                                                           >  1:67697/1‑141 (MQ=255)

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

                             aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGCTAAAAGATTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATa                                                                                                              <  1:625249/141‑1 (MQ=255)

                             aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATa                                                                                                              <  1:40471/141‑1 (MQ=255)

                             aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATa                                                                                                              <  1:40282/141‑1 (MQ=255)

                              aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAg                                                                                                             >  1:104477/1‑141 (MQ=255)

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

                                      aCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCTGtc                                                                                                     >  2:546668/1‑141 (MQ=255)

                                            cAGAAGCAAGTGACTGTGCTTGTTCTTGAACCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGt                                                                                               >  1:661663/1‑141 (MQ=255)

                                            cAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGt                                                                                               >  2:655242/1‑141 (MQ=255)

                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                            >  1:621917/15‑141 (MQ=255)

                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                           >  1:443158/14‑141 (MQ=255)

                                                 tgaagtcgatatccTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                         <  1:509801/128‑1 (MQ=255)

                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                         >  1:168329/12‑141 (MQ=255)

                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTAATAAGTGGCtgt                                                                                         >  1:368623/12‑141 (MQ=255)

                                                 gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACGTAATAGCAGCAGTCATCAGTGGCtg                                                                                          >  2:189474/1‑141 (MQ=255)

                                                 gCAAGTGACTGTGCTTGTATATTTACCGATGATATAATAGTAGAGTCGCCTATCTCGCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                          >  2:490596/1‑141 (MQ=255)

                                                   aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                       <  1:442291/132‑1 (MQ=255)

                                                   aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                        >  1:291877/1‑141 (MQ=255)

                                                    agtagatattcTTCTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                      <  2:500018/131‑2 (MQ=255)

                                                    agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                      <  2:618990/133‑2 (MQ=255)

                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                     <  1:543834/134‑3 (MQ=255)

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

                                                     gcagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                     >  2:578301/8‑139 (MQ=255)

                                                       gACTGTGCTTGTATATTTGCCGATCATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCCCATCATAGCACCAGTCAGCAGTGGCTGTGCCAt                                                                                    <  1:134979/141‑1 (MQ=255)

                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgct                                                                                   >  1:476726/6‑137 (MQ=255)

                                                         atatgagtGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                 <  1:71702/134‑7 (MQ=255)

                                                         atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                 <  1:603429/138‑7 (MQ=255)

                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCGGTGGCTGTGCTGTTGCg                                                                                >  1:117086/3‑141 (MQ=255)

                                                           ttgCTTGTATAGTTGCCGATGATATAATAGTAGCGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                <  2:101255/140‑1 (MQ=255)

                                                           gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                >  1:416956/1‑141 (MQ=255)

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

                                                            tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta                                                                              >  1:32426/1‑140 (MQ=255)

                                                              cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                            <  1:557506/141‑4 (MQ=255)

                                                               ttGTATATTCACAGATGGGATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                           <  1:594485/141‑2 (MQ=255)

                                                                 gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGtattata                                                                         >  2:374439/1‑138 (MQ=255)

                                                                  tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtgttatac                                                                        >  1:187644/11‑137 (MQ=255)

                                                                  tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatat                                                                        >  2:605751/11‑137 (MQ=255)

                                                                  tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TGGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                        >  1:451585/11‑137 (MQ=255)

                                                                  tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGACGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATAATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCTGTTGCGtattatac                                                                        >  2:605744/11‑137 (MQ=255)

                                                                   atattcacagatggtATAATAATAGAGTCGCGTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTACACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                       >  1:491760/10‑136 (MQ=255)

                                                                    tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACtt                                                                      >  1:303874/9‑141 (MQ=255)

                                                                    aatTTGCCGATGATATAATAGGAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt                                                                       <  2:262162/140‑1 (MQ=255)

                                                                     tatcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                     >  2:337992/8‑141 (MQ=255)

                                                                     attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                     >  2:121783/8‑141 (MQ=255)

                                                                     atTTGCCGATGATATAATAGTAGAGTCGACTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                      <  1:463886/141‑1 (MQ=255)

                                                                         acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGCGGAGGTGCATTAC‑T‑T‑T‑GCTAGTCCTTGTCGCTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                 >  2:640275/4‑141 (MQ=255)

                                                                         acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                 <  1:252100/138‑1 (MQ=255)

                                                                          cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                <  2:591541/139‑1 (MQ=255)

                                                                          cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                >  1:9688/3‑141 (MQ=255)

                                                                           cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                <  2:40546/141‑1 (MQ=255)

                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                            >  1:183012/4‑141 (MQ=255)

                                                                               ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATTATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                           <  2:20863/139‑1 (MQ=255)

                                                                               ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                           >  2:621136/3‑141 (MQ=255)

                                                                               ggtacaatgttagAGTCGCCTATCTCCTAGGCGGCACTT‑‑GCCGAAGAGAGGAGGTGCATTCC‑T‑T‑T‑GCTAGTCATTTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                           <  2:112381/131‑1 (MQ=255)

                                                                                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGCTCACCTCAAAACACCTGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                          >  1:126114/2‑141 (MQ=255)

                                                                                  ctaataaGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                        <  1:533308/140‑3 (MQ=255)

                                                                                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                        >  1:328466/1‑139 (MQ=255)

                                                                                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                       <  1:220128/141‑4 (MQ=255)

                                                                                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAAGCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATAATTATTGGCttggtaaac                                                       >  1:91844/1‑133 (MQ=255)

                                                                                    aataatACAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                      <  2:430551/141‑1 (MQ=255)

                                                                                     ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                     >  2:336579/5‑141 (MQ=255)

                                                                                     ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTacctcgacc                                                      >  2:429641/1‑133 (MQ=255)

                                                                                      taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATAATTTTCGTTCACATCATAGCAGCAGTCATCAAAGGATGTGCTGTTGCGTCTTATACTTATTGGCTTAGTAACCGCa                                                    >  2:281208/4‑141 (MQ=255)

                                                                                         tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                 <  1:47562/141‑3 (MQ=255)

                                                                                         tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                 <  2:299261/141‑3 (MQ=255)

                                                                                         tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                 >  1:2819/1‑139 (MQ=255)

                                                                                          agagTCGCCTCTCTCTCCGGCGTCAATT‑‑GACGACGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCAGTGCTGTTGCGTCTTATACTTATTGGCTTAGTAAACGccacaa                                                >  1:190893/1‑136 (MQ=255)

                                                                                          agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                >  1:40321/1‑138 (MQ=255)

                                                                                           gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                <  1:459516/141‑1 (MQ=255)

                                                                                            agTCCCCTATCTCTCAGGCGTCAATTTATACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                               <  1:65601/141‑1 (MQ=255)

                                                                                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAACCGCaacaaata                                             >  1:213076/1‑135 (MQ=255)

                                                                                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                             >  2:430302/1‑135 (MQ=255)

                                                                                                gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                          <  1:456075/141‑10 (MQ=255)

                                                                                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCGGTCGTCATCGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                        >  1:154206/1‑130 (MQ=255)

                                                                                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                        <  2:45839/141‑12 (MQ=255)

                                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                       >  2:17659/1‑129 (MQ=255)

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

                                                                                                   tATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TAACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGGGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTAAACGCaacaaataattaat                                       <  1:603319/141‑13 (MQ=255)

                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTACTCATTTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                      <  1:349945/141‑14 (MQ=255)

                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatccatactc                                      >  2:328789/1‑128 (MQ=255)

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

                                                                                                     tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGCCTAAGTAGCAGCCATACAACATAGGTg                                      >  2:528254/1‑141 (MQ=255)

                                                                                                     tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTTGTAAACGCaacaaataattaatag                                     <  1:426032/141‑15 (MQ=255)

                                                                                                     tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattcatag                                     >  2:195512/1‑127 (MQ=255)

                                                                                                        ctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACa                                   >  2:389607/1‑141 (MQ=255)

                                                                                                              cggcAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc                             <  2:463808/138‑1 (MQ=255)

                                                                                                                    tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                       <  1:425465/141‑1 (MQ=255)

                                                                                                                    tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATGTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                       >  2:368015/1‑141 (MQ=255)

                                                                                                                          cGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCGAGGTAGATGCAATACAAAATAGGTGCCGTTTAGCCGCACCATTaga                 >  1:17030/1‑139 (MQ=255)

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

                                                                                                                           gCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTAACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGGCTAAGTAGACGGAATACATAATAGGTGACATATAGCCGCACCAATaaaa                >  2:321984/1‑141 (MQ=255)

                                                                                                                                      gtgtA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACt     >  1:231736/1‑141 (MQ=255)

                                                                                                                                         ta‑taA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc  >  1:281535/1‑141 (MQ=255)

                                                                                                                                         ta‑taA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc  >  1:620603/1‑141 (MQ=255)

AAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACC  >  CP000730/1929724‑1930001

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